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1.0 Introduction 

1.1 Statement of Work 

t A thorough understanding of the uncertainties associated with the 

modeling and testing of the Space Shuttle Main Engine (SSME) Engine will 
i 1 greatly aid decisions concerning hardware performance and future 
f development efforts. The goals of this research effort are delineated in the 
the Statement of Work, reproduced below: 


The goal of this effort is to enhance the rocket engine steady-state performance 
computer models through the incorporation of uncertainty analysis concepts. Analytical 
tools and analysis techniques are being investigated to better support performance 
analysis requirements. These requirements include assessing vehicle/engine feed system 
interface flow characteristics, engine hardware design changes, evaluating engine 
hardware performance, predicting engine hardware operation, and supporting failure 
investigations. 

A major shortcoming within the current SSME power balance modeling scheme is that 
experimental data and the fundamental physical relationship are treated as absolutes. 
Both experimental data and the property requirements contain various sources or errors. 
The primary sources of error within the instrumentation system are calibration errors, 
si gnal processing and localized effects. The primary sources or errors within the 
fundamental relationships are uncertainties in physical approximations and fluid 
property computational predictions are forces to agree with the data at instrumented 
locations, often at the expense of physical consistency. This situation degrades the 
capability of the analytical tools and thus, reduces the amount of confidence in the 
results generated. 

A test data integration strategy was developed based upon evaluating test data with 
respect to basic fluid conservation principles (mass, energy, and momentum 
relationships). This strategy systematically transforms uncertain experimental data into 
a physically self consistent set of data. This is accomplished by forcing the minimum 
adjustment required in engine pressures, temperatures, and flowrates necessary to satisfy 
prescribed uncertainty constraints. 

This strategy incorporate uncertainty requirements explicitly. The overall success of the 
test data integration strategy is a function of determining these required uncertainty 
estimates. 

Another major shortcoming within the current modeling scheme is that the current 
performance models do not present uncertainty estimates with predictions. A general 
framework for modular uncertainty estimates with predictions. A general framework for 
modular rocket engine performance models do not present uncertainty estimates with 
predictions. A general framework for modular rocket engine performance prediction 
program is currently being developed that will integrate physical principles, rigorous 
mathematical formalism, component and system level test data, and theory-observation 
reconciliation. This development effort will allow for simple implementation of 
uncertainty estimates associated with physical relationships. Incorporation of these 
estimates within the rocket engine performance model will support two crucial 
functions. First, These uncertainty estimates will represent a confidence band associated 
with each prediction. Secondly, these estimates will provide a measure of the success of 
the performance model. 
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The research required to implement these uncertainty analysis concepts will be 
conducted within the SSME engine 3001 test program which is currently being 
conducted on the technology test bed (TTB) test facility. Engine 3001 provides a 
significantly larger number of propellant property measurements as compared to 
standard SSME modeling strategies, and ultimately, improve the use of test data in 
generating performance predictions. Phase 1 involves applying uncertainty analysis 
techniques to obtain estimates for both the bias and the precision uncertainties of engine 
measurements. Phase 2 involves incorporating uncertainty analysis techniques to 
estimate uncertainties associated with model computation. A detail description of the 
specific tasks to support these phase are described below: 

Phase 1 - Perform uncertainty analysis for engine 3001 test measurements. 

1. Examine TTB instrumentation systems 

2. Evaluate data from previous TTB testing 

3 . Identify all significant sources of errors 

4. Estimate both precision and bias uncertainties for the TTB test measurements 

5. Evaluate the use of these measurement uncertainty estimates by the PRM for 
supporting TTB test analysis 

Phase 2 - Perform uncertainty analysis on the physical relationships within the 
Performance Reconciliation Model (PRM) 

1. Identify assumptions and physical analysis on the physical approximations 
made by representing a real physical process as a mathematical model. 

2. Determine methods for quantifying the influence of such assumptions and 
physical approximations. 

3 . Estimate the modeling uncertainties. 

4. Evaluate the use of these uncertainty estimates within the PRM computations 
for supporting TTB test analysis. 






1.2 Report Overview 

This report will describe the determination of uncertainties in the 
modeling and testing of the Space Shuttle Main Engine test program at the 
Technology Test Bed facility at Marshall Space Flight Center. Section 2 will 
present a summary of the uncertainty analysis methodology used and discuss 
the specific applications to the TTB SSME test program. Section 3 will 
discuss the application of the uncertainty analysis to the test program and 
the results obtained . Section 4 presents the results of the analysis of the 
SSME modeling effort from an uncertainty analysis point of view. The 
appendices at the end of the report contain a significant amount of 
information relative to the analysis, including discussions of venturi 
flowmeter data reduction and uncertainty propagation, bias uncertainty 
documentation, technical papers published, the computer code generated to 
determine the venturi uncertainties, and the venturi data and results used in 
the analysis. 
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2.0 Uncertainty Analysis 

The vise and application of uncertainty analysis in engineering has 
evolved considerably since Kline and McClintock's classic paper 1 in 1953. 
Developments in the field have been especially rapid and significant over the 
past decade, with the methods formulated by Abemethy and co-workers 2 that 
were incorporated into ANSI/ASME Standards in 1984 3 and 1986 4 being 
superseded by a more rigorous approach 5 . Publication in late 1993 by the 
International Organization for Standardization (ISO) of the Guide to the 
Expression of Uncertainty in Measurement? in the name of ISO and six other 
international organizations has, in everything but name only, established a new 
international experimental uncertainty standard. 

The approach in the ISO Guide deals with "Type A" and "Type B" 
categories of uncertainties, not the more traditional engineering categories of 
bias and precision uncertainties, and is of sufficient complexity that its 
application in normal engineering practice is unlikely. This issue has been 
addressed by AGARD Working Group 15 on Quality Assessment for Wind 
Tunnel Testing and by the Standards Subcommittee of the AIAA Ground Test 
Technical Committee. The documents 6 - 7 produced by these groups present and 
discuss the additional assumptions necessary to achieve a less complex "large 
sample" methodology that is consistent with the ISO Guide, that is applicable to 
the vast majority of engineering testing (including most single-sample tests), and 
that retains the use of the traditional engineering concepts of bias and precision 
uncertainties. (The chapters on uncertainty methodology in the AGARD 6 and 
AIAA 7 documents were authored by the Principal Investigator of this research 
program.) 

2.1 Overview 

The word accuracy is generally used to indicate the relative closeness of 
agreement between an exp erimentally-determined value of a quantity and its 
true value. Error (5) is the difference between the experimentally-determined 
value and the truth, thus as error decreases accuracy is said to increase. Only in 
rare instances is the true value of a quantity known. Thus, one is forced to 

1 Kline, S. J., and McClinlock, F. A., "Describing Uncertainties in Single-Sample Experiments," Mechanical 
Engineering, VoL 75, 1953. 

2 Abemethy, R. B., Benedict, R. P., and Dowdell, R. B., "ASME Measurement Uncertainty,” J. Fluids 
Engineering, VoL 107, 1985. 

3 American National Standards Institute/American Society of Mechanical Engineers, Measurement Uncertainty 
for Fluid Flow in Closed Conduits, MFC-2M-1983, ASME, 1984 

4 American National Standards Institute/American Society of Mechanical Engineers, Measurement Uncertainty, 
PTC 19.1-1985 Part 1, ASME, 1986. 

5 International Organization for Standardization, Guide to the Expression of Uncertainty in Measurement, ISO, 
ISBN 92-67-10188-9, 1993 

6 Quality Assessment for Wind Tunnel Testing, AGARD-AR-304, 1994. 

7 American Institute of Aeronautics and Astronautics, Assessment of Wind Tunnel Data Uncertainty, AIAA 
Standard S-071, 1995. 
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estimate error, and that estimate is called an uncertainty, U. Uncertainty 
estimates are made at some confidence level -- a 95% confidence estimate, for 
example, means that the true value of the quantity is ejected to be within the 
±U interval about the experimentally-determined value 95 times out of 100. 


As shown in Figure 1(a), total error 5 can be considered to be composed of 
two components: a precision (random) component e and a bias (systematic) 
component p. An error is classified as precision if it contributes to the scatter of 
the data; otherwise, it is a bias error. It is assumed that corrections have been 
made for all systematic errors whose values are lcnown. The remaining' bias 

. a! 1*1 1 __ 1 TaJ 4-2 



Suppose that we are making a number of measurements of the value of a 
variable X that is absolutely steady. The k and k+1 measurements are shown in 
Figure 1(a). Since the bias is a fixed error, it is the same for each measurement. 
However, the precision error will have a different value for each measurement. 
It then follows that the total error in each measurement will be different, since 
the total error is the sum of the bias error and precision error in a measurement. 

If we continued to take measurements as previously described until we 
had a sample of N readings, more than likely as N approached infinity the data 
would behave as shown in Figure 1(b). The bias error would be given by the 
difference between the mean (average) value p of the N readings and the true 
value of X, whereas the precision errors would cause the frequency of occurrence 
of the readings to be distributed about the mean value. 

As an estimator of J3, a bias limit B is defined?. A 95% confidence estimate 
is interpreted as the experimenter being 95% confident that the true value of the 
bias error, if known, would fall within ±B. A useful approach to estimating the 
magnitude of a bias error is to assume that the bias error for a given case is a 
single realization drawn from some statistical parent distribution of possible bias 
errors. For example, suppose a thermistor manufacturer specifies that 95% of 
samples of a given model are within ±1.0 C of a reference resistance-temperature 
(R-T) calibration curve supplied with the thermistors. One might assume that 
the bias errors (the differences between the actual, but unknown, R-T curves of 
the various thermistors and the reference curve) belong to a Gaussian parent 
distribution with a standard deviation b=0.5 C. Then the interval defined by ±B 
= ± 2b = ±1.0 C would indude about 95% of the possible bias errors that could be 
realized from the parent distribution. (The bias limit is sometimes called the 
"systematic uncertainty*'.) 

As an estimator of the magnitude of the precision errors (the width of the 
distribution of readings in Figure 1(b)), a precision limit P is defined 8 . A 95% 
confidence estimate of P is interpreted to mean that the ±P interval about a 
single reading of Xi should cover p 95 times out of 100. (The precision limit is 
sometimes called the "precision uncertainty".) 

In nearly all experiments, the measured values of different variables are 
combined using a data reduction equation (DRE) to form some desired result. A 
good example is the experimental determination of mass flow rate using a 
venturi meter as discussed in Appendix II of this report. Functionally, the mass 
flow rate is given as 


8 Coleman, H. W., and Steele, W. G., Experimentation and Uncertainty Analysis for Engineers, Wiley, New 
York, 1989. 


7 


We = W e (P,T,XP,d,D,a,C D ) ( 1 ) 

One can envision that errors in the values of the variables on the right hand side 
of Eq. (1) will cause errors in the experimental result We. 

A more general representation of a data reduction equation is 

r = r(X„ X 2 , Xj) (2) 

where r is the experimental result deter mi ned from J measured variables X. 
Each of the measured variables contains bias errors and precision errors. These 
errors in the measured values then propagate through the data reduction 
equation, thereby generating the bias and precision errors in the experimental 
result, r. 

If the "large sample assumption" is made 67 then the 95% confidence 
expression for Ur becomes 

, , O) 

ui - SefBf + 2 2 j e,e,B * 

i=l i=l k=H+l 

+5>f.pf + j 6i6 k p* 

i=l i~l k=i+l 


where 


6i = 


dr 

dXi 


(4) 


and where the 95% confidence precision limit for a variable X is estimated as 


Pi = 2 Si N't 10 


and the sample standard deviation is calculated using 


St = 


1 

N-l 


t KXth - Xif 


(5) 

( 6 ) 


where the mean value is defined as 


Xt 


N 


N 


_k=l 


(7) 


and Pit is the 95% confidence estimator of the covariance of the precision errors 
in X and Xk, and B a is the 95% confidence estimator of the covariance of the bias 
errors in X and X. 


If we define the bias limit (systematic uncertainty) of the result as 


8 


( 8 ) 


Bt = 2>?5? + I diOkB* 

i-1 i~l k=i+l 

and the precision limit (precision uncertainty) of the result as 

Pr = 'ZtiP? + 2 2 Z Aft/** (9) 

i=7 i=l Jc=i+1 


then Eq. (3) can be written as 


Ur = tf+Pt ( 10 ) 

and Eqs. (8) and (9) can be viewed as propagation equations for the bias limits 
and precision limits, respectively. 

2.2 Determining Precision Limi ts 

Single Test. When the result is determined from a single test -- that is, 
at a given test condition the result is determined once using Eq. (2) 

r = r(Xi, X 2 , Xj) (2) 

and when the Xi's are considered single measurements, then Eq. (9) is used to 
find the precision limit of the result. This situation is often encountered in large 
scale engineering tests in which measurements of the variables are made at a 
given set point over a period that is small compared to the periods of the factors 
causing variability in the experiment. A proper precision limit (one indicative of 
the dispersion of the variable over several cycles of the factors causing its 
variation) cannot be calculated from readings taken over such a small time 
interval. For such data, the measurements) of a variable Xi should be 
considered a single reading -- whether the value of X is the average of 10, 10 3 or 
10 6 readings taken during the short measurement time. In such a test, the value 
for the precision limit to be associated with a single reading would have to be 
based on previous information about that measurement obtained over the 
appropriate time interval 9 . If previous readings of a variable over an 
appropriate interval are not available, then the experimenter must estimate a 
value for Pi using the best information available at that time 6 - 7 . 

For single tests in which some of the variables (X2 and X3, for instance) 
can be determined as averages from multiple readings over an appropriate time 
period but the other variables cannot be, then 

r = r(Xi, ~X 2 , X 3 Xj) (11) 


9 Steele, W. G., Taylor, R.P., Burrell, R. E., and Coleman, H. W., “The Use of Data from Previous 
Experience to Estimate the Precision Uncertainty of Small Sample Experiments,” AIAA Journal, Vol. 31, 
No. 10, 1993. 
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and Eq. (9) is used to find the precision limit of the result as follows. For the 
variables that are single readings, the Pi's are the precision limits determined 
from previous information or estimated from the best available information. For 
the averaged variables when N 2 and N 3 are equal to or greater than 10, P 2 and 
P 3 should be taken as precision limits of means, (2S2)/(N2) U2 and (2S3)/(N3) 1/2 , 
with the S's calculated using Eq. (6). When N 2 and N 3 are less than 10, it is the 
authors' recommendation that the precision limits used in Eq. (9) for the 
averaged variables be taken as (P 2 )/(N 2 ) 1/2 and (Ps)/^) 172 , where P 2 and P 3 are 
determined from previous information, as is done for the single reading 
variables. 


For tests in which multiple readings of all of the variables can be obtained 
over an appropriate period, the following method is recommended. 


Multiple Tests. If a test is performed so that M multiple sets of 
measurements (Xi, X 2 , ... Xj)k at the same test condition are obtained, then M 
results can be determined using Eq. (2) and an average result r can be 
determined using 


r 


1 M 
= — Y 
M ^ 


rt 


k=l 


( 12 ) 


If the M sets of measurements were obtained over an appropriate time period, 
the precision limit that should be associated with a single result would be 

Pr = tSr (13) 


where t is determined with M-l degrees of freedom and is taken as 2 for M>10 
and & is the standard deviation of the sample of M results 


Sr 


i 

M-l 


M _ 

S (r k - rf 

k=l 


\l/2 


(14) 


The precision limit that should be associated with the average result is given by 


= _A_ 

r 4m 


(15) 


with Pr given by Eq. (13). Using the large sample assumption, the uncertainty 
that should be associated with a single result would be 


and with an average result r 


U 2 r =B?+(2S r ) 2 



(16) 

(17) 


with Br given by Eq. (8). 
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Correlated Precision Uncertainties. The Pik terms in Eq. (3) take into 
account the possibility of precision errors in different variables being correlated. 
These terms have traditionally been neglected 1 - 3 - 4 - 5 - 7 , although precision errors in 
different variables caused by the same uncontrolled factor(s) are certainly 
possible and can have a substantial impact on the value of the precision limit 10 . 
In such cases, one would need to acquire sufficient data to allow a valid 
statistical estimate of the precision covariance terms to be made if using Eq. (3). 
Note, however, that the multiple tests approach using Eq. (14) implicitly 
includes the correlated error effect -- a definite advantage when multiple sets of 
measurements over an appropriate time period are available. 


2.3 Estimating Bias Limits 

Bias Limits of Individual Variables. When attempting to estimate 
the bias limits Bi of the individual variables in Eq. (8), one might separate the 
bias errors which influence the measurement of a variable into different 
categories: calibration errors, data acquisition errors, data reduction errors, test 
technique errors, etc. Within each category, there may be several elemental 
sources of bias. For instance, if for the Jth variable, Xj, there are M elemental 
bias errors identified as significant and whose bias limits are estimated as (Bj)i, 
(Bj> 2 , ..., (Bj)m, then the bias limit for the measurement of Xj is calculated as the 
root-sum-square (RSS) combination of the elemental limits 


Bj — 


t(Bj A 


(18) 


The elemental bias limits, (Bi)k, must be estimated for each variable X 
using the best information one has available at the time. In the design phase of 
an experimental program, manufacturer's specifications, analytical estimates 
and previous experience will typically provide the basis for most of the estimates. 
As the experimental program progresses, equipment is assembled, and 
calibrations are conducted, these estimates can be updated using the additional 
information gained about the accuracy of the calibration standards, errors 
associated with the calibration process and curvefit procedures, and perhaps 
analytical estimates of installation errors. 


As Moffat 11 suggests, there can be additional conceptual bias errors 
resulting from not measuring the variable whose symbol appears in the data 
reduction equation. An example would be a point temperature measurement 
interpreted to be indicative of a cross-section averaged temperature, but there 
may be a cross-sectional variation of temperture, which may or may not have a 


10 Hudson, S. T., Bordelon, W., and Coleman, H. W., “Effect of Correlated Precision Errors on the 
Uncertainty of a Subsonic Venturi Calibration,” AIAA-95-0797, 1995. 

11 Moffat, R. J., “Describing the Uncertainties in Experimental Results,” Experimental Thermal and Fluid 
Science, Vol. 1, 1988. 
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predictable profile, causing the “average” value to be different than the point 
value. Hence, the inclusion of an elemental bias term for the conceptual error 
would be appropriate. 

Correlated Bias Limits. Correlated bias limits are those that are not 
independent of each other, typically a result of different measured variables 
sharing some identical elemental error sources. It is not unusual for the 
uncertainties in the results of experimental programs to be influenced by the 
effects of correlated bias errors in the measurements of several of the variables. 
A typical example occurs when different variables are measured using the same 
transducer, such as multiple pressures sequentially ported to and measured 
with the same transducer or temperatures at different positions in a flow 
measured with a single probe that is traversed across the flow field. Obviously, 
the bias errors in the variables measured with the same transducer are not 
independent of one another. Another common example occurs when different 
variables are measured using different transducers all of which have been 
calibrated against the same standard, a situation typical of the electronically 
scanned pressure (ESP) measurement systems in wide use in aerospace test 
facilities. In such a case, at least a part of the bias error arising from the 
calibration procedure will be the same for each transducer, and thus some of the 
elemental bias error contributions in the measurements of the variables will be 
correlated. 

The Bik terms in Eq. (8) must be approximated -- there is in general no 
way to obtain the data with which to make a statistical estimate of the 
covariance of the bias errors in X and the bias errors in Xj. The approximation 
of such terms was considered in detail in Ref. 12, where it was shown that the 
approach that consistently gives the most satisfactory approximation for the 
correlated bias limits was 


B* - t, (Bt) a (Bk) a (19) 

a=l 

where L is the number of elemental systematic error sources that are common 
for measurements of variables X and X. 

If, for example, 


r = r(Xi, Xi) ( 20 ) 

and it is possible for portions of the bias limits Bi and B 2 to arise from the same 
source(s), then Eq. (8) gives 


12 Brown, K. K., Coleman, H. W., Steele, W. G., and Taylor, R. P., "Evaluation of Correlated Bias 
Approximations in Experimental Uncertainty Analysis," AIAA 94-0772, 1994. 
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Bf = e 2 ,B 2 1 +02B 2 2 +2e I e 2 Bi2 ( 21 ) 

For a case in which the measurements of Xi and X 2 are each influenced by 4 
elemental error sources and sources 2 and 3 are the same for both Xi and X 2 , Eq. 
(18) gives 


= (B,f 1 + (B,f 2 +(B,f 3 + (B,f 4 (22) 


and 


B I = (B 2 f 1 + (B 2 f 2 +(B 2 f 3 + (B 2 f 4 (23) 


while Eq. (19) gives 


B, 2 = (B,) 2 (B 2 ) 2 +(B,) 3 (B 2 ) 3 


(24) 


2.4 Application of Uncertainty Analysis to TTB Testing 

The focus in this effort is to identify the uncertainty that should be 
associated with a measured variable such as temperature or pressure or with a 
determined result such as flowrate that is calculated using a number of 
measured variables. The uncertainty given by Eq. (16) -- that associated with a 
single result -- is the appropriate uncertainty to use when data from a single 
TTB test are compared with the output of a predictive model. 


Desired Variable 

Environment 
(Installation & conceptual 
biases; unsteadiness) 

Sensor 

(Calibration biases) 


Data Acquisition System, DAS 
— ■ — ^ ■- I (Calibration biases) 

Measured Value of 
V ariable 

Figure 2.2 Schematic of TTB Measurement Process from an Error Sources Viewpoint 
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Figure 2.2 shows a schematic of the viewpoint used in identifying error 
sources that contribute to the overall uncertainty. The desired variable is taken 
to be the one with which a model output will be compared -- a cross-section 
averaged temperature, for example, that would usually be referred to as "the 
temperature" of the flow at a particular location in the engine. If the sensor 
responds to temperature at a point, then an installation or conceptual bias exists 
due to the sensor not actually responding to the desired variable (the average 
temperature). This is an elemental source that must be included in Eq. (18), and 
it is potentially one of the dominant elemental sources in temperature and 
pressure data in TTB testing. The estimation of these sources is a task in the 
follow-on effort to the work documented in this report. The traditional 
"measurement uncertainty" sources are shown as biases in the sensor calibration 
and biases in the calibration of the data acquisition system (DAS). Additionally, 
the effect of unsteadiness in, and due to, the operating environment must be 
considered since the sensor calibrations and DAS pre-test calibration checks are 
not done with the engine operating, and the unsteadiness certainly can have an 
effect on the final system output - the measured value of the variable. 

Choice of the appropriate precision limit to use with TTB data needs to be 
carefully done. A precision limit determined using a standard deviation from a 
time slice during one test gives information about the steadiness of the "steady 
state" at that operating condition during that particular test, but includes no 
effects of the test-to-test variation of the variable at that operating condition. As 
discussed in Section 3, the authors believe that computing a standard deviation 
of a variable or result from multiple tests, all of which were at the same 
operating condition, gives the appropriate precision estimate for use in 
discussing the uncertainty in a measured TTB variable. It is also the 
appropriate precision limit to consider when comparing the results from one test 
to results from another test in an effort to determine if a change in component, 
for instance, had any discernible effect on the value of the result. 

Detailed discussion of the uncertainty estimates associated with TTB 
flowrate measurements is given in Section 3 of this report. 

2.5 Application of Uncertainty Analysis in SSME Modeling 

When comparing output of a model with experimental data, the 
imcertainties that should be associated with the model predictions must be 
considered for proper conclusions to be drawn. In the past, most of the work 
reported in this area has simply considered the sensitivity of the model output to 
uncertainties in the input data. This obviously does not include any 
uncertainties in the model itself and thus is not a satisfactory approach. In this 
research effort, we have divided the sources that cause uncertainty in the model 
output into three categories: (1) imcertainties due to assumptions and 
approximations in the model, (2) uncertainties due to the incorporation of 
previous experimental data into the model, and (3) imcertainties due to the 
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numerical solution algorithm. Consideration of the third category is not within 
the scope of this program. 

The first category, uncertainties due to assumptions and approximations 
in the model, does not include the installation and/or conceptual bias source 
shown in Figure 2.2 and discussed above since that uncertainty is associated 
with the measured value. Consider the temperature at a particular position in 
the flow. The uncertainty associated with the measured value of the 
temperature includes the effect of malting a point measurement but desiring a 
cross-sectional averaged value. The inability of the model to calculate a correct 
average temperature at a particular location because the one-dimensional flow 
approximation has been made results in an uncertainty in the predicted 
temperature. (Stated another way, if the model predicts the correct average 
temperature at a particular location, then the one-dimensional flow 
approximation has caused no uncertainty in the model output.) 

The uncertainties due to the incorporation of previous experimental data 
in the model arise when material property data is used, when valve resistance 
characteristics are used, when pump maps are used, etc. These are all instances 
in which previous experimental data has been used by replacing the data with 
curvefits. The original data contained uncertainties, but the curvefit equations 
used in the predictive models have been treated as the "truth" in most previous 
considerations of uncertainty in model outputs. Adding further complication, 
there is no accepted way of estimating the influence of systematic uncertainties 
on the uncertainty that should be associated with a regression. This aspect has 
been investigated in this program and is discussed in Section 4. An AIAA paper 
reporting the progress of this effort is included as Appendix IV. 
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3.0 Results of Application of Uncertainty Analysis to TTB Testing 

An investigation to determine the experimental uncertainties 

associated with the test measurements from the SSME Engine 300 1 installed 
in the Technology Test Bed facility was conducted. This investigation 
consisted of reviewing existing documents, discussions with NASA personnel, 
review of other technical literature, and new analyses. Since the 
thermodynamic performance analysis of the SSME was the motivation 
behind this contractual effort, the pressure, differential pressure, 
temperature, and mass flow rate measurements were the focus of the 
investigation. Initial discussions concluded to initially focus on the 
determination of uncertainty in the flowrate measurements, with particular 
emphasis on the venturi flowmeter determinations. This section will discuss 
the information obtained upon which the assessment of individual 
uncertainty source estimates were made. 

3.1 Measurement System 

The Technology Test Bed test measurement system is described here 
as all components between the phenomena being measured and the final 
computer data file in engineering units, including the sensors, transducers, 
data acquisition systems, and data reduction routines. A previous study by 
Sverdrup Technology 13 studied the MFSC test facilities to assess the 
uncertainty in the measurement systems and to assess if any significant 
discrepancies existed between test areas. That report identified and 
quantified the uncertainties between the sensor and the final data file. 
However, several important points about this study need to be made. First, 
this study did not assess the uncertainty in a given measurement during the 
engine test, therefore any additional uncertainty due to the operating 
environment was not assessed. Secondly, no installation or conceptual 
uncertainties were considered. Finally, the study was performed in 1992 and 
prior to the adoption of new standards for the assessment of uncertainties 5 , 
and some of the specific procedures used to combine uncertainty sources were 
not in accordance with the current standard. 

The specific aspects of the TTB measurement system were reviewed 
with the TTB data acquisition personnel. This review showed the procedures 
and techniques being used are self-consistent, with a pre-test procedure 
conducted which recalibrates the data acquisition system prior to each test. 
This ensures that measurement system drift and gross errors do not go 
undetected. 

3.2 Analysis of Previous TTB Test Data 


13 Fish, James E., NASA/MSFC Test Area Measurement System Uncertainty Study, Sverdrup 
Technology, Inc., MSFC Group, Report No. 335-002-92, October 1992. 
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The review of previous test data was cumbersome due to the format of 
the data being incompatible with commercially available software and the 
initial difficulties in utilizing MSFC resources for data analysis. To achieve a 
set of data to review which could be defined as the “same” hardware, tests 
TTB039 through TTB 051 were chosen. These tests were conducted with 
Engine 3001 with the large throat combustion chamber and a consistent set 
of other hardware. Another difficulty in the data analysis to assess 
uncertainties was the lack of repetition in the test profiles, with each test 
being conducted for the analysis of specific performance aspects. The data 
was reviewed in the full sample-rate format (25 or 50 samples/sec) and at the 
reduced sample-rate format (1 sample/sec). Because of inconsistencies in the 
way test data is stored in the computer systems, all of the measurements for 
these tests are not readily accessible. For example, most of the venturi mass 
flowrate calculations are not available from the computer systems being 
used. This created problems in assessing the precision uncertainties for the 
venturi flowrate measurements, particularly trying to assess test-to-test 
precision uncertainty behavior. To alleviate this situation a new computer 
program was developed by the COTR to access the test data directly and 
compute the mass flowrates. This program was then modified by the 
researchers to calculate the uncertainties associated with the given test data. 
A discussion of this new software tool is given in Section 3.5. 

3.3 Determination Of Mass Flow Rate Uncertainties 

The mass flow rate uncertainties for the venturi flowmeters were 
determined using the methodology previously discussed in Section 2. The 
data reduction equation is the equation for mass flow rate presented in 
Appendix III and the expression for the uncertainty in the mass flow rate 
and the necessary partial derivatives are presented in Appendix IV. This 
expression shows the density of the fluid within the square root, however the 
density of the fluid cannot be explicitly measured and must be determined 
within the data reduction by using the measured pressure and temperature 
and some equation of state. 

The bias limits used in the uncertainty propagation equations were 
estimated based upon the information gathered from the available 
documentation, discussions with TTB personnel, and engineering experience. 
The values used in the calculation of the flowrate uncertainties are shown in 
Appendix V. 

Potentially significant bias uncertainties to be addressed in the TTB 
measurements are the conceptual bias uncertainties discussed in Section 2. 
The conceptual bias uncertainties in the temperature and pressure 
measurements are particularly important in this effort because of the 
interest in comparing the experimental results with the analytical 
predictions. In many of the SSME measurements the flowfield is highly 
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complex due to the sharp turns and bends, valves, pump and turbine inlets 
and discharges, and other complicating factors. These factors accentuate the 
difference between the physical quantity at the sensor and the quantity for 
which the measurement is desired, typically an average value at a cross- 
section. Assessment of these bias uncertainties has not been accomplished. 
These assessments will require extensive review of the measurement, the 
sensor and its installation, the thermodynamic and fluid dynamic flowfield, 
and their interaction. This detailed analysis of the measurements will take 
place during the next contract period and based upon a prioritized list 
developed with the COTR. 

The precision limits for the mass flowrate uncertainties are dependent 
upon the question being asked, or rather, what is purpose for the 
information. Precision limits can be calculated in many different ways, but 
the interpretation of the precision limit and the use of it depends upon data 
used to calculate it. The variables which must be considered for the precision 
limit calculation include: 
engine number 

specific engine component configuration 
engine test(s) 
power level 
test profile 

specific engine adjustments 
time slice within the test 

data sample rate (data points used for standard deviation calculation) 

The precision limits for the flowrate uncertainties were based upon 
review of the flowrate data for the chosen time slice. Precision limits were 
estimated in two primary ways. First, one was based upon the full sample- 
rate data within each test. The second precision limit was estimated based 
upon averaging the full sample-rate data over a given time slice, for a given 
test, at a chosen power level to provide a single data point for that test 
condition and using similar points from other tests to estimate a precision 
limit. 


The uncertainty in the thermophysical property data must be included 
in the uncertainty analysis since there was some experimental uncertainty in 
the original experiments upon which the property tables were developed. If 
the property data is represented by a curve-fit of the experimental data an 
uncertainty associated with the data and with the curve-fit must be included. 
The venturi flowrate data reduction utilizes the thermophysical property 
routine GASP 14 . The GASP documentation shows the uncertainty in density 


14 Hendricks, Robert C., Baron Anne K., and Peller, Ddiko C., “GASP - A Computer Code for 
Calculating the Thermodynamic and Transport Properties for Ten Fluids: Parahydrogen, Helium, Neon, 
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for liquid hydrogen and liquid oxygen to be within 0.25%. The GASP 
program is based upon National Bureau of Standards data for the 
thermodynamic properties of hydrogen and oxygen. 

3.4 Calibration of Mass Flowmeters, Determination of Discharge 
Coefficient, Cd 

The accurate determination of the mass flow rate for cryogenic rocket 
propellants is difficult because of the special problems presented. These 
problems arise because of the extremely cold temperatures, the property data 
uncertainty, the difference in density and viscosity with respect to the 
calibration fluid, and limitations on the calibration procedures because of 
safety considerations Each of these can introduce an uncertainty in the mass 
flow rate determination. 

The accurate deter min ation of a mass flow rate depends upon the 
ability to trace the output from a given device back to a standard certified by 
the National Institute of Standards and Technology (NIST) or some other 
respected standard. The basic problem with cryogen flowmeters is that there 
are very limited facilities which can produce an accurate standard using the 
actual cryogenic fluid. The facilities which do exist are limited to the 
calibration of liquid nitrogen flows, or have limited capacity. A few facilities 
exist which can perform calibration with the appropriate cryogen, but the 
reference standard in these systems use a meter which was calibrated with 
water and then adjusted for the particular cryogenic fluid. This procedure 
relies on the concept of dynamic similarity, matching the Reynolds number of 
the test fluid with the Reynolds number of the calibration fluid, and making 
corrections for the dimensional changes due to thermal contraction and other 
miscellaneous effects. A literature survey indicates that very little work has 
taken place to assess the accuracy of this procedure, primarily due to the cost 
and complexity involved with developing the necessary experiments with 
liquid oxygen and liquid hydrogen. Most of the documented work describes 
the calibration and use of turbine-type flowmeters. Since the facility 
flowmeters are turbine meters this is an important aspect to study. One of 
these studies conducted with liquid oxygen which used liquid nitrogen as the 
calibration fluid indicated that there is a difference on about the same order 
of magnitude as the uncertainty in the physical property data, about 0.25% 15 . 
This study also used liquid argon and the researcher concluded that turbine 
flowmeters have a dependence on the fluid properties. However, that study 
was too limited to make a firm conclusion. Since a property dependence was 


Methane, Nitrogen, Carbon Monoxide, Oxygen, Fluorine, Argon, and Carbon Dioxide,” NASA TN D- 
7808, February 1975. 

15 Brennan, J.A. LaBrecque, J.F., and Kneebone, C.H., ’’Progess Report on Cryogenic Flowmetering at 
the National Bureau of Standards,” Instrumentation in the Cryogenic Industry, Volume 1, Proceedings of 
the first Biennial Symposium on Cryogenic Instrumentation, October 1 1-14, 1976, Held in Conjunction 
with the 3 1st Annual ISA Conference and Exhibit 
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noted between the cryogenic argon and nitrogen in comparison with the 
cryogenic oxygen, one would expect a difference between water at room 
temperature and cryogenic oxygen. Another study 16 investigated the use of 
calibrating the liquid hydrogen turbine flowmeters with high-pressure 
nitrogen because: “Comparisons have been made between water and liquid 
hydrogen calibrations. Results of these comparisons indicated that for 
inaccuracies less than 1 percent at full-scale flow, water calibrations are 
inadequate in predicting the meter constant for liquid-hydrogen flow in 
turbine-type flowmeters.” 

Very little literature is available which discusses attempts to assess 
the uncertainty due to differences in the calibration fluids for differential 
pressure producer mass flowmeters, such as venturis and orifices, 
particularly for the operating Reynolds numbers of the SSME venturis. 
However, it seems obvious from a review of the fluid and thermodynamic 
processes going on in the venturi that a non-negligible difference exists. The 
description of the venturi flowmeter data reduction equation shown in 
Appendix HI and the Rocketdyne venturi calibration report 17 show that 
si gnific ant assumptions and approximations are made to achieve the data 
reduction equations. Two of the fundamental assumptions are that the flow 
through the venturi is adiabatic and one-dimensional, which is not true 
during engine operation. Another significant factor is that none of the 
venturis are installed with the recommended length of upstream straight 
duct. Rocketdyne attempted to account for this problem by calibrating the 
venturis with the “as-installed” ducts connected, however they could not 
reproduce the other flow characteristics, such as flow swirl, turbulence, 
oscillations, etc. 

The Rocketdyne venturi calibration report develops a polynomial 
curve-fit to obtain the value of the discharge coefficients where the curve-fit 
is based upon the Reynolds number of the flow through the venturi. They 
recommended using an iterative method to determine the discharge 
coefficient to use to determine flowrate through the venturi. For all of the 
venturi flowmeters calibrated with water in Rocketdyne’s facility dynamic 
similarity could not be obtained, that is, the flowmeters could not be 
calibrated at the Reynolds number expected during an engine test. In fact, 
most of the calibrations were performed at a fraction of the operating 
Reynolds number, from 1.7% to 10%, and the curve-fits were extrapolated out 
to the operating condition. 


16 Szaniszlo, Andrew J., and Krause, Lloyd N., “Simulation of Liquid-Hydrogen Turbine-Type Flowmeter 
Calibrations using High-Pressure Gas,” NASA TN D_3773. 

17 Lepore, Frank A., Rocketdyne Division, Space Shuttle Main Engine No. 3001, Technology Test Bed, 
Differential Flowmeters Calibration Final Report, Contract No. NAS8-27980, March 1980. 
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The extrapolation of the calibration curve-fits and the use of a 
different fluid for the calibration introduce an uncertainty which must be 
associated with the determination of the mass flowrate. In the absence of 
perfor min g detailed experiments to assess these uncertainties, an estimate 
must be made based upon the existing information. The information upon 
which to base these estimates are the calibration data figures in the 
Rocketdyne report and the literature review previously discussed. A review 
of the calibration data figures shows that a curve was fit through the data, 
but for which a different curve would produce a substantially different result 
at the operating conditions. Figure 3.1 can be used as an example. The 
curve generated by Rocketdyne is used in the data reduction, however it can 
be observed that a significantly different calibration curve would have been 
obtained if a higher order curvefit was used, or if the last few data points 
were more heavily weighted. 

VENTURI h /M CALIBRATION DATA 


TTBE/ITEM 296 



REYNOLDS NUMBER - Re(dx) 
li TEST »1 + TEST «2 

v TEST o3 CXPNL.CURVEFIT 

Figure 3.1 Venturi #296 Calibration Curve 


3.5 Venturi Flowmeter Uncertainty Computer Program 

A computer program was developed to calculate the uncertainty in the 
mass flowrate determinations. This program was based upon a program 
initially developed by the COTR to calculate the flowrates for the venturi 
flowmeters and expanded to perform the uncertainty calculations. This 
program resides within the EADS 10 (Engineering Acquisition and Data 
System) computer system and assesses the raw, full sample-rate data for all 
of the engine measurements . The program uses the measured pressure, 
differential pressure, and temperature for each venturi, combined with the 
specific dimensional data and other constants (such as discharge coefficient), 
averages the data for the chosen time slice and calculates the mass flowrate. 
The bias limits and precision limi ts for the uncertainty sources for each 
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venturi are stored in the program as constants and propagated through data 
reduction equation using the uncertainty propagation equation presented in 
Appendix IV. The only correlated bias uncertainty included in the 
uncertainty propagation equation is the uncertainty in the venturi 
dimensions. All other possible bias correlations were considered negligible, 
particularly with respect to the magnitude of the discharge coefficient bias 
uncertainty. The partial derivatives are determined numerically using a 
finite difference technique, often referred to as a jitter routine. The data 
used for the bias limits and precision limits and a short description of how 
the estimates were obtained is shown in Appendix V. 

The TTB venturi data reduction utihV.es the methodology 
recommended by Rocketdyne which includes an iterative technique to 
determine the value of the discharge coefficient. The purpose behind the 
iterative technique is to provide a value for the discharge coefficient from the 
polynomial curvefit generated during the calibration. However, the value of 
the discharge coefficient does not change much within the operating range, so 
that using a constant value of Cd is more appropriate. The iterative 
technique tends to encourage a higher level of confidence than is warranted, 
based upon the previously discussed uncertainties in the discharge 
coefficients. The program developed to determine the uncertainty in the 
mass flowrates does not use this iterative technique and provides essentially 
the same values for the calculated flowrate. 

This program will allow the uncertainties in the venturi flowrates to 
be determined from any TTB test and allows the recalculation of the 
uncertainties when the bias and precision limits for the individual venturis 
are updated. 

3.6 Flowrate Uncertainty Results 

Tests TTB039 through TTB050 were chosen as the basic data to 
conduct the uncertainty analysis upon. These tests were conducted with the 
large-throat combustion chamber and represented engine configurations 
upon which comparisons to the analytical model were desired. A time slice at 
each of three power levels (100% RPL, 104% RPL, and 109% RPL) was 
chosen, and the full-sample rate (50 samples/sec) data was averaged over the 
time slice to obtain a single data point for each PID at each power level. The 
time slices used at each power level are shown in Table 3.1. 

The averaged data was used in the venturi mass flowrate and mass 
flowrate uncertainty program to obtain the venturi flowrate and the bias 
limit for each power level for each test. It was observed that the bias limit, as 
a percent of the flowrate, was constant for all tests and all power levels. The 
venturi mass flowrate, the bias limit, and the percent bias uncertainty 
calculated for each of the tests considered at each power level and for each 
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Test# 

100% RPL 

i 104% RPL 

109% RPL 

TTB039 

80-85 sec 

95-100 sec 

195-200 sec 

TTB040 

60-65 sec 

80-85 sec 

120-125 sec 

TTB041 

80-85 sec 

95-100 sec 

195-200 sec 

TTB042 

80-85 sec 

95-100 sec 

195-200 sec 

TTB043 

50-55 sec 

75-80 sec 

34-39 sec ! 

TTB044 

45-50 sec 

100-105 sec 

34-39 sec i 

TTB045 

100-105 sec 

150-155 sec 

185-190 sec i 

TTB047 

10-15 sec 

35-40 sec 

50-55 sec 

TTB048 

115-120 sec 

85-90 sec 

50-55 sec 

TIB049 

145-150 sec 

50-55 sec 

25-30 sec 

TTB050 

16-21 sec 

25-30 sec 

34-38 sec 


Table 3.1 Time Slices used for Precision Limit Determinations 
* TTB046 not included because test profile did not include a steady-state 
time slice at chosen power levels. 

venturi are shown in the tables in Appendix VII. The bias limits, standard 
deviations, and uncertainties for the venturis are shown in Table 3.2, for 
100% rated power level. Appendix VII contains the mass flowrate data for 
the tests from which the standard deviations were calculated as well as the 
calculated bias limits for each test. It is observed that the bias limit is 
primarily a function of the bias uncertainty in the discharge coefficient. 

The preferred method of determining the precision limit for a result is 
to calculate the standard deviation of the results (mass flowrates in this case) 
instead of propagating the precision limits of the individual variables. This 
allows for any correlated precision error sources to be inherently accounted 
for. Thus a set of tests which represents the nominal engine configuration is 
desired. The series of tests between TTB039 and TTB050 was chosen 
because the engine configuration was based upon the large-throat 
combustion chamber and design changes. However, during this series of 
tests a Taguchi test matrix was conducted and engine components were 
varied simultaneously. The changes between the standard Rocketdyne 
designed high pressure oxidizer turbopump (HPOTP) and the Pratt and 
Whitney HPOTP, and changes in the high pressure fuel turbopump create 
problems with data analysis. Since these engine components have different 
performance the grouping is not appropriate. If the components changed had 
all been “line-replacement” units, where each unit is expected to provide the 
same performance as the replaced unit the grouping would be meaningful. 
The specific engine component configuration for each test is shown in Figure 
A-VI.l in Appendix VI. Review of this table showed that four tests were 
conducted with Rocketdyne HPOTP’s, with the other major components being 
the same or being line replacement units, thus this series of four was 
grouped. Other specific engine test variable changes are potentially 
responsible for some of the precision uncertainty, but the lack of test profile 
repetition in the TTB program makes specific distinctions diffi cult. 


23 



















































A standard deviation for the flowrates for the ten tests, as well as the 
four-test group, designated as RD HPOTP, was obtained, as shown Appendix 
VII. When these standard deviations are combined with the calculated bias 
limits using the “large-sample” approximation the results for 100% RPL are 
shown in Table 3.2, below. 


PID 


w 

(lb/sec) 

Bw 

(%) 

Sw(%) 
10 Tests 

S w (%) 

RD HPOTP 

U w (%) 

10 Tests 

U w (%) 

RD HPOTP 

8801 

LPFT Inlet 

28.0 

2.1 

2.3 

1.9 

5.1 

4.3 

8818 

CCV Inlet 

73.8 

3.0 

3.4 

0.6 

7.4 

3.2 

8815 

Nozzle Clnt #1 

13.8 

2.1 

5.5 

4.1 

11.2 

8.5 

8816 

Nozzle Clnt #2 

13.5 

2.1 

4.7 

1.2 

9.6 

3.2 

8817 

Nozzle Clnt #3 

13.6 

2.1 

18.5 

7.5 

37.1 

15.1 

8802 

LPOT Inlet 

188. 

2.1 

2.4 

1.8 

5.2 

4.2 

8819 


947. 

2.1 

1.2 

0.8 

3.2 

2.6 

8804 

OPBLOX 

26.1 

2.1 

8.7 

1.6 

17.5 

3.8 

8805 

OPBLH 2 

35.7 

2.0 

3.6 

1.5 

7.5 

2.5 

8810 

FPBLOX 

61.3 

2.1 

4.6 

3.9 

9.4 

8.1 


Table 3.2 - Venturi Mass Flowrates and Flowrate Uncertainties @100% RPL 


PID #8817, Nozzle Coolant #3 is shown to have a very large precision 
uncertainty at all three power levels. Futher review of the test data shows 
that the large standard deviation is caused by exceptionally low 
measurements during tests TTB043 and TTB045. This could be an actual 
result due to a specific hardware characteristic, a random phenomena, or due 
to an instrumentation malfunction. From the information available a specific 
cause cannot be identified and thus the data point cannot be eliminated. 

It is interesting to note from Tables A-VII.l - A-VII.6 in Appendix VII 
that the standard deviation of the total of the 10 venturi flowmeters is much 
lower than the standard deviation of any of the individual flowmeters, 0.8% 
at 100% RPL, for the 10 tests and 0.5% for the four tests. Since the primary 
engine control point is the main combustion chamber pressure, and since the 
chamber pressure is directly proportional to the total flowrate (for a set 
mixture ratio), this result is not suprising. It also indicates the overall 
performance variation from test-to-test is small, however the variation within 
the engine changes from test-to-test as the engine re-balances due to specific 
hardware changes. 

The input data for the venturi flowrate calculations are shown in 
Table A-VI.2 and show that the pressure and temperature measurements 
have a small standard deviation and the differential pressure measurements 
have a large test-to-test variation. This verifies the opinion of the TTB 
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instrumentation engineers that the differential pressure measurements were 
not as accurate as the other measurements. 
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4.0 Application of Uncertainty Analysis in SSME Modeling 

Assessing the uncertainties in the output of the SSME model is 

important to making decisions based upon those results. Uncertainties in 
the model output can be separated into three general categories, 
uncertainties due to modeling assumptions and approximations, 
uncertainties due to the numerical solution, and uncertainties due to using 
previous experimental information. 

The primary methodology identified in the literature for assessing 
uncertainties in analytical models is a perturbation technique, perturbing 
the model inputs and checking the resultant output. This technique cannot 
adequately account for all of the identified modeling uncertainties and 
another methodology is needed. 

The model of the SSME being used in this effort is a model based upon 
the rocket engine modeling platform developed by Pratt and Whitney, 
ROCETS (Rocket Engine Transient Simulations). This model provides a 
more structured format wherein the individual engine components are 
modeled in individual modules and the modules are solved simultaneously to 
provide the performance prediction. 

4. 1 Modeling Assumptions and Approximations 

When a mathematical or engineering model of a physical system is 
developed certain assumptions and approximations about the system are 
made to simplify the system to one for which mathematical expressions can 
describe. By making these simplifications an error is introduced and the 
model cannot exactly describe the physical system. 

Some of the primary assumptions and approximations made within 
the SSME model include: 

• 1-dimensional 

• fully developed 

• steady-state 

• adiabatic 

• ideal gas 

The primary problem with attempting to assess uncertainties with 
respect to these assumptions and approximations is that if the uncertainty to 
associate with a particular assumption or approximation can be estimated, 
then the model could be improved to include this estimate instead of trying to 
estimate the uncertainty. For example, it was determined by a researcher 
working on a complementary effort that a turbine exit temperature was being 
predicted using an ideal gas, constant specific heat approximation, which for 
the specific temperature range of interest was a poor approximation. Instead 
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of trying to estimate an uncertainty to associate with that approximation the 
model is being altered to include a better thermodynamic description of the 
process. Hence, extensive effort within this research program was not 
warranted in this area. 

4.2 Numerical Solution Uncertainties 

When the system of equations is solved numerically the exact solution 
will not be obtained. The error in the numerical solution is a combination of 
the round-off error specific to the computer system and the truncation error 
in the numerical solution scheme. The magnitude of the numerical solution 
uncertainties are assumed to be negligible with respect to the other 
uncertainties. Hence, extensive effort was not warranted in this area under 
this contract. 

4.3 Uncertainties from using previous experimental information 

In all of the component modules, information from previous testing is 

used. For example, the model of the liquid oxygen flow through a given duct 
or through a given valve is based upon its component testing, which provides 
an equation for the resistance through the duct as function of the flowrate. 
This test information is often reduced to the form of a line or curve. 

The value of the discharge coefficient used in the test data reduction or 
the polynomial constants in the thermodynamic property routines are 
examples of using previously obtained uncertain test information in a model. 

4.4 Linear Regression Uncertainty 

The methodology to assess the uncertainty in the coefficients of a 
linear regression were developed as part of this effort. The uncertainty 
analysis methodology discussed in Section 2 was applied to the expressions 
for the regression coefficients to develop the technique. The details of this 
methodology were presented at a recent conference 18 . This technical paper is 
attached in Appendix IV. The work presented in this paper demonstrates 
that this technique provides a method to determine the uncertainty in linear 
regression coefficients, and includes the effect of some correlated bias 
uncertainties. This technique is the first methodology which considers bias 
uncertainties and correlated bias uncertainties in determining regression 
coefficient uncertainties. 

The methodology is fully developed at this time for linear regressions, 
but extention to the more general regression forms used in most models 
requires additional effort. 


18 Brown, Kendall K., Coleman, Hugh W., and Steele, Glenn W., “Estimating Uncertainty Intervals for 
Linear Regression,” AIAA Paper 95-0796, 33rd Aerospace Sciences Meeting and Exhibit, January 9-12, 
1995, Reno, NV. 
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5.0 Summary 

An uncertainty analysis of the SSME test program at the Technology 
Test Bed facility has shown that non-negligible uncertainties exist in the 
determination of the mass flowrates, ranging from 2% to up to 10%, and even 
larger in some cases. The bias uncertainty is dominated by the bias limit in 
the experimentally determined venturi discharge coefficient. The bias limi ts 
for the venturis can be refined as more, and more specific, information about 
the instrumentation and the installation bias uncertainties is obtained. 

The uncertainty intervals provided in this report can be used for 
comparison to performance prediction models; however the engine 
configuration being modeled must be representative of the engine 
configuration for which the precision uncertainty was determined. 

A primary source of uncertainty in the analytical modeling of the 
SSME is from the use of previous experimental information. This 
information is usually utilized in the form of regressions or curvefits. Initial 
research efforts to determine a methodology to properly account for the 
uncertainty in the regression coefficients has shown to be promising, and 
work is continuing to develop the methodology. 
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Appendix I 

Venturi Flowmeter Data Reduction 

The data reduction methodology for the determination of mass flow rate from differential 
pressure venturi flowmeters is reproduced from Fluid Meters: Their Theory and 
Application by ASME 19 and the Rocketdyne Final Report [ref 17], 

The mass flow rate through a differential pressure venturi mass flow meter is designated as 
We, with units of [lbm/sec] and the data reduction equation is: 

W, = KYf m C w Jt& 


where K is a unit conversion factor 


= 1 \ 2g e 

4 V12 2 


0.5250204 


Y a is an expansion factor 



where k is ratio of specific heats, Cp/Cy and x is the pressure ratio, defined as 


x = 



F a is a thermal expansion factor, defined as 

F a = l + 2a(r-528) 


G\y is the calibration coefficient for the particular venturi, defined as 


, £Cp 


where J3 is the ratio of throat diameter to the inlet diameter, f3=d/D, and Cj) is the 
discharge coefficient. The discharge coefficient must be determined experimentally, 
Rocketdyne developed curvefit equations for Cq as a function of inlet Reynolds number. 
For all the venturis except No. 247 the curvefit equation is 


19 Fluid Meters: Their Theory and Application, 6th Ed., The American Society of Mechanical Engineers, 
New York, NY, 1971. 
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c d = a- 


B. 


l Re, J 


where the Reynolds number is determined from 


Re, = 


4W e 

%\xd 


(R e d)t ls the transition Reynolds number and the curvefit for venturi item No. 247 is 


C D = A + £(Re d ) + C(Re d ) 2 +D( Re d ) 3 +E(Re d ) 4 

Further information and details on the venturi calibration is found in the Rocketdyne 
Flowmeter Final Report, reference [17], 
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Appendix II 

Mass Flow Rate Uncertainty Determination 


When the equations in Appendix I are combined for incompressible fluids (Y a = 1) the 
equation for the mass flow rate can be written as 

, pc , 

W e = (0.525)(1 + 2 <x(T- 528)) ■ ' 4 JpAP 

fW 

and performing a general uncertainty analysis the equation for the uncertainty in the mass 
flow rate can be written as 


uir = e 2 a ui + e 2 T u 2 T + e 2 u 2 + e 2 u 2 D + e 2 u c 




ul 


and the partial derivatives are written as 


el 
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el = (0.525)(l + 2a(T - 528)) 
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0* = (0. 525)(l + 2a(T - 528)) ■ 


d 2 C T 


'-(U 



=(0.525)(l+2cr(T-528)) 




Since the density of the fluid is not a measured variable, but must be determined as a 
function of the measured pressure and temperature, it is replaced in the data reduction 
equation with its functional relationship. 


W e = (0525)(1 + 2<z(T- 528)) 



MpjTw 


Using this data reduction equation and applying the uncertainty propagation equation in a 
slightly different form we obtain the expression for the uncertainty as 



KdbPj 




dW, 
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where it is assumed that the uncertainty in the determination of the density from the 
property subroutine is absorbed within the bias limits estimated for the temperature and 
pressure. From review of the thermophysical property subroutine used in the data 
reduction, the uncertainty in the density given by the routine is approximately 0.25% at 
95% confidence. Note that the precision limit of the flowrate is combined directly and not 
propagated with the partial derivatives, as are the bias limits. This method provides a 
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better estimate of the precision limit and automatically accounts for any correlated 
precision error behavior within the measured variables. 


If the bracketed term is multiplied and divided by W 2 , terms such as 

i f lV 

w 2 U cj cj 


are obtained. This term for Cd is important since it indicates that the square of the bias 
limit for the discharge coefficient is linear with respect to the square of the uncertainty in 
the flowrate. This indicates that the uncertainty in the discharge coefficient is an important 
parameter. Similarly the other terms of this nature can be obtained, (except for the 
pressure and temperature used in the density function) but are more complicated and are 
not as easily interpreted. 
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ABSTRACT 

The Space Shuttle Main Engine (SSME) tests conducted 
at Marshall Space Flight Center’s (MSFC) Technology 
Test Bed (TTB) facility are being used to assess the 
performance of new SSME hardware components. The 
experimental data from the engine tests is analyzed and 
compared with predictions from the engine s performance 
prediction model. A research effort is in progress to 
quantify experimental uncertainties and analytical 
uncertainties within the SSME performance prediction 
model. This paper will discuss some of the unique 
problems encountered in quantifying these uncertainties. 

INTRODUCTION 

The Space Shuttle Main Engine (SSME) tests conducted 
at Marshall Space Flight Center’s (MSFC) Technology 
Test Bed (TTB) facility are being used to assess the 
performance of new SSME hardware components. These 
new components are being tested in the SSME Engine 
3001 program because it provides a significantly higher 
number of instrumentation points than a normal flight 
engine. The data is being used to guide the development 
of new hardware for SSME flight engines as well as 
future advanced liquid hydrogen/liquid oxygen engines. 
The test data is also used with the SSME Power Balance 
Model (PBM), the engine’s steady-state performance 
prediction model, to assess the predictions the model 
provides and in some cases to make alterations to the 
model based upon the test data. The fact that the test 
data contains uncertainties has not been taken into 
account, and no assessment has been made of the 
uncertainty in the model's predictions due to the 


assumptions, approximations, and uncertain data used in 
the development and solution of the model. The goal of 
the current effort is to enhance rocket engine steady-state 
performance models and rocket engine performance 
analysis through the incorporation of uncertainty analysis 
techniques. This paper discusses this on-going effort and 
some of the unique problems being investigated. 

TECHNOLOGY TEST BED FACILITY 

The Technology Test Bed (TTB) facility "was established 
to validate new propulsion technology hardware advances 
for large liquid oxygen and liquid hydrogen rocket 
engines.” (NASA, 1991) The need for a facility in which 
detailed knowledge of the internal environment of large 
rocket engines could be obtained was recognized as being 
necessary for the technology development required for 
advanced hydrogen/oxygen engines. Planning for the TTB 
began in 1982 to use a highly instrumented Space Shuttle 
Main Engine and the test stand at MSFC which was 
originally constructed for the testing of the first stage of 
the Saturn V, a five engine cluster of F-l RP-1/L0 2 rocket 
engines. The facility modifications were completed in 
1984 and the first tests were conducted in 1989 with a 
SSME with standard flight instrumentation(NASA, 19911. 
Testing of the highly instrumented SSME, identified as 
Engine 3001, began in December of 1990 and has 
continued since then. The technology development items 
which have been tested in the TTB include the following 
(NASA, 1991): 

* Large throat combustion chamber 

* Combustion stability studies 

* Engine start transient modifications 
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* Engine shutdown evaluation 

* Test turnaround time reduction 

* Continuous fuel system purge 

* LH 2 umbilical leak tests 

* Advanced fuel turbopump 

* Advanced oxidizer turbopump 

Engine 3001 provides "over 750 special measurements, 
including flowrates, pressures, temperatures, and strains” 
(McConnaughey, et al., 1992). 

SSME PERFORMANCE ANALYSIS 

The major requirements of the SSME performance 
analysis include assessing vehicle/engine feed system 
interface flow characteristics, engine hardware design 
changes, evaluating engine hardware performance, 
predicting engine hardware operation, and supporting 
failure investigations. In much of this analysis the test 
data is compared to the predictions provided by the SSME 
performance prediction model, commonly referred to as 
the Power Balance Model (PBM). The PBM was 
originally formulated during the initial development of the 
SSME and has been continuously modified and updated. 

In its current form, the PBM utilizes information from 
previous tests and property data for the propellants to 
develop a model prediction. A problem with this approach 
is that both the experimental data and the physical 
relationships in the model are treated as absolutes. The 
data from previous tests and the property data contain 
uncertainties. The physical relationships in the model are 
only approximations (1-D, etc) of reality and therefore 
have their own uncertainties. Previously, experimental 
data has been treated as absolute and computational 
predictions were forced to agree with the data at 
instrumented locations, often at the expense of physical 
consistency (conservation of energy, for example). 

The highly empirical nature of the Power Balance Model 
provides significant areas for uncertainties to be 
introduced; however, there is no mechanism in the model 
to propagate those uncertainties and provide an 
uncertainty associated with the models prediction. Also, 
because of its empirical nature the model's capability to 
provide meaningful predictions for engines which contain 
significant new hardware is in question. 

In an effort prior to the one discussed in this paper, a 
test data/model reconciliation strategy was developed 
based upon evaluating test data with respect to basic fluid 
conservation principles (mass, energy, and momentum 
relationships). This strategy systematically transformed 
uncertain experimental data into a physically self 
consistent set of data. This was accomplished by forcing 
the minimum adjustment required in engine pressures, 
temperatures, and flowrates necessary to satisfy 
prescribed uncertainty constraints. This strategy 


explicitly incorporated uncertainty estimates for 
experimental data as well as the physical properties, but 
did not account for model uncertainties. 

In a research effort parallel to the one discussed in this 
paper, a new performance model is being developed. This 
model is based upon a general rocket engine performance 
analysis computational tool developed by Pratt & Whitney 
known as ROCETS (ROCket Engine Transient 
Simulation). While the model was developed with the 
capability of performing transient analyses, starts, 
throttling, and shutdowns, it is also useful as a steady 
state prediction model. ROCETS is a modular platform in 
which modules are put together, along with the specific 
hardware and thermodynamic cycle characteristics, to 
generate a computational model for most types of rocket 
engines. This is advantageous for the SSME Engine 3001 
program since the specific hardware characteristics for a 
particular new component can be easily changed within 
the ROCETS model. These types of changes would be 
extremely difficult to accomplish within the PBM. The 
development of the reconciliation strategy for 
incorporation with the ROCETS based model (similar to 
that discussed in the preceding paragraph) is a concurrent 
research effort (Santi, 1994). 

CONSIDERING UNCERTAINTY IN TESTING AND 
MODELING 

Since both experimental results and "predictions” from 
a model contain uncertainties, one needs to consider how 
these uncertainties should be taken into account when 
evaluating a comparison of test results and model output. 
The general case is shown schematically in Figure 1. A 
physical system is studied both experimentally and 
analytically. In most situations, the analytical model is 
"calibrated” using experimental data from the physical 
system being studied or from a similar system - friction 
factor vs. Reynolds number correlations, for example. The 
model also uses property data, which originally was found 
experimentally and contains uncertainties. 

The methodology for estimating the experimental 
uncertainties has received considerable attention, and the 
most current approaches are presented and discussed in 
Coleman and Steele (1989), ISO (1993) and AGARD 
(1994). To the authors' knowledge, no equivalent 
methodology exists for estimating the uncertainty that 
should be associated with the "predictions" of a model. 

In the analytical modeling approach the physical system 
is studied from the aspect of what physical laws can be 
applied and which laws of physics will help answer the 
posed questions. In a fluid system the conservation laws 
(continuity, momentum, energy) are important; however, 
these laws are generally not applied in an exact sense in 
a model. To obtain a solution which will answer the 


2 


BROWN, ET. AL. 


posed question assumptions and approximations to the 
physical laws are typically made— for example, making a 
one-dimensional approximation for a three-dimensional 
physical process or assuming a structural material 
behaves linearly with a constant modulus of elasticity. 

There are many other assumptions and approximations 
which are commonly used to reduce general equations to 
solvable forms, including assumptions such as homogeous, 
isotropic, steady, incompressible, isothermal, adiabatic, 
isentropic, constant acceleration, rigid body, constant 
rotation, etc. It should be realized that these assumptions 
and approximations introduce an uncertainty into the 
output of the model. As mentioned above, often in the 
analytical approach results from an experiment or 
physical property data must be introduced along with the 
assumptions and approximations, thus introducing more 
elements of uncertainty. With the assumptions and 
approximations reducing the general physical laws to a 
solvable form and the necessary other data included, the 
model is complete and a mathematical solution can be 
found, that is, a prediction obtained from the model. This 
solution contains some uncertainty, but there is not any 
recognized methodology to quantify the uncertainty. In 
computational approaches, such as computational fluid 
dynamics, there are methods to estimate the error 
involved in the numerical solution of the differential 
equations. But little attention has been paid to how the 
other assumptions and approximations affect the degree 
of goodness of the final solution. 

This can be taken further by asking how should the two 
answers, the experimental solution and the analytical 
solution, be compared. It is extremely unlikely that the 
two answers will be the same and since each has an 
uncertainty interval associated with it what can be said 
about how the two uncertainty intervals interact? Do the 
intervals overlap? If so, what does it mean and what 
conclusions can or should be drawn? If the uncertainty 
intervals do not overlap and a significant difference exists 
between the analytical solution and the experimental 
solution, is there justification for adjusting the analytical 
solution with a "fudge-factor”? 

UNCERTAINTY IN SSME TESTING 

The aspect of the Engine 3001 test program that is of 
primary interest in' this effort is the engines 
thermodynamic performance. As such, the variables of 
primary interest are the mass flowrates, pressures, and 
temperatures in various points of the engine. When 
Engine 3001 was developed 9 flowmeters, 15 temperature 
transducers, and 12 pressure transducers were added. 
Both differential pressure venturi flowmeters and orifice 
plate flowmeters are used for the flowrate determinations, 
depending on the location and type of flow to be 


measured. The temperature probes used are mainly 
thermocouples with a few RTD's in special locations, and 
the pressure transducers include absolute transducers, 
differential pressure transducers, and some high 
frequency pressure transducers. 

The operating environment in a Space Shuttle Main 
Engine is very harsh, much harsher than that in most 
systems in which an uncertainty analysis has been 
applied. The combustion chamber operates at greater 
than 3000 psia, the high pressure fuel turbopump 
pressurizes the hydrogen fuel to over 5500 psia, and the 
oxygen is pressurized to over 3800 psia. At these 
pressures both of the propellants are in the supercritical 
region while at cryogenic temperatures. The size 
constraints necessary for the SSMEs to fit in the Space 
Shuttle dictate that the ducting be very compact. The 
extreme nature of the SSME environment creates some 
unique problems in instrumentation and measurement 
and in the interpretation of those measurements. 

One of the major tasks in this research effort is to assess 
the major contributors to the uncertainty in the 
measurements taken in the TTB Engine 3001 program. 
Some of the more significant uncertainty sources that 
have been identified and are being investigated are 
discussed below. These include the venturi discharge 
coefficient uncertainty and conceptual uncertainties. 

In the analysis of the SSME performance the system 
flowrates are key parameters. The determination of the 
uncertainty in the discharge coefficient for the differential 
pressure venturi flowmeters appears to be a major 
contributor to the flowrate uncertainty. The SSME s are 
manufactured by Rocketdyne Corporation, and Engine 
3001 was a specially built engine with modifications to the 
piping and ducting to install the venturi nozzles, orifice 
plates, and pressure and temperature probe taps. 
Because one of purposes of the Engine 3001 program was 
to characterize the internal flow environment in the 
SSME, Engine 3001 needed to be as similar to flight 
engines as possible, but with the added instrumentation. 
This created a number of problems for the design and 
installation of the venturi flowmeters, the most critical of 
these being encroachment on the required length of 
straight duct upstream of the nozzle entrances. The 
venturis were also designed with consideration of 
minimizing the pressure drop to help ensure that the 
cryogenic fluids would not be prone to cavitation and to 
minimize the added flow resistance. 

In an attempt to account for some of these issues, 
Rocketdyne designed the calibration system to calibrate 
the venturis with the associated connecting ducting. 
They then performed calibrations on the venturis in their 
facility with ambient temperature water as the 
calibration fluid since a cryogenic calibration facility 
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capable of handling similar flowrates was not, and is not, 
available. The dynamic similarity between the calibration 
fluid and the operating fluid conditions was relied upon to 
provide appropriate values for the discharge coefficients 
for the test venturis. In most of the engine venturis the 
maximum Reynolds number tested during the calibration 
process was less than the Reynolds numbers expected 
during engine testing, in some cases by over an order of 
magnitude. The orifice plates were calibrated by Colorado 
Engineering Experiment Station, Inc. and were calibrated 
with the associated ducting from the engine, but were 
calibrated with ambient air as the calibration fluid and 
good Reynolds number correlations were able to be 
obtained.. 

The potential difference between calibrating the venturis 
and orifices with water or air and testing with 
supercritical cryogenic fluids is an effect we believe might 
be a significant source of uncertainty. From the lack of 
published literature and from comments from individuals 
working in the field, at this time it appears this effect has 
not been fully addressed. 

Another important area of uncertainty to be investigated 
is the conceptual uncertainty associated with the meaning 
of a particular measurement. Moffat (1988) introduced 
the concept of conceptual bias error as the error due to 
the difference in what the instrument or probe was 
reading and the use of that measurement in the data 
reduction routine. For example, if the flowrate of a fluid 
in a duct is desired, and the velocity of that fluid is 
measured at a single point in the flow field, there will a 
conceptual bias error due to the probe measuring the 
velocity at a single point instead of measuring the average 
velocity of the bulk flow. The complicated and complex 
nature of the flows in the piping and ducting of the 
engine, due to sharp bends, valves, pump exits, and the 
like, assure that the readings of the pressure tap or 
thermocouple located at the inner surface of the pipe are 
not representative of the bulk flow in the duct. 

Difficulties arise in quantifying the value of the 
conceptual bias errors. If one is able to determine an 
estimate for the conceptual error, that estimate could be 
used to adjust the measured value and an uncertainty due 
to the correction included. The error can be estimated by 
either analytical methods or by additional 
experimentation. In the Engine 3001 program additional 
experimentation for the purpose of assessing conceptual 
bias errors is prohibitive, with respect to the overall 
testing goals and the competition for resources. Thus the 
process of making analytical estimates for the conceptual 
bias errors is being undertaken. In some cases, CFD 
modeling of the flows in particular ducts has been 
previously completed for other purposes and that 
information will be used in formulating estimates. 


The uncertainty introduced with the use of equations of 
state to represent the physical properties- of the 
propellants also needs to be addressed, particularly since 
much of the system operates much above the critical 
point. Experimental physical property data for hydrogen 
and oxygen at these conditions is limited and that 
available has a non-negligible uncertainty, about 0.5% for 
liquid hydrogen. 

UNCERTAINTY IN SSME MODELING 

As previously mentioned, the ROCETS SSME model is 
based upon each hardware component having its own 
module in the program which contains the mathematical 
model for that component. Each of the modules is a one- 
dimensional thermodynamic and fluid dynamic model for 
the particular piece of hardware. The ROCETS program 
contains the equations for the hardware components 
typically found in rocket engines and inputs for the 
specific hardware characteristics for the engine being 
modeled adjust the general equations to complete the 
model. The modules then form a set of equations which 
have to be solved. The ROCETS platform is a more 
structured platform which is more conducive to the 
development of the modeling uncertainties. 

In order to understand why the modelling of the SSME 
is such a difficult task, a description of the engine is 
presented. The Space Shuttle Main Engine operates in a 
staged combustion cycle. The staged combustion cycle is 
the highest performance rocket engine cycle available and 
accordingly the most complicated. The maximum amount 
of available energy is extracted from the propellants by 
using the propellants to provide the power to operate the 
turbopumps and then ensuring that all of the propellants 
enter the main combustion chamber. The combustion 
processes occurring in the oxidizer preburner and the fuel 
prebumer occur at approximately 5500 psia, and the main 
combustion chamber operates at 3160 psia, at 100% rated 
power. 

In trying to assess the uncertainties which are 
introduced in the modelling process a comparison between 
the real operation of the engine and the operation of the 
engine as it is being modelled is necessary. For example, 
consider the flow of liquid oxygen through the main 
oxidizer valve. In the real engine the flow is three- 
dimensional, non-uniform, highly turbulent, and heat 
transfer is occurring. In the modeling process it is 
assumed that the axial velocity is the only important 
velocity and the flow is completely uniform, the effect of 
the turbulence can be characterized by the Reynolds 
number or the velocity in conjunction with the valve 
characteristics, the flow is adiabatic and the process is 
isenthalpic. A method to estimate the uncertainty 
propagated through the model due to each of these 
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assumptions and approximations must be formulated. 

The modeling of the combustion processes is another 
uncertainty contributor. It models the combustion 
processes with an adiabatic flame temperature calculation 
using pressure corrected ideal gas assumptions. It does 
not account for real thermodynamic effects such as 
dissociation and non-equilibrium, which are potentially 
significant effects at the high chamber pressures involved. 
What is the uncertainty propagated through the model 
due to these assumptions and approximations? Again, a 
method to quantify these uncertainties must be developed 
and their influence on the overall result assessed. 

After the conditions for which the model is to be solved 
and reasonable initial guesses are input, the ROCETS 
model uses a combination of a Newtonian and a Broyden 
numerical solution scheme to obtain a solution for the 
conditions throughout the engine. This solution scheme 
is relatively robust - it converges to a solution quickly and 
reliably. However, the solution scheme only provides a 
result within some preassigned tolerance, thus 
introducing an uncertainty due to the mathematical 
solution process and necessitating a methodology to assess 
this uncertainty. 

SUMMARY 

The ultimate goal of this research effort is to improve 
the development of rocket propulsion systems by providing 
information about the uncertainty of the experimental 
data and analytical predictions upon which decisions are 
based. Assessment of the experimental uncertainties in 
the SSME test program is presenting some interesting 
challenges due to the severe environment of the engine. 
Assessment of the uncertainties in the analytical 
predictions is even more challenging because there is not 
an accepted methodology to follow. 

As of August 1994, several tasks are being pursued in 
support of the two research efforts described above. The 
progress in support of quantifying experimental 
uncertainties include the gathering of necessary 
information about the SSME and test program, a review 
of the instrumentation systems, identification of some 
significant uncertainty sources, and assessment of the 
uncertainties in the venturi flow measurements. 
Determination of the final values for flow rate, pressure, 
and temperature uncertainties will be completed by the 
end of CY 94. The progress in support of quantifying 
analytical uncertainties within the SSME performance 
prediction model includes identification of the key 
questions to be investigated. Further work on the 
modeling uncertainties will be pursued upon completion 
of the ROCETS SSME test data reductions model (Santi, 
1994). 
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FIGURE 1. EXPERIMENTAL AND ANALYTICAL ANALYSIS FLOWCHART 
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Abstract 

The best straight line through a set of experimental 
data is often obtained through the application of 
linear regression analysis, which provides the values 
of the slope and intercept for this line. The 
uncertainty intervals that should be associated with 
the values of the slope and intercept are also 
important and needed information. Standard 
statistical techniques to estimate the uncertainties in 
the slope and intercept are of limited use, primarily 
due to their assumptions which preclude their use 
with bias uncertainties. The approach to 
determining the uncertainty in the values of slope 
and intercept presented in this paper is based upon 
applying the uncertainty propagation equations to 
the regression analysis equations for the slope and 
intercept. This approach provides for the inclusion 
of precision uncertainties, bias uncertainties, and 
correlated bias uncertainties. Using a Monte Carlo 
type simulation technique, it is shown that this 
approach provides appropriate estimates of the 
uncertainty intervals in cases with bias uncertainties 
and precision uncertainties in both the X and Y 
measurements. 

Nomenclature 

Br=bias limit 
Bik=covariance estimator 
c=Y -intercept 
m=slope of line 

M=number of dependent variable measurements per 

independent variable 

N=number of data points 

Pi=precision limit 

^experimental result 
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Spstandard deviation 
X=independent variable 
X,=measured variables 
Y=dependent variable 
Up=uncertainty interval 
factual bias error 
9i=sensitivity, Eq.8 
Pbik=correlation coefficient 

Introduction 

In many experimental programs, the quantity of 
interest is not the data itself, it is the relationship 
within the data For example, a tensile test is often 
conducted to determine the modulus of elasticity for 
a specimen, where the modulus of elasticity is the 
slope of the stress-strain relationship that is assumed 
linear in the elastic region. This relationship can be 
expressed in the general form 

Y = mX + c (1) 

where m is the slope of the line and c is the intercept 
on the Y-axis. 

What are the uncertainties in the calculated 
values of m and c and why are they needed? A use 
which partially motivated this woric is an on going 
effort to determine the uncertainty from a Space 
Shuttle Main Engine prediction model which uses 
regression coefficients from experimental data. This 
requires knowledge of the uncertainty in the 
regression coefficients. 

The uncertainties in the slope, m, and the Y- 
intercept, c, are functions of the uncertainties in the 
deteiminations of the X and Y variables. The values 
of m and c are obtained by minimizing the sum of 
the squares of the deviations between the line and 
the data points, commonly known as the method of 
least squares. The development of the equations to 
calculate m and c can be found in statistics books 1,2 ' 5 
and only the equations will be presented here. For N 
(Xi, Yj) data pairs, the slope of the line, m, is 
determined from 
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and the intercept, c, is determined from 
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In this paper we present a methodology to 
determine the uncertainties in linear regression 
coefficients. The effectiveness of this methodology 
was analyzed using a Monte Carlo-type simulation 
assuming the true relationship between the X and Y 
variables is in fact linear. 


Experimental Uncertainty Analysis 

The background of the methodology to 
obtain uncertainty estimates and how they propagate 
through a given data reduction equation is not 
presented in this paper - only a brief overview is 
given. The reader is referred to references [4], [5], 
and [6] for a more thorough discussion. 

The word accuracy is generally used to 
indicate the relative closeness of agreement between 
an experimentally-determined value of a quantity 
and its true value. Error is the difference between 
the experimentally-determined value and the truth, 
thus as error decreases accuracy is said to increase. 
Only in rare instances is the true value of a quantity 
known. Thus, it is necessary to estimate error, and 
that estimate is called an uncertainty, U. 
Uncertainty estimates are made at some confidence 
level — a 95% confidence estimate, for example, 
means that the true value of the quantity is expected 
to be within the ±U interval about the 
experimentally-determined value 95 times out of 
100 . 

Total error can be considered to be 
composed of two components: a precision (random) 
component e and a bias (systematic) component (J. 
An error is classified as precision if it contributes to 
the scatter of the data; otherwise, it is a bias error. 
As an estimator of (3, a systematic uncertainty or bias 
limit, B, is defined. A 95% confidence estimate is 
interpreted as the experimenter being 95% confident 
that the true value of the bias error, if known, would 
fall within ±B. A useful approach to estimating the 


bias limit is to assume that the bias error for a given 
case is a single realization drawn from some 
statistical parent distribution of possible bias errors. 
For example, suppose a thermistor manufacturer 
specifies that 95% of samples of a given model are 
within ±1.0 C of a reference resistance-temperature 
(R-T) calibration curve supplied with the 
thermistors. One might assume that the bias errors 
(the differences between the actual, but unknown, 
R-T curves of the various thermistors and the 
reference curve) belong to a Gaussian parent 
distribution with a standard deviation b=0.5 C. 
Then the interval defined by ±B = ± 2b = ±1.0 C 
would include about 95% of the possible bias errors 
that could be realized from the parent distribution. 

As an estimator of the magnitude of the 
precision errors, a precision uncertainty or precision 
limit P is defined. A 95% confidence estimate of P 
is interpreted to mean that the ±P interval about a 
single reading of X* should cover the (biased) parent 
population mean 95 times out of 100. 

In nearly all experiments, the measured 
values of different variables are combined using a 
data reduction equation (DRE) to form some desired 
result. A general representation of a data reduction 
equation is 


r = r(Xi, X 2 , ■■■, Xj) (4) 

where r is the experimental result determined from J 
measured variables X t . Each of the measured 
variables contains bias errors and precision errors. 
These errors in the measured values then propagate 
through the data reduction equation, thereby 
generating the bias and precision errors in the 
experimental result, r. 

If the "large sample assumption” is made 6 , 
then the 95% confidence expression for U r becomes 

Ul = B 2 r + P 2 r (5) 

where we define the bias limit (systematic 

uncertainty) of the result as 
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and the precision limit (precision uncertainty) of the 
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and where B* is the 95% confidence estimate of the 
covariance appropriate for the bias errors in X s and 
X k and Pfc is the 95% confidence estimate of the 
covariance appropriate for the precision errors in X 4 
and X b and 
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The 95% confidence precision limit for a 
variable X s can be estimated as 

Pxi = 2 Sxi ( 9 ) 


where the sample standard deviation for X, is 
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and the mean value for X, is defined as 
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Typically, correlated precision uncertainties 
have been neglected so that the Pi’s in Eq. (7) are 
taken as zero, and that is assumed in the work 
reported here. For a thorough discussion of the 
estimation of the B&'s in Eq. (6), the reader is 
referred to Reference [7]. For the assumptions used 
in the work reported in this paper. 


Bik = Bi Bk (12) 


Regression Uncertainty 

The standard statistical techniques for the 
determination of uncertainties in regression 
coefficients have been mainly restricted to random 
uncertainties. Recently it has been accepted that bias 
uncertainties and particularly correlated bias 
uncertainties play an important role in engineering 
measurements and must be properly accounted for. 6 
As such, a methodology which properly accounts for 
the propagation of both precision and bias 
uncertainties through linear least squares regression 
equations has been the subject of recent work. 

Price 8 presents a suggested methodology which 
provides an uncertainty interval associated with a 
given predicted value of Y. Price’s method uses the 
traditional statistical techniques for the random 
uncertainties (standard error of the estimate 
equations) and a new technique to propagate bias 
uncertainties from elemental error sources in X and 
Y into the predicted T value. The ability of the 


methodology to handle correlated bias uncertainties 
is not explicitly discussed, nor is it shown that the 
methodology implicitly accounts for correlated bias 
uncertainty sources. 

Clark'' again utilizes the traditional statistical 
techniques to calculate the standard deviations of the 
slope and intercept due to the precision uncertainties. 
He then uses a separate technique to propagate the 
elemental bias uncertainties into a bias limit for the 
slope and intercept, but the technique does not 
account for correlated bias uncertainties. 

Montgomery 1 is a standard statistical linear 
regression analysis text and discusses the Ridge 
Regression method as an appropriate technique to 
use with biased data and is generally used with 
higher order regressions. The Ridge technique is still 
mainly interested in providing the “best fit” of the 
data and does not allow for the inclusion of bias 
limits estimated by non-statistical methods. 

The approach presented in this paper is an 
application of the uncertainty analysis methodology 
presented above to the regression equations for slope 
and intercept. 

Considering Eq.s (2) and (3) to be data 
reduction equations of the form 

m— Tri^X^yX^,..., Xf/,Y^, lj,...,iy) (13) 


and 


c —c(X^, X ^ , . . . , X^ , Y} , Y2 , . . , , Y n ) (14) 

and applying the uncertainty analysis equations, Eq.s 
(5)-(12), the most general form of the expression for 
the uncertainty in the slope of the line, m, is 
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where B y is the systematic uncertainty for the Y t 
variable, B x is the systematic uncertainty for the 
X i variable, B YYt is the covariance estimator for 
the correlated bias uncertainties in the Y- and K 
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variables, Bx y x t is the covariance estimator for 
correlated bias uncertainties in the X t and X k 
variables, B XY is the covariance estimator for the 
correlated bias uncertainties between X t and Y , . 
P y is the random uncertainty for the Y i variable, 
and P x is the random uncertainty for the X t 
variable. 

A similar expression for the uncertainty in the 
intercept is 
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The partial derivatives are 
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The equauons above show the most general form of 
the equations for the uncertainty in the slope and the 
intercept, allowing for correlation of bias errors 
among the different X s, among the different T s and 
also among the JTs and Ts.. If none of the 
systematic error sources are common between the X 
variables and the Y variables, the last term of the 
equations, the X-Y covariance estimator, is zero. 
This was the form of the equations used in the 
simulations. 


Monte Carlo Simulations 

Monte Carlo-type simulations are often used in 
uncertainty analysis to determine the effectiveness erf 
a particular uncertainty methodology. For this work, 
what is referred to as a Monte Carlo-type simulation 
simply means generating numbers to represent 
experimental data with some amount of error 
randomly obtained from a predefined error 
distribution population. Figure 1 is a schematic 
flowchart of the Monte Carlo simulation technique 
used. 

The Monte Carlo simulations were conducted in 
the following manner. “True" values for data from a 
linear relationship with specified coefficients were 
determined. The word true is emphasized to 
indicate that it represents the actual physical 
quantity of the parameter if it could be measured 
without any bias error or precision error, which is 
always an unobtainable value. The two-sigma (2 
standard deviation or 95% confidence) bias limits 
and precision limits for each X and Y were then 
specified. The errors in each variable were assumed 
to come from these normally distributed error 
populations with the specified standard deviations. 
A random value for each bias error and precision 
error was found from a Gaussian random deviate 
generator subroutine using the specified standard 
deviations. The Gaussian deviates have a mean erf 
zero and an equal probability of being positive or 
negative. In the cases presented in this paper all erf 
the X bias errors were defined as being correlated 
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and of the percent of full scale type, so the same 
random deviate was used for each X,. The same held 
true for the 7, variables, but with a different bias 
error than that in the Xr s. Precision errors were 
obtained by sampling the precision error populations 
repeatedly to obtain independent random deviates for 
each Xi and 7*. For each X„ 7, pair, the individual 
bias errors and precision errors were then summed 
and added to the true value to obtain a data point 
with errors from the specified error populations. 
These data points were then used in the linear least 
squares equations to obtain the value of the 
regression coefficients. These values represent the 
regression coefficients of the experiment when the 
bias and precision errors are present. 

A 95% confidence uncertainty interval for the 
result was calculated from the uncertainty 

propagation analysis equations for m and c. A ±U m 
interval was placed around the slope coefficient 
value, m , and if the true value of the slope was found 
to be within the interval a counter was incremented. 
A similar procedure was used with c. This 
procedure was repeated 10,000 times and the percent 
coverage, or number of times the true result was 
within the estimated interval, was determined. 
Using this procedure, the effectiveness of the 
uncertainty propagation equation could be 

investigated by checking whether or not the true 
value is within the 95% confidence uncertainty 
interval about the measured result 95% of the time. 

A useful statistic from the simulation is the 
uncertainty ratio, the ratio of the average uncertainty 
intervals for the regression coefficients from the 
10,000 iterations divided by the true 95% confidence 
intervals. The true 95% uncertainties are calculated 
as twice the sample standard deviations, S m and 5 C> , 
from the 10,000 samples of the regression 
coefficients. The sample standard deviations from 
the 10,000 sample population can be expected to be 
good representations of * the actual standard 
deviations of the infinite population of the 
coefficients with the elemental uncertainty sources as 
defined. An uncertainty ratio of or near unity shows 
that the uncertainty methodology works for the 
particular case, with values greater than one 
meaning an overprediction and values less than one 
meaning an underprediction. 

RESULTS 

Simulations were performed with errors only in 
7, multiple readings of 7 at each X setpoint, and with 
errors in both X and 7 for single and multiple 
readings per setpoint. The coefficients m and c were 


then determined with the standard linear regression 
equations, equations (2) and (3), and the uncertainty 
interval was calculated. 

The simulations progressed from simple to 
complex, with the increasingly more complex 
simulations being a closer representation of an actual 
experimental program. Table 1 shows the input data 
used as the truth , Table 2 shows a summary of the 
simulations for the number of data points and the 
method used for the precision limits, and Table 3 
shows the 2a (95%) intervals used to determine the 
bias and precision errors from the Gaussian random 
deviate routines. The percent coverage and the 
uncertainty ratios were determined both with the 
propagation equations including the effect of 
correlated bias uncertainties and neglecting them. 

General 

A few comments and conclusions can be drawn 
from the results of all the simulations. 

The first and third simulations are essentially 
verifications of the methodology. The actual 2a 
(95%) confidence intervals for both the precision 
uncertainties and the bias uncertainties were used in 
the uncertainty propagation equations, including 
correlated bias effects. Since they yielded 
approximately 95% coverage the proposed 
methodology is demonstrated as appropriate. 

In the simulations presented in this paper all of 
the bias uncertainties are a fixed value, or “percent 
of full scale.” With the bias uncertainties being the 
same value for each data pair, the linear least 
squares line through the data simply translates 
vertically or horizontally from the true line. The 
value of the slope remains the same and there is no 
systematic uncertainty in the slope. The value of the 
intercept will change and will accordingly have an 
uncertainty. In the numerical calculation of the 
uncertainty of the slope when all of the bias 
uncertainties are correlated and there are no 
precision uncertainties, an error due to truncation 
and round-off was encountered. This problem was 
avoided by including a very small precision 
uncertainty during the cases of dominant bias 
uncertainties. 

When the bias uncertainties are of a “percent of 
reading” nature, so that the bias limit is a function of 
the magnitude of the variable, the slope MU have an 
associated non-zero systematic uncertainty. This 
case is not considered in this paper. 

Simulation 1 
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The first simulation included a systematic error 
and a random error for Y only. The systematic errors 
were totally correlated, that is the same value of 
systematic error was used for each Yj. A random 
error, from the same error population distribution, 
but with a different value was used at each data pair. 
The uncertainty intervals were determined using 
equations (15) and (16) with all X uncertainties 
defined as zero. The simulation was repeated for the 
cases of dominant systematic error, dominant 
random error, and for systematic error proportional 
to random error. 

Figure 2 shows the effectiveness of the 
regression coefficient uncertainty propagation 
equation for the slope as a function of type of 
dominant uncertainty (bias or precision) and as a 
function of the number of data pairs. This figure 
clearly shows that the propagation equation provides 
the appropriate 95% coverage for all cases. Similar 
results yielding 94% to 96% coverage for all types of 
dominant uncertainties was obtained for the intercept 
uncertainty propagation equation. Thus, when the 
correct values for the precision limit and bias limit 
are used in the propagation equations the correct 
uncertainty interval is obtained. Obtaining the 
correct value for the precision limit is only possible 
if an infinite amount of previous information is 
available upon which to determine the precision 
limits. 

Simulation 2 

The second simulation still only included errors 
in the Y variable, however M additional readings of 
the Y variable were generated at each X variable 
setpoint. The same systematic error was used for all 
of the readings within each iteration, but a different 
precision error was used for each Y. The m anc c 
coefficients for the N times M data pairs was found 
with the linear regression equations. The 
uncertainty interval was then calculated with the 
precision limit, P y , determined by calculating the 

standard deviation of the M Y readings at each X 
setpoint and using the large-sample approximation 
(t=2). 

Figures 3 through 6 show the effectiveness of 
the regression coefficient propagation equations as a 
function of the number of points and for the 
dominant types of uncertainties. Figure 3 shows that 
approximately 94% to 96% coverage is obtained 
after the total number of data points reaches about 25 
or 30. A similar result was obtained for the intercept 
uncertainty for dominant precision uncertainties. 
Figures 4 and 5 show the percent coverage for the 


case of comparable magnitude bias and precision 
errors for the slope uncertainty and intercept 
uncertainty, respectively. These figures show that 
for even very small numbers of data points that the 
percent coverage is about 93%. In reference [6] it is 
discussed that the difference between 93% and 95% 
is essentially irrelevant since estimates for the bias 
limits cannot be made to that accuracy. When the 
systematic uncertainty is dominant (and the bias 
limit is estimated correctly) about 95% coverage is 
obtained, as shown in figure 6 for the intercept. 

Simulation 3 

The third simulation is similar to the first, but 
errors in both the X and Y variables are included. A 
single experimental (X,Y) data pair is determined at 
each X variable setpoint. The uncertainty intervals 
were determined from Equations (15) and (16). All 
systematic errors in the X variables are correlated 
and all systematic errors in the Y variables are 
correlated, but there is no correlation of systematic 
errors between the X and the Y variables. The true 
95% confidence intervals for both the bias and 
precision uncertainties are used in the uncertainty 
propagation equations. 

With the same three variations of dominant 
uncertainty types, bias, comparable, and precision, 
the uncertainty in the slope and the uncertainty in 
the Y-intercept the coverage is essentially the desired 
95% and the plot appears identical to Figure 2. 

Simulation 4 

The fourth simulation is the most general. It 
includes systematic and random errors in both X and 
Y and there are M sets of (X,Y) experimental data at 
each nominal X setpoint. The uncertainty intervals 
were determined from Equations (15) and (16). As 
in the third simulation, the systematic errors in X are 
correlated and the systematic errors in Y are 
correlated, but there is no correlation of systematic 
errors between X and Y. The precision uncertainty is 
determined by calculating the standard deviation of 
the M readings of both variables at each X setpoint 
and using the large sample approximation at each 
setpoint. 

Figure 7 shows the percent coverage for 
dominant precision uncertainties, for greater than 
about 15 total data points (3 Y data points at 5 X 
setpoints or 5 Y data points at 3 X setpoints) the 
“large-sample approximation” in the regression 
propagation equations provides coverage in excess of 
about 93%. A similar plot is obtained for the 
intercept. Plots essentially identical to Figures 4 and 
5 for comparable bias and precision uncertainties 
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and Figure 6 for dominant bias uncertainty were 
obtained for this simulation, but not shown here. 


Simulation 5 

The final simulation has the same form of 
experimental data generation as the fourth 
simulation, but the linear regression is applied to the 
mean of the M readings at each setpoint 
Accordingly, the precision uncertainty is determined 

using the standard deviation of the mean, 5, / Va 7, 
and the large-sample approximation so that the 
precision limits are 



( 21 ) 


and 



( 22 ) 


where the index i represents the i ,h X setpoint. 

This methodology produces a higher coverage 
for a smaller number of data pairs, as shown in 
Figures 8 through 10. In the case of dominant 
precision uncertainties, Figures 8 and 9, the desired 
coverage of about 95% is reached with only a few 
measurements at each X setpoint. When bias 
uncertainties are included and of comparable 
magnitude to the precision uncertainties, and the 
correct bias limits have been used, the percent 
coverage is always greater than 92% and essentially 
in the desired 94% to 96% range, as shown in Figure 
10. When the bias uncertainties are the dominant 
uncertainty source the coverage is again essentially 
95% for both m and c and a plot similar to Figure 6 
is obtained. 


Effect of Correlation 

The effect of correlated bias uncertainties has 
been a neglected part of uncertainty analysis until 
relatively lately, and a new method of accounting for 
correlated bias uncertainties was recently presented. 
In many, if not most situations where a regression is 
performed the test data will come from the same test 
apparatus. Since the data is taken with the same 
equipment the bias errors in the experimental data 
will be from the same sources and will therefore be 
correlated. As part of the Monte Carlo simulation 
the percent coverage and the uncertainty ratio were 
determined with the regression propagation 
equations both including and excluding the 
correlated bias uncertainty terms. Figures 11 and 12 


show the uncertainty interval ratio (the calculated 
uncertainty interval divided by the actual 95% 
uncertainty interval) for the slope and intercept, 
respectively. The data from all of the simulations 
and all of the experiments within a given simulation 
type were grouped in categories of types of dominant 
uncertainties: bias dominant, precision dominant, or 
bias and precision of comparable magnitudes. 
Within each of the three dominance categories the 
results were randomized to confound the trends 
within the data. Figure 11 shows the uncertainty 
ratio for the slope, m. In the experiments where the 
bias uncertainties were dominant ignoring the effect 
of correlation provided a dramatic overestimate of 
the uncertainty interval, by several orders erf 
magnitude, since the actual uncertainty is of very 
small magnitude in these cases. In the experiments 
where the bias uncertainties and the precision 
uncertainties are of the same order of magnitude the 
uncertainty interval calculated ignoring the effect of 
correlated biases is still dramatically overestimated. 
A similar result is seen for the uncertainty ratio for 
the Y-intercept, although the overestimate is not 
quite as dramatic. Both figures also show that when 
precision uncertainties are dominant the appropriate 
uncertainty interval uncertainty interval is obtained 
even if the effect of correlation is ignored in the 
propagation equation. 

These results show quite vividly the impact of 
not accounting for correlated bias uncertainties when 
they actually exist within the data. This also 
demonstrates the potential error involved with using 
the traditional, statistically based regression 
uncertainty methods which ignore bias uncertainties. 

Conclusion 

The techniques developed in this paper provide 
an engineering method to determine the uncertainty 
in linear regression coefficients for data containing 
both random uncertainties, systematic uncertainties, 
and correlated systematic uncertainties. There are 
many applications of linear regression analysis 
uncertainty which were not studied as a part of this 
effort. Work is continuing in several areas. 
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(60,60) 
(70,70) 
(80, 80) 
(90, 90) 
(100, 100) 
(110, 110) 
(120, 120) 
(130, 130) 
(140, 140) 


Table 1. True (X,Y) data for all simulations (true 
line:Y=X+0) 
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1 

2c 
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2S y 
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1 
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3,5, 7,9 

3,5,7, 9, 
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2S y , 2 Sx 
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2 Sy 2 S x 

>fii’ 4 m 
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Sim. 

Bias 

Dominant 
By, Py, 
Bx» Px 

Biass 
Precision 
By, Py, 
BxlPx 

Precision 
Dominant 
By, Py, 
Bx* Px 

1 

5, 0.05, 

2.5, 2.5, 

0,5, 


n/a*, n/a* 

n/a*, n/a* 

n/a*, n/a* 

2 

5, 0.05, 

2.5, 2.5, 

0,5, 


n/a*, n/a* 

n/a*, n/a* 

n/a*, n/a* 

3 

5, 0.05, 

2.5, 2.5, 

0, 5, 


5, 0.05 

2.5, 2.5 

0,5 

4 

5, 0.05, 

2.5, 2.5, 

0, 5, 


5, 0.05 

2.5, 2.5 

0,5 

5 

5, 0.05, 

2.5, 2.5, 

0, 5, 


5, 0.05 

2.5, 2.5 

0,5 


Table 3. Simulation Uncertainty Input Data. 

n/a* • These simulations did not include any uncertainty in 


X. 
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Figure 1. Monte Carlo Simulation Flowchart 









Figure 2. Simulation 1. Plot of Coverage of Slope vs 
Type of Dominant Uncertainty and Number of Data 
Uncertainties only in Y. 



Figure 3. Simulation 2. Plot of Coverage of Slope for 
Dominant Precision Uncertainties vs Number of Data. 
Uncertainties only in Y. 



Figured Simulation 2. Plot of Coverage of Slope for 
Comparable Bias and Precision Uncertainties vs 
Number of Data. Uncertainties only in Y. 


Figure 5. Simulation 2. Plot of Coverage of Intercept 
for Comparable Bias and Precision Uncertainties vs 
Number of Data. Uncertainties only in Y. 
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Figure 6. Simulation 2. Plot of Coverage of Intercept 
for Dominant Bias Uncertainties vs Number of Data. 
Uncertainties Only in Y. 


Figure 8. Simulations. Plot of Coverage of Intercept 
for Dominant Precision Uncertainties vs Number of 
Data. Uncertainties in X and Y 


Figure 7. Simulation 4. Plot of Coverage of Slope for 
Dominant Precision Uncertainties vs Number of Data. 
Uncertainties in X and Y. Plot identical for Intercept 



»<rfy«p«x 

Figure 9. Simulation 5. Plot of Coverage of Slope for 
Dominant Precision Uncertainties vs Number of Data. 
Uncertainties in X and Y 
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FigurelO Simulations. Plot of Coverage of Slope for 
Comparable Bias and Precision Uncertainties vs 
Number of Data. Uncertainties in X and Y. 
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Figure 11. Plot of uncertainty ratios vs dominant uncertainty for slope uncertainty. 
With and without considering correlated bias uncertainties. 
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Figure 12. Plot of uncertainty ratios vs dominant uncertainty for Y-intercept uncertainty. 
With and without considering correlated bias uncertainties. 
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Appendix V 

Elemental Bias Limits Used in Venturi Uncertainty Determination 


Venturi # 139 (PE 


Source 

Description 

Bp=l% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bap=1% 

Engineering estimate 

Bcd=002 

~2%, Rocketdyne report, Cd Extrapolation, Cryogenic application 

Bdi= 0.0005 in 

Machining tolerance, Rocketdyne report 

Bd 2=0.0005 in 

Machining tolerance, Rocketdyne report 

B a =6.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.l LPOTP Inlet Venturi Bias Limit Inputs 


Venturi # 20 (PID 8801) LPFT Inlet 

Source 

Description 

B P =1% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bai^I^ 

Engineering estimate 

Bcd=0.02 

~2%, 0.8% from RD report, but only cal’d at 7.6% Re 

Bdi= 0.0005 in 

Machining tolerance, Rocketdyne report 

Bd 2=0 0005 in 

Machining tolerance, Rocketdyne report 

B a =4.2e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.2 LPFT Inlet Venturi Bias Limit Inputs 


Venturi #397 (PIE 


Source 

Description 

Bp=l% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bap=1% 

Engineering estimate 


~3%, 2% from RD report, but only cal’d at 1 .7% Re 

Bdi= 00005 in 

Machining tolerance, Rocketdyne report 

Bd 2=0.0005 in 

Machining tolerance, Rocketdyne report 

B a =5.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.3 CCV Inlet Venturi Bias Limit Inputs 
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Venturi # 426 (PID 8819) LPOP Discharge 

Source 

Description 

Bp=l% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bap=1% 

Engineering estimate 

Bcm = 0.020 

~2%, 1% from RD report, but only cal’d at 4.8% Re 

Bdi= 0.0005 in 

Machining tolerance, Rocketdyne report 

Bd2=0 0005 in 

Machining tolerance, Rocketdyne report 

B a =5.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.4 LPOP Discharge Venturi Bias Limit Inputs 


Venturi # 268 (PID 8804) OPB LOX 

Source 

Description 

Bp=l% 

Sverdrup Report i 

Bt=1% 

Sverdrup Report 

B^p=1% 

Engineering estimate 

Bcd=0.021 

~2%, .8% from RD report, but only cal’d at 18% Re 

Bdi == 0.0005 in 

Machining tolerance, Rocketdyne report 

Bd 2 = 0.0005 in 

Machining tolerance, Rocketdyne report 

B o =6.7e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.5 OPB LOX Venturi Bias Limit Inputs 


Venturi #271 (PID 8805) OPB LH 2 

Source 

Description 

Bp=l% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bap = 1% 

Engineering estimate 

Bca^OOS 1 

~3%, .2.1% from RD report, but only cal’d at 1.86% Re 

Bdi^O. 0005 in 

Machining tolerance, Rocketdyne report 

Bd 2 == 0*0005 in 

Machining tolerance, Rocketdyne report 

B a =3.2e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.6 OPB LH 2 Venturi Bias Limit Inputs 
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Venturi #296 (PIP 8810) 


FPBLOX 


Source 


Description 


Bp=l% 


Sverdrup Report 


Bt=1% 


Sverdrup Report 


Bap=1% 


Engineering estimate 


Bcd=0.021 


~2%, .1.25% from RD report, but only cal’d at 20% Re and 
poor curve-fit 


Bdi=O.OOQ 5 in 


Machining tolerance, Rocketdyne report 


Bp 2 ~~ 0.0005 in 


Machining tolerance, Rocketdyne report 


B a — 6.7e-07 m/in-R 


~10%, Engineering estimate 


Table A-V.7 FPB LOX Venturi Bias Limit Inputs 



L5) Nozzle Coolant #1 

Source 

Description 

Bp=l% 

Sverdrup Report 

Bt=1% 

Sverdrup Report 

Bap=1% 

Engineering estimate 

Bcd = 0020 

~2%, .1.15% from RD report, but only cal’d at 3.40% Re 

B D i=0.0005 

Machining tolerance, Rocketdyne report 

Bd2 == 00005 

Machining tolerance, Rocketdyne report 

B o =5.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.8 Nozzle Coolant #1 Venturi Bias Limit Inputs 


Venturi #3542 (PID8816) Nozzle Coolant #2 

Source 

Description 

Bp=l% 

Sverdrup Report 

Bx=l% 

Sverdrup Report 

Bap=1% 

Engineering estimate 

Bcd=0.020 

~2%, .1.15% from RD report, but only cal’d at 3.40% Re 

Bdi= 0.0005 in 

Machining tolerance, Rocketdyne report 

Bd 2=0.0005 in 

Machining tolerance, Rocketdyne report 

B a =5.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.9 Nozzle Coolant #2 Venturi Bias Limit Inputs 
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Venturi # 3543 (PID 8817) Nozzle Coolant #3 

Source 

Description 

Bp = 1 % 

Sverdrup Report 

Bf=l% 

Sverdrup Report 

Bap == 1% 

Engineering estimate 

Bcd=0.020 

~2%, .1.15% from RD report, but only cal’d at 3.40% Re 

Bdi= 0.0005 in 

Machining tolerance, Rocketdyne report 

Bdz^ OOOS in 

Machining tolerance, Rocketdyne report 

B a =5.5e-07 in/in-R 

~10%, Engineering estimate 


Table A-V.10 Nozzle Coolant #3 Venturi Bias Limit Inputs 
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Appendix VI 

Engine Configuration and Averaged Test Data 
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ENGINE 

TEST_NO 

TESTDATE 

ADUR 

PWRHEAD 

MAININJ 

MCC 

NOZZLE 

LPFTPUNIT 

HPFTPUNIT 

LPOTPUNIT 

HPOTPUNIT 

LPFTPSERIAL 

HPFTPSERIAL 

LPOTPSERIAL 

HPOTPSERIAL 

3001 

a50039 

1/14/93 

210 

9 

2008 

6001 

2026 

4004R2 

4406R2 

2218 

2722 

4873079 

4108675 

4876205 

4871900 

3001 

a50040 

2/11/93 

180 

9 

2008 

6001 

2026 

4004R2 

4406R2 

2218 

2722 

4873079 

4108675 

4876205 

4874500 

3001 

a50041 

4/8/93 

203 

9 

2008 

6001 

2026 

4004R2 

4406R2 

2218 

8105R5 

4873079 

4108675 

4876205 

BUF500 

3001 

a50042 

4/30/93 

203 

9 

2008 

6001 

2026 

4004R2 

8006R1 

2218 

8105R5 

4873079 

XFPW06 

4876205 

BUF500 

3001 

850043 

7/8/93 

190 

9 

2008 

6001 

2026 

4004R2 

4406R2 

2218 

8202R4 

4873079 

4108675 

4876205 

XOPW02 

3001 

a50044 

9/16/93 

197 

9 

2008 

6001 

2026 

4004R2 

4406R2 

2217 

8202R5 

4873079 

4108675 

4876390 

XOPW02 

3001 

a50045 

12/22/93 

200 

9 

2008 

6001 

2026 

4004R2 

2227R2 

2217 

231 5R3 

4873079 

4866524 

4876390 

4877218 

3001 

a50046 

1/27/94 

200 

9 

2008 

6001 

2026 

4004R2 

2227R2 

4204 

231 5R3 

4873079 

4866524 

4874193 

4877218 

3001 

a50047 

2/15/94 

185 

9 

2008 

6001 

2026 

4004R2 

2227R2 

4204 

231 5R3 

4873079 

4866524 

4874193 

4877218 

3001 

a50048 

3/18/94 

185 

9 

2008 

6001 

2026 

4004R2 

2227R2 

4401 R1 

231 5R3 

4873079 

4866524 

4876580 

4877218 

3001 

850049 

4/12/94 

180 

9 

2008 

6001 

2026 

4004R2 

2227R2 

4401 R1 

231 5R3 

4873079 

4866524 

4876580 

4877218 

3001 

a50050 

7/6/94 

79 

9 

2008 

6001 

2026 

4004R2 

8106 

4204 

8107R10 

4873079 

XFPW06 

4874193 

XOPW07 

3001 

a50051 

7/20/94 

40 

9 

2008 

6001 

2026 

4004R2 

8003 

4204 

8107R10 

4873079 

XFPW03 

4874193 

XOPW07 


Table A-VL1 Engine Hardware Component Configuration 
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TTB045 

TTB047 

TTB048 

TTB049 

TTB050 

4022 

4022 

4027 

4023 

3999 

254 

262 

276 

251 

256 

463 

463 

465 

466 

474 

5309 

5308 

5318 

5307 

5224 

396 

398 

395 

387 

386 

90 

91 

90 

91 

93 

5460 

5439 

5468 

5461 

5384 

112 

128 

105 

110 

120 

89 

90 

89 

89 

92 

5440 

5430 

5447 

5440 

5385 

108 

114 

109 

108 

118 

91 

91 

91 

91 

93 

5430 

5429 

5445 

5440 

5392 

99 

131 

141 

135 

121 

88 

89 

88 

j 88 

91 

3687 

3659 

3646 

3653 

3796 

192 

203 

190 

189 

207 

192 

192 

191 

193 

188 

3739 

3718 

3692 

3696 

3718 

249 

248 

254 

246 

”245 

194 

194 

193 

194 

191 

6856 

6850 

6802 

6812 

6621 

114 

121 

113 

114 

134 

216 

216 

215 

216 

215 " 

4948 

4945 

4958 

4953 

4888 

231 

231 

247 

233 

261 

277 

278 

278 

278 

281 

6908 

6897 

6837 

6843 

6653 

117 

109 

101 

121 

97 

206 

206 

205 

207 

205 



AVG 

S 

ggg 


4039.2 

40.4 

1.0 


263.1 

11.5 



468.2 

6.1 

1.3 


5336.9 

85.4 

1.6 


377.2 

27.5 

ME 


92.4 

1.9 

2.1 


5479.4 

74.7 

1.4 


121.8 

14.6 

12.0 


97.1 

10.4 

10.7 


5466.2 

74.5 

1.4 


117.9 

12.1 

10.3 


93.0 

2.3 



5464.6 

73.6 

1.3 


145.1 

25.8 

17.8 


156.2 

204.1 

130.7 


3708.8 

67.8 

1.8 


190.9 

8.9 

ME 


190.7 

1.9 

1.0 


3705.3 

23.6 

0.6 


245.9 

6.3 

2.5 


192.3 

1.6 

0.9 


6704.2 

138.9 

2.1 


105.8 

72.1 

68.1 

mwsmwKEwm 

IH 

52.0 

1.0 


245.8 

18.4 

mm 


280.9 

3.9 

1.4 


6760.0 

134.8 

2.0 


115.0 


9.0 


205.5 

ME 

0.6 






































































































































































































































































































































































































































Appendix VII 

Calculated Venturi Mass Flowrate Data 
and Test-to-Test Standard Deviations 
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Table A-VH1 Venturi Mass Flowrates for Tests TTB 039-050 
and Standard Deviations @ 100% RPL 



Calcuated Venturi Mass Flow Rates @ 100% RPL data (same engine 
components; Rocketdyne HPOTP) 



ttb049 

ttb048 

ttb047 

ttb045 

8801 

26.6 

27.8 

27.2 

26.8 

8818 

71.9 

72.6 

72.8 

72.7 

8815 

12 . 6 " 

12.4 

13.6 

12.8 

8816 

12.3 

12.4 

12.6 

12.3 

8817 

13.8 

14.0 

13.6 

11.8 

8802 

174 

175 

181 

176 

8819 

907 

924 

913 

914 

8804 

22.2 

22.1 

22.8 

22.2 

8805 

32.9 

33.8 

32.8 

32.9 

8810 

64.0 

58.7 

60.9 

63.1 

total 

1337 

1352 

1350 

1344 




Table A-VEL2 Venturi Mass Flowrates for Tests TTB 045, 0048-049 
and Standard Deviations @ 100% RPL 
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Calcuated Venturi Mass Flow Rates (Ib/sec) @ 104% RPL 
data 


PID 



ttb048 

ttb047 


ttb044 

ttb043 

ttb042 

8801 

28.0 

28.1 

29.2 

28.5 

28.3 

28.4 

29.3 

27.4 

8818 

73.8 

74.0 

74.8 

75.0 

75.0 

74.1 

74.2 

72.1 

8815 

13.8 

13.9 

13.3 

13.8 

13.1 

13.4 

12.7 

13.4 

8816 

13.5 

13.5 

13.1 

13.3 

12.7 

13.0 

13.3 

13.3 

8817 

13.6 

14.5 

14.6 

14.0 

6.0 

14.3 

WKE 

15.0 

8802 

188 

178 

179 

185 

180 

181 

182 

177 

8819 

948 

943 

962 

951 

949 

926 

932 

933 

8804 

26.1 

23.9 

24.0 

24.5 

23.9 

27.2 

27.4 

27.4 

8805 

35.7 

33.7 

34.8 

33.5 

34.0 

36.1 

36.2 

32.3 

8810 

61.3 

67.6 

63.2 

65.4 

67.6 

69.0 

66.6 

65.6 

total 

1402 

1390 

1408 

1404 

1389 

1382 

1381 

1377 


Table A-VIL3 Venturi Mass Flowrates for Tests TTB 040-050 


and Standard Deviations @ 104% RPL 


Calcuated Venturi Mass Flow Rates (Ib/sec) @ 104% RPL data (same 
engine components; Rocfcetdyne HPOTP) 



ttb049 

ttb048 

ttb047 

ttb045 


avg 

S 

3 (%) 

8801 

28.1 

29.2 

28.5 

28.3 


28.5 

0.50 

1.7 

8818 

74.0 

74.8 

75.0 

75.0 


74.7 

0.47 

0.6 

8815 

13.9 

13.3 

13.8 

13.1 


13.5 

0.38 

2.8 

8816 

13.5 

13.1 

13.3 

12.7 


13.1 

0.32 

■EB 

8817 

14.5 

14.6 

14.0 

6.0 


12.3 

4.18 

34.0 

8802 

178 

179 

185 

180 


180.4 

3.07 

1.7 

8819 

943 

962 

951 

949 


951.1 

7.77 

0.8 

8804 

23.9 

24.0 

24.5 

23.9 


24.1 

0.27 

1.1 

8805 

33.7 

34.8 

33.5 

34.0 


34.0 

0.57 

1.7 

8810 

67.6 

63.2 

65.4 

67.6 


65.93 

2.11 


total 

1390 

1408 

1404 

1389 

1398 

9.58 

0.7 


Table A-VIL2 Venturi Mass Flowrates for Tests TTB 045, 048-049 
and Standard Deviations @ 100% RPL 
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Calcuated Venturi Mass Flow Rates (Ib/sec) @ 109% RPL data 


PID 



ttb048 

ttb047 


ttb044 

ttb043 

ttb042 




avg 

S 


8801 

29.8 

29.9 

31.0 

30.3 

30.0 

30.1 

31.1 

28.6 

30.4 



30.1 

IKS) 

mm 

8818 

76.6 

76.6 

77.4 

77.7 

77.9 

76.9 

76.9 

75.8 

76.5 



76.9 

0.65 

0.9 

8815 

14.7 

14.8 

14.4 

14.5 

13.8 

14.6 

13.9 

13.0 

14.0 



14.2 

0.57 

4.0 

8816 

14.4 

14.1 

14.2 

14.0 

13.4 

14.2 

14.3 

13.0 

13.7 



13.9 

0.45 

3.3 

8817 

14.3 

15.4 

15.5 

14.9 

8.3 

16.2 

8.4 

15.4 

16.0 



13.8 

3.15 

22.8 

8802 

195 

182 

184 

190 

185 

187 

186 

184 

181 



186 

4.47 

mm 

8819 

993 

989 

1005 

992 

991 

970 

973 

972 




987 

12.1 

1.2 

8804 

28.6 

26.6 

26.3 

27.0 

26.2 

29.2 

29.2 

29.8 

25.3 



27.6 

1.65 

6.0 

8805 

37.0 

34.7 

36.2 

34.5 

35.3 

38.0 

38.5 

33.2 

36.7 



36.0 

1.73 

4.8 

8810 

66.0 

72.6 

68.5 

71.0 

72.9 

74.6 

72.6 

70.0 

76.3 



71.6 

3.13 

mm 

total 

1470 

1455 

1472 

1466 

1453 

1450 

1443 

1435 

1466 



1457 

12.68 

0.9 

Table A-V115 Venturi 

MassF 

lowrates for Tests TTB 04 

0-050 






and Standard Deviations @ 109% RPL 








Calcuated Venturi Mass Flow Rates (Ib/sec) @ 109% RPL data 





(same engine components; Rocket 

dyne H 

POTP) 








ttb049 

ttb048 

ttb047 

ttb045 


avg 

S 

S(%) 





8801 

29.9 

j 31.0 

30.3 

30.0 


30.3 

0.47 

1.5 





8818 

76.6 

77.4 

77.7 

77.9 


77.4 

0.55 

0.7 





8815 

14.8 

14.4 

14.5 

13.8 


14.4 

6.41 

2.9 





8816 

14.1 

14.2 

14.0 

13.4 


13.9 

0.35 

2.5 





8817 

15.4 

15.5 

14.9 

8.3 


13.5 

3.47 

25.6 





8802 

182 

184 

190 

185 


185 

3.56 

1.9 





8819 

989 

1005 

992 

991 


994 

7.37 

0.7 





8804 

26.6 

26.3 

27.0 

26.2 


26.5 

0.36 

1.3 





8805 

34.7 

36.2 

34.5 

35.3 


35.2 

0.74 

2.1 





8810 

72.6 

68.5 

71.0 

72.9 


71.2 

2.00 

2.8 





total 

1455 

1472 

1466 

1453 


1462 

8.95 

0.6 






Table A-VEL6 Venturi Mass Flowrates for Tests TTB 045, 048-049 
and Standard Deviations @ 109% RPL 
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Appendix VULl 

Venturi Mass Flowrate Data and Calculated Bias Limits 
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ttb050/100% 

W 

Ib/sec 

Bw 

Ib/sec 

26.61 

0.55 

71.41 

2.17 

13.08 

0.27 

12.78 

0.27 

12.98 

0.27 

183.29 

3.8 

910.39 

18.69 

24.05 

0.5 

34.39 

0.7 

57.47 

1.19 


ttt>050/104% 


Bw 

Ib/sec 


■ 

W 

Ib/sec 

8801 

28.02 

8818 

73.78 

8815 

13.77 

8816 

13.48 

8817 

13.58 

8802 

188.29 

8819 

947.63 

8804 

26.09 

8805 

35.68 

8810 

61.29 


0.291 2.1 


3.9 2.07 


19.43 2.05 


0.54 2.07 


0.72 2.02 


1.27 2.07 


Table A-VIEL1 Venturi Flowrates and Bias Limits for test TTB050 


ttb050/109% 

W 

Ib/sec 

Bw 

Ib/sec 

29.78 

0.61 

76.57 

2.33 

14.67 

0.3 

14.37 

0.3 

14.33 

0.3 

195.4 

4.05 

993.1 

20.35 

28.55 

0.59 

37.04 

0.75 

66.03 

1.37 


Bw (%) 


ttb049/100% 


W 

Bw 

Ib/sec 

Ib/sec 

26.63 

0.55 

71.89 

2.19 

12.61 

Eli] 

12.31 

KE3 

13.75 

0.29 

173.9 

3.6 

907.18 

18.65 

22.16 

0.46 

32.92 

0.67 

64 

1.33 

Tabic 

ttb048/100% 

W 

Bw 

Ib/sec 

Ib/sec 

27.79 

0.57 

72.63 

2.21 

12.35 

0.26 

12.37 

0.26 

14.03 


174.63 

3.62 

923.88 

18.98 

22.05 

0.46 

33.84 

0.68 

58.65 

1.22 



ttb049/104% 

W 

Bw 

Ib/sec 

Ib/sec 

28.12 

0.58 

73.99 

2.25 

13.91 

0.29 

13.47 

0.28 

14.45 

0.3 

177.65 

3.68 

943.16 

19.36 




23.92 


33.74 


67.59 


Venturi Flowrates and Bias Limits for test TTB049 


ttb049/109% 

W 

Ib/sec 

Bw 

Ib/sec 

29.91 

0.61 

76.63 

2.33 

14.78 

0.31 

14.12 

0.3 

15.37 

0.32 

182.2 

3.77 

988.5 

20.28 

26.56 

0.55 

34.73 

0.7 

72.56 

1.5 



ttb048/104% 


W 

Ib/sec 


Bw Bw (%) 
Ib/sec 


8801 

29.23 

■EE 

8818 

74.76 

2.27 

8815 

13.31 

0.28 

8816 

13.11 

0.28 

8817 

14.64 

0.31 

8802 

179.17 

3.71 

8819 

961.66 

19.73 

8804 

23.97 

KEI 

8805 

34.81 

0.7 

8810 

63.16 

1.31 



ttb048/109% 


W 

Ib/sec 

Bw 

Ib/sec 

Bw (%) 

30.95 

0.63 

2.04 

77.37 

2.35 

3.04 

14.38 

0.3 

2.07 

14.15 

0.3 

2.11 

15.52 

0.33 

2.1 

184 

3.81 

2.07 

1005 

20.6 

2.05 

26.26 

0.54 

2.07 

36.19 

0.73 

2.02 

68.5 

1.42 

2.07 


Table A-VHL3 Venturi Flowrates and Bias Limits for test TTB048 
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Table A-Vm.4 Venturi Flowrates and Bias Limits for test TTB047 


ttb045/100 I ttb045/104 

W I Bw Bw (%) W I Bw Bw (%) 

Ib/sec Ib/sec Ib/sec Ib/sec 


8801 

26.84 

0.55 

2.05 


8801 

28.25 

0.58 

2.05 

8818 

72.72 

2.21 

3.04 


8818 

75.02 

2.28 

3.04 

8815 

12.75 

0.26 

2.07 


8815 

13.08 

0.27 

2.07 

8816 

12.33 

0.26 

2.11 


8816 

12.72 

0.27 

2.11 

8817 

11.83 

0.25 

2.11 


8817 

6.03 

0.13 

2.11 


175.52 

3.64 

2.07 


8802 

179.91 

3.73 

2.07 

8819 

913.64 

18.77 

2.05 


8819 

948.5 

19.47 

2.05 

EEJ 

22.18 

0.46 

2.07 


8804 

23.94 

0.5 

2.07 


32.86 

0.66 

2.02 


8805 

33.95 

0.69 

2.02 

8810 

63.05 

1.31 

2.07 


8810 

67.55 

1.4 

2.07 


Table A-VIILS Venturi Flowrates and Bias Limits for test TTB045 




)0 




14 



ttb044/109 

rn 

■ 

W 

Ib/sec 

Bw 

Ib/sec 

Bw (%) 


W 

Ib/sec 

Bw 

Ib/sec 

Bw (%) 


W 

Ib/sec 

Bw 

Ib/sec 

R31 

8801 

27.13 

0.56 

2.05 


8801 

28.42 

0.58 

2.05 


8801 

30.1 

0.62 

2.04 

8818 

71.88 

2.19 

3.04 


8818 

74.13 

2.25 

3.04 


8818 

76.91 

2.34 

3.04 

8815 

13.05" 

0.27 

2.07 


8815 

13.44 


2.07 


8815 

14.64 

0.3 

2.07 

8816 

12.52 

0.26 

2.11 


8816 

13.01 

0.27 

2.11 


8816 

14.15 

0.3 

2.11 

8817 

14.62 

0.31 

2.1 


8817 

14.32 

■ESI 

2.1 


8817 

16.21 

0.34 

2.11 

EES 

176.13 

" 3.65 

2.07 


8802 

181.06 

3.75 

2.07 


8802 

186.5 

3:86 

2.07 

8819 

891.37 

18.33 

2.06 


8819 

925.54 

19 

2.05 


8819 

970 

19.89 

2.05 

8804 

25.55 

0.53 

2.07 


8804 

27.19 

0.56 

2.07 


8804 

29.19 

0.6 

2.07 

EES 

34.61 

0.7 

2.02 


8805 

36.12 

0.73 

2.02 


8805 

37.95 

0.77 

2.02 

1 881 0 

64.48 

1.34 

2.07 


8810 

69 

1.43 

2.07 


8810 

74.61 

1.55 

2.07 


Table A-VHL6 Venturi Flowrates and Bias Limits for test TTB044 
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ttb043/100 

W 

Ib/sec 

Bw 

Ib/sec 

27.94 

0.57 

71.92 

2.19 

12.18 

0.25 

12.59 

0.27 

7.18 

0.16 

176.47 

3.65 

896.57 

18.43 

25.67 

0.53 

34.97 

0.71 

61.85 

1.28 


ttb043/104 

W 

Ib/sec 

Bw 

Ib/sec 

29.27 

0.6 

74.19 

2.26 

12.74 

0.27 

13.29 

0.28 

7.27 

0.16 

182.08 

3.77 

932.3" 

19.13 

27.38 

0.57 

36.22 

0.73 

66.59 

1.38 


Bw (%) 


ttb043/1 09 

W 

Bw 

Ib/sec 

Ib/sec 

31 .05 

0.63 

76.87 

2.34 

13.89 

0.29 

14.26 

0.3 

8.37 

0.18 

185.5 

3.84 

973.1 

19.95 

29.22 

0.6 

38.51 

0.78 

72.6 

1.51 


Table A-VIIL7 Venturi Flowrates and Bias Limits for test TTB043 



Table A-VHL7 Venturi Flowrates and Bias Limits for test TTB042 


8801 

26.15 

0.54 

2.05 

8818 

68.57 

2.08 

3.04 

8815 

13.17 

0.27 

2.07 

8816 

13.11 

0.28 

2.11 

8817 

14.59 

0.31 

2.11 

8802 

172.59 

3.57 

2.07 

8819 

898.76 

18.46 

2.05 

8804 

25.67 

■ESI 

2.07 

8805 

30.98 

0.63 

2.03 

8810 

62.11 

1.29 

2.07 


8801 

27.37 

0.56 

2.05 

8818 

72.12 

2.19 

3.04 

8815 

13.39 

0.28 

2.07 

8816 

13.33 

0.28 

2.11 

8817 

14.96 

0.31 

2.11 

8802 

177.16 

3.67 

2.07 

8819 

933.2 

19.14 

2.05 

8804 

27.38 

0.57 

2.07 

8805 

32.26 

0.65 

2.02 

8810 

65.58 

1.36 

2.07 


8801 

28.58 

0.59 

2.05 

8818 

75.79 

2.3 

3.04 

8815 

13 

0.27 

2.07 

8816 

12.99 

0.27 

2.11 

8817 

15.35 

0.32 

2.11 

8802 

183.9 

3.81 

2.07 

8819 

972.3 

19.93 

2.05 

8804 

29.83 

0.62 

2.07 

8805 

33.19 

0.67 

2.02 

8810 

69.96 

1.45 

2.07 



27.95 

0.57 

2.04 

65.15 

1.98 

3.04 

14.75 

0.31 

2.07 

14.3 

0.3 

2.11 

16.07 

0.34 

2.11 

170.03 

3.52 

2.07 

920.51 

18.9 

2.05 

20.96 

0.43 

2.07 

34.23 

0.69 

2.03 

66.52 

1.38 

2.07 



28.61 

0.58 

2.04 

73.7 

2.24 

3.04 

13.34 

0.28 

2.07 

13 

0.27 

2.11 

15.22 

0.32 

2.11 

175 

3.62 

2.07 

953.27 

19.55 

2.05 

22.5 

0.47 

2.07 

35.84 

0.72 

2.02 

70.23 

1.46 

2.07 



30.44 

0.62 

2.04 

76.47 

2.33 

3.04 

13.98 

0.29 

2.07 

13.71" 

0.29 

2.11 

16 

0.34 

2.11 

180.8 

3.75 

2.07 

996 

20.42 

2.05 

25.27 

0.52 

2.07 

36.69 

0.74 

2.02 

76.29 

1.58 

2.07 


Table A-VIIL8 Venturi Flowrates and Bias Limits for test TTB040 
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Appendix IX 

Test-to-Test Plots of Venturi Mass Flowrate 
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Plot of Total Flow Rate per Test 
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72 



56 




(Ib/sec) 



5 ' 















(Ib/sec) 


HPOP 



6 


vs Test # 
















Appendix X 

Example Output from Venturi Uncertainty Program 

The output provided is from test TTB050 for 3 -5 second time slices. The 
definitions of the variables listed are as follows: 


PID 

PI 

DP 

T 

BPl 

BDP 

BTl 

BCD 

BALP 

BDl 

BD2 

PW 

W 

UW 

%UW 

DWDP1= 


DWDDP= 


DWDT= 

DWDCD= 


DWALP= 


DWDD1= 


DWDD2= 


Measurement identification 

Inlet pressure 

Differential pressure 

Temperature 

Inlet pressure bias limit 

Differential pressure bias limit 

Temperature bias limit 

Discharge coefficient bias limit 

Thermal expansion coefficient bias limit 

Venturi throat diameter bias limit 

Venturi inlet diameter bias limit 

Flowrate precision limit 

Flowrate 

Uncertainty in flowrate 
Percent uncertainty in flowrate 


dW 


8P 1 

BWP= 

aw 


8AP 

BWDP= 

aw 


8T 

BWDT= 

aw 


dC D 

aw 

BWCD= 

da 


aw 

BWALP= 

am 


aw 

BWD1= 

dDl 

- - 


BWD2= 


f 


8W 

dP 

aw 

KaAP 

' ' aw 

^aAT 

f aw 


b b 


v 


B 


hp 


B, 


AT 


B r 


0C D 
dW_ n 

da 


aw_ 
. am 

aw 

\8D2 


B 


D\ 


B 1 


D 2 
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TEST#: 8010050. ENGINE #: 3001. #OFPIDS: 43 # OF SLICES: 3 


SLNUM >*♦*** SLICE START TIME > 16.0 SLICE END TIME > 21.0 

ILTTB= 1 DTTB(ILTTB)= 3998.69 

ILTTB = 2 DTTB(ILTTB)= 256.49 

ILTTB = 3 DTTB(ILTTB)= 473.62 

ILTTB = 4 DTTB(ILTTB)= 5223.84 

ILTTB = 5 DTTB(ILTTB) = 386.29 

ILTTB = 6 DTTB(ILTTB)= 92.91 

ILTTB = 7 DTTB(ILTTB)= 5383.70 

ILTTB = 8 DTTB(ILTTB)= 119.84 

ILTTB- 9 DTTB(ILTTB) - 91.82 

ILTTB = 10 DTTB(ILTTB) = 5385.07 

ILTTB = II DTTB(ILTTB) = 117.56 

ILTTB = 12 DTTB(ILTTB) = 92.69 

ILTTB = 13 DTTB(ILTTB) = 5392.00 

ILTTB = 14 DTTB(ILTTB) = 121.34 

ILTTB = 15 DTTB(ILTTB) = 90.92 

ILTTB = 16 DTTB(ILTTB)= 3796.21 

ILTTB = 17 DTTB(ILTTB)= 207.17 

ILTTB = 18 DTTB(ILTTB)= 188.41 

ILTTB = 19 DTTB(ILTTB)= 3717.76 

ILTTB = 20 DTTB(ILTTB)= 245.44 

ILTTB =21 DTTB(ILTTB) = 190.50 

ILTTB = 22 DTTB(ILTTB)= 6621.04 

ILTTB = 23 DTTB(ILTTB) = 134.31 

ILTTB = 24 DTTB(ILTTB)= 215.41 

ILTTB = 25 DTTB(ILTTB) = 4887.78 

ILTTB = 26 DTTB(ILTTB) = 261.37 

ILTTB = 27 DTTB(ILTTB)= 280.75 

ILTTB = 28 DTTB(ILTTB)= 6653.16 

ILTTB = 29 DTTB(ILTTB)= 97.34 

ILTTB = 30 DTTB(ILTTB)= 204.71 

LPFT INLET FLOWRATE - P(8801) > 26.61 LB/S 
CCV INLET FLOWRATE - P(8818) > 71.41 LB/S 
NOZ CLNT FLOWRATE - P(8815) > 13.08 LB/S 
NOZ CLNT FLOWRATE - P(8816) > 12.78 LB/S 
NOZ CLNT FLOWRATE - P(8817) > 12.98 LB/S 
LPOT INLET FLOWRATE -P(8802)> 183.29 LB/S 
HPOP DISCH FLOWRATE - P(8819) > 91039 LB/S 
OPB LOX FLOWRATE -P(8804)> 24.05 LB/S 
OPB FUEL FLOWRATE - W(8805) > 3439 LB/S 
FPB LOX FLOWRATE -W(8810)> 57.47 LB/S 
PID 8801 LPFT INLET 

VENTURI <20> 

PI 3998.69 DP 256.49 T 473.62 

BP1 37.9621 BDP 2.0717 BT1 1.8841 

BCD 0.01992 BALP .417E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3128E-02 DWDDP .4582B01 DWDT -.243E^)1 
DWDCD .2624E+02 DWALP ^.293E+04 DWDD1 -.120E-K)2 DWDD2 .5076E+02 
BWP .1410E-01 BWDP .9010B02 BWDT .2100E^02 
BWCD .2731E+00 BWALP .1490B05 BWD1 .3589E4M BWD2 .6442E-03 
W 26.6139 UW 0.54654 %UW 2.0536 
********* ***** ******** ***** 

PID 8818 CCV INLET 

VENTURI <397> 

PI 5223.84 DP 38639 T 92.91 

BP1 52.2384 BDP 3.8629 BT1 0.9291 

BCD 0.03063 BALP .548E4)6 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .2260E-02 DWDDP .8257E-01 DWDT -.114E+00 
DWDCD .6990E+02 DWALP -.612E+05 DWDD1 -.156E+02 DWDD2 .1069E-K)3 
BWP .1394E-01 BWDP .1017E+00 BWDT .1127E-01 
BWCD .4584E+01 BWALP .1126E-02 BWD1 .6050E-04 BWD2 .2855E-02 
W 71.4057 UW 2.17125 %UW 3.0407 
***** ********* ***** ******** ***** 
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PID 8815 NOZ CLNT 

VENTURI <3541> 

PI 5223.84 DP 386.29 T 92.91 

BP1 53.8370 BDP 1.1984 BT1 0.9182 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3235E-03 DWDDP .5311E-01 DWDT -.217E-01 

DWDCD .1287E+02 DWALP -.110E+05 DWDD1 -.328E+01 DWDM .3210E+02 
BWP .3033E-03 BWDP .4051E-02 BWDT .3978E-03 
BWCD .6842E-01 BWALP .3607E-04 BWD1 .2688E-05 BWM .2575E-03 
W 13.0809 UW 0.27095 %UW 2.0713 

***** ********* ***** ******** ***** 


PID 8816 NOZ CLNT 

VENTURI <3542> 

PI 5385.07 DP 117.56 T 92.69 

BP1 53.8507 BDP 1.1756 BT1 0.9269 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BM 0.00050 

PW 0.0000 

DWDP1 .3153E-03 DWDDP .5293E-01 DWDT -.215E-01 
DWDCD .1283E+02 DWALP -.110E+05 DWDD1 -.325E+01 DWDM .3133E+02 
BWP .2882E-03 BWDP .3872E-02 BWDT .3977E-03 
BWCD .6799E-01 BWALP .3607E-O4 BWD1 .2639E-05 BWM .2453E-03 
W 12.7825 UW 0.26978 %UW 2.1106 


PID 8817 NOZ CLNT 

VENTURI <3543> 

PI 5392.00 DP 121.34 T 90.92 

BP1 53.9200 BDP 1.2134 BT1 0.9092 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BM 0.00050 

PW 0.0000 

DWDPl 3074E-03 DWDDP .5182E-01 DWDT -.220E-01 
DWDCD .1299E+02 DWALP -.134E+05 DWDD1 -331E+0I DWDM .3190E+02 
BWP .2747E-03 BWDP .3954E-02 BWDT .3993E-03 
BWCD .6978E-01 BWALP .5388E-04 BWD1 .2740E-05 BWM .2544E-03 
W 12.9831 UW 0.27326 %UW 2.1047 

***** ********* ***** •****••• ***** 


PID 8802 LPOT INLET 

VENTURI <139> 

PI 3796.21 DP 207.17 T 188.41 

BP1 37.9621 BDP 2.0717 BT1 1.8841 

BCD 0.02066 BALP .6S1E-06 BD1 0.00050 BM 0.00050 

PW 0.0000 

DWDPl .1137E-02 DWDDP .4420E400 DWDT -.183E+00 
DWDCD .1774E+03 DWALP -.124E+06 DWDD1 -.446E+02 DWDM .2986E+03 
BWP .1863E-02 BWDP .8386E-H30 BWDT .1I95E-K)0 
BWCD .1343E+02 BWALP .6565E-02 BWD1 .4981E-03 BWM .2229E-01 
W 183.2910 UW 3.79655 %UW 2.0713 
***** ********* ***** ******** ***** 


PID 8819 HPOP DISCH 

VENTURI <426> 

PI 3717.76 DP 245.44 T 190.50 

BP1 37.1776 BDP 1.9050 BT1 2.4544 

BCD 0.01964 BALP .549E-06 BD1 0.00050 BD2 0.00050 

PW 0 0000 

DWDPl .5711E-02 DWDDP .1853E-K)1 DWDT -.930E-KJ0 
DWDCD .9270E+03 DWALP -.124E+06 DWDD1 -.326E+03 DWDD2 .9771E+03 
BWP .4508E^)1 BWDP .1247E+02 BWDT .5211E+01 
BWCD 3315E+Q3 BWALP .1167E+QG BWDl .26S8BOI BWD2 .23X7E+Q0 
W 9103882 UW 18.69283 %UW 2.0533 
***** ********* ***** ******** ***** 
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PID $804 OPB LOX 

VENTURI <268> 

PI 6621.04 DP 134.31 T 215.41 

BP1 66.2104 BDP 1.3431 BT1 2.1541 

BCD 0.02066 BALP 0.671E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .1475E-03 DWDDP .8961E-01 DWDT -.218E-01 
DWDCD .2318E+02 DWALP -.136E+05 DWDD1 -.635E401 DWDD2 .8418E-K)2 
BWP .9539E-04 BWDP .1448E-01 BWDT .2207E-02 
BWCD .2294E+00 BWALP .8374E-04 BWD1 .1009E-04 BWD2 .1772E-02 
W 24.0457 UW 0.49778 %UW 2.0702 

***** ********* ***** ******** ***** 


PID 8805 OPB FUEL 
VENTURI <27 1> 

PI 4887.78 DP 261.37 T 280.75 

BP1 48.8778 BDP 2.6137 BT1 2.8075 

BCD 0.01992 BALP 0.323E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .2638E-02 DWDDP .1117E-03 DWDT -.480E-01 
DWDCD .3366E+02 DWALP -.170E+05 DWDD1 -.135E+02 DWDD2 .6442E+02 
BWP .1662E-01 BWDP .8523E-07 BWDT .1812E-01 
BWCD .4495E+00 BWALP .3017E-04 BWD1 .4587E-04 BWD2 .1037E-02 
W 343923 UW 0.69639 %UW 2.0248 
***** ********* ***** ******** ***** 


PID 8810 FPB LOX 
VENTURI <296> 

PI 6653.16 DP 97.34 T 204.71 

BP1 66.5316 BDP 0.9734 BT1 2.0471 

BCD 0.02052 BALP .671E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 3165E-03 DWDDP .2953E+00 DWDT -.508E-01 
DWDCD .5605E+02 DWALP -.371E+05 DWDD1 -.584E+01 DWDD2 .1143E+03 
BWP .4435E-03 BWDP .8260E-01 BWDT .1082E-01 
BWCD .1323E+01 BWALP .6186E-03 BWD1 .8539E-05 BWD2 .3264E-02 
W 57.4704 UW 1.19193 %UW 2.0740 
***** ********* ***** ******** ***** 


SLNUM >♦*♦** SUCE START TIME > 16.0 SLICE END TIME > 21.0 

8801 26.61 0.55 2.05 

8818 71.41 2.17 3.04 

8815 13.08 0.27 2.07 

8816 12.78 0.27 2.11 

8817 12.98 0.27 2.10 

8802 183.29 3.80 2.07 

8819 910.39 18.69 2.05 

8804 24.05 0.50 2.07 

8805 34.39 0.70 2.02 
8810 57.47 1.19 2.07 
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TEST#: 8010050. ENGINE #: 3001. #OFPIDS: 43 # OF SLICES: 3 


SLNUM >♦**** SLICE START TIME > 25.0 SUCE END TIME > 30.0 

ILTTB- 1 DTTB(ILTTB)= 4185.40 

ILTTB= 2 DTTB(ILTTB)= 271.89 

ILTTB = 3 DTTB(ILTTB)= 468.97 

ILTTB- 4 DTTB(ILTTB)= 5490.61 

ILTTB = 5 DTTB(ILTTB)= 410.88 

ILTTB = 6 DTTB(ILTTB)= 94.94 

ILTTB- 7 DTTB(ILTTB)= 5657.44 

ILTTB- 8 DTTB(ILTTB)= 132.33 

ILTTB = 9 DTTB(ILTTB) = 93.74 

ILTTB = 10 DTTB(ILTTB)= 5660.28 

ILTTB- 11 DTTB(ILTTB) = 130.22 

ILTTB = 12 DTTB(ILTTB) = 94.90 

ILTTB- 13 DTTB(ILTTB)= 5666.14 

ILTTB- 14 DTTB(ILTTB)= 132.12 

ILTTB- 15 DTTB(ILTTB) = 92.95 

ILTTB- 16 DTTB(ILTTB)= 3981.13 

ILTTB- 17 DTTB(ILTTB) = 218.58 

ILTTB- 18 DTTB(ILTTB)= 189.46 

ILTTB- 19 DTTB(ILTTB) = 3895.58 

ILTTB- 20 DTTB(ILTTB) = 265.92 

ILTTB = 21 DTTB(ILTTB) = 191.59 

ILTTB- 22 DTTB(ILTTB)= 6917.97 

ILTTB- 23 DTTB(ILTTB) = 157.95 

ILTTB- 24 DTTB(ILTTB)= 216.81 

ILTTB - 25 DTTB(ILTTB) - 5138.72 

ILTTB- 26 DTTB(ILTTB) = 272.94 

ILTTB- 27 DTTB(ILTTB)= 281.84 

ILTTB- 28 DTTB(ILTTB)= 6954.25 

ILTTB- 29 DTTB(ILTTB) = 110.65 

ILTTB- 30 DTTB(ILTTB) = 206.33 

LPFT INLET FLOWRATE * P(8801) > 28.02 LB/S 
CCV INLET FLOWRATE - P(8818) > 73.78 LB/S 
NOZ CLNT FLOWRATE - P(8815) > 13.77 LB/S 
NOZ CLNT FLOWRATE - P(8816) > 13.48 LB/S 
NOZ CLNT FLOWRATE - P(8817) > 13.58 LB/S 
LPOT INLET FLOWRATE - P(8802) > 188.29 LB/S 
HPOP DISCH FLOWRATE - P(8819) > 947.63 LB/S 
OPBLOX FLOWRATE -P(8804)> 26.09 LB/S 
OPB FUEL FLOWRATE - W(8805) > 35.68 LB/S 
FPB LOX FLOWRATE- W(8810)> 61.29 LB/S 
PID 8801 LPFT INLET 
VENTURI <20> 

PI 4185.40 DP 271.89 T 468.97 

BP1 39.8112 BDP 2.1858 BT1 1.8946 

BCD 0.01992 BALP .417E.06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDPl .3124E42 DWDDP .4541E4)1 DWDT -.256E^)1 
DWDCD .2763E+02 DWALP -.329E+04 DWDD1 -.126E-K)2 DWDD2 .5341E+02 
BWP ,1547Er01 BWDP .9854E4)2 BWDT .2347E-02 
BWCD 3030E+00 BWALP .1886E-05 BWD1 .3970E-04 BWD2 .7131E4)3 
W 28.0166 UW 0.57543 %UW 2.0539 

***** •****•*•* ***** ******** ***** 

PID 8818 CCV INLET 

VENTURI < 397 > 

PI 5490.61 DP 410.88 T 94.94 

BP1 54.9061 BDP 4.1088 BT1 0.9494 

BCD 0.03063 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDPl .2265E-02 DWDDP .7986E-01 DWDT -.114E400 
DWDCD .7222E+02 DWALP -.668E+05 DWDD1 -.160E+02 DWDD2 .1104E+03 
BWP .1546E-01 BWDP .1077E+00 BWDT .1163E-01 
BWCD .4893E+01 BWALP .1340E-02 BWD1 .6425E-04 BWD2 .3044E-02 
W 73.7762 UW 2.24304 %UW 3.0403 

***** ********* ***** ******** ***** 
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PED 8815 NOZ CLNT 

VENTURI <3541> 

PI 5490.61 DP 410.88 T 94.94 

BP1 56.5744 BDP 1.3233 BT1 0.9374 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3287E-03 DWDDP .5051E-01 DWDT -.223E-01 
DWDCD .1355E+02 DWALP -.117E+05 DWDD1 -.345E+01 DWDD2 .3378E+02 
BWP .3459E-03 BWDP .4468E-02 BWDT .4359E-03 
BWCD .7583 E*01 BWALP .4079E-04 BWDl .2971E-05 BWD2 .2854E-03 
W 13.7743 UW 0.28523 %UW 2.0707 

***** ********* ***** ******** ***** 


PED 8816 NOZ CLNT 

VENTURI <3542> 

PI 5660.28 DP 130.22 T 94.90 

BP1 56.6028 BDP 1.3022 BT1 0.9490 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3267E-03 DWDDP .5005E-01 DWDT -.219E-01 
DWDCD .1352E+02 DWALP -.115E+05 DWDD1 -.343E+01 DWDD2 .3300E-K)2 
BWP .3420E-03 BWDP .4247E-02 BWDT .4308E-03 
BWCD .7555E-01 BWALP 3958E-04 BWDl .2937E-05 BWD2 .2722E-03 
W 13.4754 UW 0.28430 %UW 2.1098 

***** ********* ***** ******** ***** 


PID 8817 NOZ CLNT 
VENTURI <3543> 

PI 5666.14 DP 132.12 T 92.95 

BP1 56.6614 BDP 1.3212 BT1 0.9295 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 3195E-03 DWDDP .4951E-01 DWDT -.218E-01 
DWDCD .1358E+02 DWALP -.137E+05 DWDD1 -346E+01 DWDD2 .3335E+02 
BWP 3278E-03 BWDP .4278E-02 BWDT .4123E-03 
BWCD .7624E-01 BWALP .5671E-04 BWDl .2985E-05 BWD2 .2781E-03 
W 13.5752 UW 0.28555 %UW 2.1034 
***** ********* ***** ******** ***** 


PID 8802 LPOT INLET 
VENTURI <139> 

PI 3981.13 DP 218.58 T 189.46 

BP1 39.8112 BDP 2.1858 BT1 1.8946 

BCD 0.02066 BALP .651E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .1165E-02 DWDDP .4306E+00 DWDT -.188E+00 
DWDCD .1822E-K)3 DWALP -.128E+06 DWDD1 -.459E402 DWDD2 .3068E+03 
BWP .2150E-02 BWDP .8860E-K)0 BWDT .1264E-K)0 
BWCD .1417E-H)2 BWALP .6915E-02 BWDl .5264E-03 BWD2 .2353E-01 
W 188.2875 UW 3.90044 %UW 2.0715 
***** ********* ***** ******** ***** 


PID 8819 HPOP DISCH 

VENTURI <426> 

PI 3895.58 DP 265.92 T 191.59 

BP1 38.9558 BDP 1.9159 BT1 2.6592 

BCD 0.01964 BALP .549E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .5951E-02 DWDDP .1781E-K)1 DWDT -.958E+00 
DWDCD .9650E+03 DWALP -.128E+06 DWDD1 -.339E+03 DWDD2 .1017E+04 
BWP .5375E-01 BWDP .1165E+02 BWDT .6488E+01 
BWCD 3592E+03 BWALP .1340E+00 BWDl .2880E-01 BWD2 .2586E+00 
W 947.6333 UW 19.43365 %UW 2.0508 
***** ********* ***** ******** ***** 
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PID 8804 OPB LOX 
VENTURI <268> 

PI 6917.97 DP 157.95 T 216.81 

BP1 69.1797 BDP 1.5795 BT1 2.1681 

BCD 0.02066 BALP 0.671E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .1544E*03 DWDDP .8257E-01 DWDT -.234E-01 

DWDCD .2516E+02 DWALP -.159E+05 DWDDI -.689E+01 DWDD2 .9133E+02 
BWP .1141E-03 BWDP .1701E-01 BWDT .2565E-02 
BWCD .2701E+00 BWALP .1140E-03 BWD1 .1188E-04 BWD2 .2085E-02 
W 26.0900 UW 0.54008 %UW 2.0701 

***** ********* ***** ******** ***** 


PID 8805 OPB FUEL 
VENTURI <27 1> 

PI 5138.72 DP 272.94 T 281.84 

BP1 51.3872 BDP 2.7294 BT1 2.8184 

BCD 0.01992 BALP 0.323E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .2554E-02 DWDDP .1151E-03 DWDT -.485E-01 
DWDCD .3493E+02 DWALP -.175E405 DWDDI -.141E+02 DWDD2 .6684E+02 
BWP .1722E-01 BWDP .9864E-07 BWDT .1869E-01 
BWCD .4841E+00 BWALP .3187E-04 BWD1 .4949E-04 BWD2 .1117E-02 
W 35.6798 UW 0.72162 %UW 2.0225 
***** ********* ***** ******** ***** 


PID 8810 FPB LOX 

VENTURI <296> 

PI 6954.25 DP 110.65 T 206.33 

BP1 69.5425 BDP 1.1065 BT1 2.0633 

BCD 0.02052 BALP .671E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3290E-03 DWDDP .2769E+00 DWDT -.535E-01 
DWDCD .5971E+02 DWALP -396E+05 DWDDI -.623E-KJ1 DWDD2 .1218E+03 
BWP .5236E-03 BWDP .9390E-01 BWDT .1219E01 
BWCD .1501E+01 BWALP .7049E-03 BWD1 .9708E^)5 BWD2 .3709E-02 
W 61.2905 UW 136982 %UW 2.0718 


SLNUM >•***• SUCE START TIME > 25.0 SUCE END TIME > 30.0 


8801 

28.02 

0.58 

2.05 

8818 

73.78 

2.24 

3.04 

8815 

13.77 

0.29 

2.07 

8816 

13.48 

0.28 

2.11 

8817 

13.58 

0.29 

2.10 

8802 

188.29 

3.90 

2.07 

8819 

947.63 

19.43 

2.05 

8804 

26.09 

0.54 

2.07 

8805 

35.68 

0.72 

2.02 

8810 

61.29 

1.27 

2.07 
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TEST#: 8010050. ENGINE #: 3001. #OFPIDS: 43 # OF SLICES: 3 

SLNUM >**•** SUCE START TIME > 34.0 SLICE END TIME > 38.0 

ILTTB - 1 DTTB(ILTTB)= 4426.71 

ILTTB = 2 DTTB(ILTTB)= 291.82 

ILTTB * 3 DTTB(ILTTB) = 464.59 

ILTTB = 4 DTTB(ILTTB) = 5838.52 

ILTTB = 5 DTTB(ILTTB) = 440.25 

ILTTB = 6 DTTB(ILTTB) = 97.59 

ILTTB = 7 DTTB(ILTTB) = 6016.47 

ILTTB = 8 DTTB(ILTTB)= 149.35 

ILTTB = 9 DTTB(ILTTB) = 96.34 

ILTTB « 10 DTTB(ILTTB) = 6020.90 

ILTTB = 11 DTTB(ILTTB)= 147.22 

ILTTB = 12 DTTB(ILTTB)= 97.43 

ILTTB = 13 DTTB(ILTTB)= 6026.65 

ILTTB « 14 DTTB(ILTTB) = 146.61 

ILTTB = 15 DTTB(ILTTB)= 95.76 

ILTTB - 16 DTTB(ILTTB)= 4233.10 

ILTTB ~ 17 DTTB(ILTTB)= 235.31 

ILTTB = 18 DTTB(ILTTB) = 190.95 

ILTTB = 19 DTTB(ILTTB)= 4136.17 

ILTTB = 20 DTTB(ILTTB)= 292.10 

ILTTB = 21 DTTB(ILTTB)= 193.23 

ILTTB = 22 DTTB(ILTTB) = 7323.87 

ILTTB = 23 DTTB(ILTTB)= 189.11 

ILTTB = 24 DTTB(ILTTB) = 219.42 

ILTTB - 25 DTTB(ILTTB)= 5467.02 

ILTTB = 26 DTTB(ILTTB)= 283.85 

ILTTB = 27 DTTB(ILTTB)= 284.05 

ILTTB = 28 DTTB(ILTTB)= 7369.33 

ILTTB = 29 DTTB(ILTTB) = 12839 

ILTTB = 30 DTTB(ILTTB)= 208.63 

LPFT INLET FLOWRATE - P(8801) > 29.78 LB/S 
CCV INLET FLOWRATE - P(8818) > 76.57 LB/S 
NOZ CLNT FLOWRATE - P(8815) > 14.67 LB/S 
NOZ CLNT FLOWRATE - P(8816) > 1437 LB/S 
NOZ CLNT FLOWRATE - P(8817) > 1433 LB/S 
LPOT INLET FLOWRATE - P(8802) > 19537 LB/S 
HPOP DISCH FLOWRATE - P(8819) > 993.06 LB/S 
OPB LOX FLOWRATE - P(8804) > 28.55 LB/S 
OPB FUEL FLOWRATE - W(8805) > 37.04 LB/S 
FPB LOX FLOWRATE -W(8810)> 66.03 LB/S 
PID 8801 LPFT INLET 

VENTURI <20> 

PI 4426.71 DP 291.82 T 464.59 

BP1 423310 BDP 23531 BT1 1.9095 

BCD 0.01992 BALP .417E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 3099E-02 DWDDP .4479E-01 DWDT -.271E-01 
DWDCD .2938E+02 DWALP -.403E+04 DWDD1 -.134E+02 DWDD2 .5674E+02 
BWP .1721E-01 BWDP .1111E-01 BWDT .2682E-02 
BWCD 3425E+00 BWALP .2817E-05 BWD1 .4467E-04 BWD2 .8049E-03 
W 29.7824 UW 0.61157 %UW 2.0534 

***** ********* ***** ******** ***** 

PID 8818 CCV INLET 

VENTURI <397> 

PI 5838.52 DP 440.25 T 97.59 

BP1 583852 BDP 4.4025 BT1 0.9759 

BCD 0.03063 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .2265E-02 DWDDP ,7734E^)1 DWDT -.113E+00 
DWDCD .7498E+02 DWALP -.640E+05 DWDD1 -.166E+02 DWDD2 .1145E+03 
BWP .1749E-01 BWDP .1159E-KK) BWDT .1213E-01 
BWCD .5275E+01 BWALP .1231E-02 BWD1 .6884E-04 BWD2 .3277E02 
W 76.5676 UW 232890 %UW 3.0416 

***** ********* ***** ******** ***** 
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PID 8815 NOZ CLNT 
VENTURI <3541> 

PI 5838.52 DP 440.25 T 97.59 

BP1 60.1647 BDP 1.4935 BT1 0.9634 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DW DPI .3286E-03 DWDDP .4741E-01 DWDT -.233E-01 

DWDCD .1443E+02 DWALP -.148E+05 DWDD1 -.366E+01 DWDD2 .3596E+02 
BWP 3909E-03 BWDP .5014E-02 BWDT .5045E-03 
BWCD .8602E-01 BWALP .6565E-04 BWD1 .3343E-05 BWD2 .3233E-03 
W 14.6686 UW 030373 %UW 2.0707 

***** ********* ***** ******** ***** 


PID 8816 NOZ CLNT 

VENTURI <3542> 

PI 6020.90 DP 147.22 T 97.43 

BP1 60.2090 BDP 1.4722 BT1 0.9743 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3289E-03 DWDDP .4701E-01 DWDT -.222E-01 
DWDCD .1441E+02 DWALP -.124E+05 DWDD1 -.364E+01 DWDD2 .3517E-K)2 
BWP .3921E-03 BWDP .4789E-02 BWDT .4692E-03 
BWCD .8588E-01 BWALP .4581E-04 BWD1 .3312E-05 BWD2 .3093E-03 
W 14.3658 UW 0.30303 %UW 2.1094 

***** ********* ***** ******** ***** 


PID 8817 NOZ CLNT 

VENTURI <3543> 

PI 6026.65 DP 146.61 T 95.76 

BP1 60.2665 BDP 1.4661 BT1 0.9576 

BCD 0.02033 BALP .548E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3237E-03 DWDDP .4706E-01 DWDT -.222E-01 
DWDCD .1435E+02 DWALP -.124E-KJ5 DWDD1 -.363E401 DWDD2 .3521E-K)2 
BWP 3805E-03 BWDP .4760E-02 BWDT .4507E-03 
BWCD .8512E-01 BWALP .4581E-04 BWD1 3302E-05 BWD2 .3099E-03 
W 143339 UW 0.30167 %UW 2.1046 

***** ********* ***** ******** ***** 


PID 8802 LPOT INLET 

VENTURI <139> 

PI 4233.10 DP 235.31 T 190.95 

BP1 423310 BDP 23531 BT1 1.9095 

BCD 0.02066 BALP .651E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .1179E-02 DWDDP .4150E+00 DWDT -.193E-H)0 
DWDCD .1891E+03 DWALP -.134E+06 DWDD1 -.476E+02 DWDD2 .3183E+03 
BWP .2489E-02 BWDP .9535E+00 BWDT .1353E-K)0 
BWCD .1526E+02 BWALP .7591E-02 BWD1 .5676E-03 BWD2 .2533E-01 
W 1953723 UW 4.04706 %UW 2.0715 
***** ********* ***** ******** ***** 


PID 8819 HPOP DISCH 
VENTURI <426> 

PI 4136.17 DP 292.10 T 193.23 

BP1 41.3617 BDP 1.9323 BT1 2.9210 

BCD 0.01964 BALP .549E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .6197E-02 DWDDP .1701E-H)1 DWDT -.993E+00 

DWDCD .1012E+04 DWALP -.134E+06 DWDD1 -.356E+03 DWDD2 .1066E+04 

BWP .6570E-01 BWDP ,1081E*K)2 BWDT .8408E+01 

BWCD 3947E+03 BWALP .1340E-KK) BWD1 3161E4)1 BWD2 .2841E-KK) 

W 993.0620 UW 2035188 %UW 2.0494 
***** ********* ***** ******** ***** 


72 


PID 8804 OPB LOX 
VENTURI <268> 

PI 7323.87 DP 189.11 T 219.42 

BP1 73.2387 BDP 1.8911 BT1 2.1943 

BCD 0.02066 BALP 0.671E4D6 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .1687E-03 DWDDP .7550E-01 DWDT -.250E-01 

DWDCD .2754E+02 DWALP -.159E+05 DWDD1 -.755E+01 DWDD2 .9995E+02 
BWP .1527E-03 BWDP .2038E01 BWDT .3015E4J2 
BWCD .3237E+00 BWALP .1140E-03 BWD1 .1424E-04 BWD2 .2497E-02 
W 28.5502 UW 0.59121 %UW 2.0708 

***** ********* ***** ******** ***** 


PID 8805 OPB FUEL 
VENTURI <271> 

PI 5467.02 DP 283.85 T 284.05 

BP1 54.6702 BDP 2.8385 BT1 2.8405 

BCD 0.01992 BALP 0.323E-06 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .2442E-02 DWDDP .1226E-03 DWDT -.490E-01 
DWDCD .3627E+02 DWALP -.184E+05 DWDD1 -.146E+02 DWDD2 .6944E+02 
BWP .I782E-01 BWDP .1211E-06 BWDT .1941E-01 
BWCD .5221E+00 BWALP 3541E-04 BWD1 .5352E-04 BWD2 .1206E-02 
W 37.0400 UW 0.74839 %UW 2.0205 
***** ********* ***** *••*••*• ***** 


PID 8810 FPB LOX 
VENTURI <296> 

PI 7369.33 DP 128.39 T 208.63 

BP1 73.6933 BDP 1.2839 BT1 2.0863 

BCD 0.02052 BALP .671E416 BD1 0.00050 BD2 0.00050 

PW 0.0000 

DWDP1 .3540E-03 DWDDP .2572E+00 DWDT -.565E-01 
DWDCD .6434E+02 DWALP -.425E+05 DWDD1 -.675E+01 DWDD2 .1312E+03 
BWP .6805E-03 BWDP .1091E+00 BWDT .1390E-01 
BWCD .I743E+01 BWALP .8142&03 BWD1 .1I38E-04 BWD2 .4305E-02 
W 66.0350 UW 136820 %UW 2.0719 


SLNUM >****» SUCE START TIME > 34.0 SUCE END TIME > 38.0 

8801 29.78 0.61 2.05 

8818 76.57 233 3.04 

8815 14.67 0.30 2.07 

8816 1437 030 2.11 

8817 1433 0.30 2.10 

8802 19537 4.05 Z07 

8819 993.06 2035 2.05 

8804 28.55 0.59 2.07 

8805 37.04 0.75 2.02 
8810 66.03 137 2.07 
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Appendix XI 


Venturi Uncertainty Analysis Computer Program 

This program is written is FORTRAN and resides in the EADS 10 computer system as 
program VENJB2 in the following directory, SIMPSSP.MSFC.VENKB. 
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PROGRAM VENJB 
C 

C THIS PROGRAM HAS BEEN MODIFIED BY K. BROWN TO CALCULATE THE 
C UNCERTAINTIES IN THE VENTURI CALCULATIONS. 

C 

DOUBLE PRECISION DTTB(40),DH 
INTEGER PNSLC,PBTIME,PETIME 

DIMENSION PIDTTB(40),WW( 10),UNW( 10),PER( 1 0),NPID( 1 0) 

C 

NPID(1)=8801 

NPID(2)=8818 

NPID(3)=8815 

NPID(4)=8816 

NPID(5)=8817 

NPID(6)=8802 

NPID(7)=8819 

NPID(8)=8804 

NPID(9)=8805 

NPED(10)=8810 

C 

C READ IN TEST DATA 
C 

READ(36,100) TESTJENGNUMJPIDNJSLCN 
READ(36,1 1 0) PN SLC^PBTIMEyPETIME, 

& (PIDTTB(ITTB), ITTB=1,40) 

1 READ(36,120JENI>=999) NSLC,BTIME3TIME, 

& (DTTB(ITTB)4TTB= 1 ,40) 

100 FORMAT(yTllJ8.0,T32,F5.0,T51,D,T70J3/) 

110 FORMAT(l 5(1X44)) 

120 F0RMAT(5(1X£14.8E2)) 

C 

WRITE (6,150) TESTJZNGNUMJPIDNJSLCN 
WRITE (6, 160) P NSLCJ>BTIMEJ>ETIME, 

& (PIDTTB(ITTByTTB== 1 ,40) 

WRITE(55,150) TESTJENGNUMJPIDN,ISLCN 
WRITE ( 6,170) NSLC3TIME3TIME, 

& (DTTB(TTTB)4TTB= 1 ,40) 

150 FORMAT(/2X,TEST#: , 4 ? 8.03X, , ENGINE#: , ^5.03X, 

& *# OF PIDS: '433X,*# OF SUCES: ’43/) 

160 FORMAT(15(lX44)) 

170 FORMAT(5(lX£14.8E2)) 

C DO7Q0,ILTTB=U0 

C700 READ(10,702) DTTB(ILTTB) 

C READ(10,*) 

C801 FORMATO 
C702 FORMAT(4X4 J 9.2) 

C 

WRITE(55,205) NSLC,BTIME£TIME 
205 FORMAT(2X,*SLNUM >'45,2X,*SLICE START TIME > VF5.1, 

$ 2X/SLICE END TIME > ^5.1) 

DO209,ILTTB=U0 
WR1TE(55,210) ILTTB43TTB(ILTTB) 

210 FORMAT(2X,UTTB^45,4XT>TTB(ILTTB) = , 4 : 9.2) 

209 CONTINUE 
WRITE(55,*) 

C 

H=1.0 

CALL FLOWOR(20,DTTB(l)4DTTB(2)43TTB(3)^4iDH,W,WG4tO,VDPT) 
WRITE(55,320)WG 

320 FORMAT(2X,LPFT INLET FLOWRATE - P(8801) >'J9.2,' IBIS') 

CALL FLOWOR(3974)TTB(4)4)TIB(5)43TTB(6X94iDH,W,WG4lO,VDPT) 
WRTTE(55330)WG 

330 FORMAT(2X,'CCV INLET FLOWRATE - P(8818) IBIS') 

CALL FLOWOR(354143TIB(7)43TTB(8)43TTB(9)^4LDH,W,WG4tO,VDPT) 
WRITE(55340)WG 

340 FORMAT(2X,'NOZ CLNT FLOWRATE - P(8815) >'?9.2; IBIS') 

CALL FLOWOR(35424)TTB( 10XDTTB(1 l)4>TTB(12X94LDH,W,WG4lO,VDPT) 
WRTTE(55,350)WG 

350 FORMAT(2X, , NOZ CLNT FLOWRATE - P(881 6) >'f9.2,' IBIS') 

CALL FLOWOR(35434)TTB(13XDTTB(14)43TTB(15X94LDH,W,WG,RO,VDPT) 
WRITE(55,360)WG 
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360 FORMAT(2X,'NOZ CLNT FLOWRATE - P(8817) >',F9.2,' LB/S') 

CALL FLOWOR( 1 39,DTTB(1 6),DTTB( 1 7),DTTB( 1 8),9,H,DH, W.WG.RO, VDPT) 
WRITE(55,370) W 

370 FORMAT(2X/LPOT INLET FLOWRATE - P(8802) >',F9.2,' LB/S') 

CALLFLOWOR(426J3TTB(19XDTTB(20),DTTB(21X9,HJDH,W,WG,RO,VDPT) 
WRITE(55380) W 

380 FORMAT(2X,HPOP DISCH FLOWRATE - P(88 19) >',F9.2,' LB/S) 

CALLFLOWOR(268 > DTTB(22XDTTB(23XDTTB(24),9 > H,DH,W,WG > RO,VDPT) 

WRITE(55,390)W 

390 FORMAT(2X,’OPB LOX FLOWRATE - P(8804) >',F9.2,’ LB/S') 

CALL FLOWOR(27 1 ,DTTB(25),DTTB(26),DTTB(27X9,H,DH,W,WG,RO, VDPT) 
WRITE(55,400)WG 

400 FORMAT(2X,'OPB FUEL FLOWRATE - W(8805) >\F9.2,' LB/S') 

CALLFLOWOR(296,DTTB(28),DTTB(29),DTTB(30),9,H,DH,W,WG,RO,VDPT) 
WRITE(55,410) W 

410 FORMAT(2X,TPB LOX FLOWRATE- W(88 10) >\F9.2,' LB/S') 

C 

C THE TIME STEP - H- FOR THE PARTIAL DERIVATIVES IS 
H=1.001 

C CALCULATE VENTURI FLOW FOR VENTURI 20 AND THE 
C ASSOCIATED VENTURI UNCERTAINTY, THE UNCERTAINTY DATA IS 
BP1=0.01*DTTB(16) 

BT=0.01*DTTB(18) 

BDP=0.01*DTTB(17) 

BCD=0.01992 
BALPHA=4. 1 7E-07 
BD1=0.0005 
BD2=0.0005 
PW=0.0 

CALL FLOWOR(20,DTTB(1),DTTB(2),DTTB(3),9,ILDH,W,WG,RO,VDPT) 

C WRTTE(55,720) WG 

C720 FORMAT(2X,LPFT INLET FLOWRATE - P(8801) >'JF9.2,' LB/S 1 ) 

C CALCUALTE THE PARTTIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(1)*H 

CALL FLOWOR(20,HP1,DTTB(2)J>TTB(3X9,ILDH,W,HW,RO,VDPT) 
DWDPl=(HW-WGy(HPl-DTTB(l)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(2)*H 

CALL FLOWOR(20DTTB(1)4IDPDTTB(3)^ADH,W > HW,RO,VDPT) 
DWDDP=(HW-WGy(HDP-DTTB(2)) 

C PARITIAL DERIVATIVE FOR T1 
HT1=DTTB(3)*H 

CALL FLOWOR(20,DTTB(1)J)TTB(2),HT1,9,HJDH,W,HW,RO,VDPT) 
DWDT=(HW-WGy(HTl-DTTB(3)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(20J>TTB(iyDTTB(2XDTTB(3),l JIJDH,W3W4IO,VDPT) 
DWDDl=(HW-WGyDH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(20J)TTB(l)J)TTB(2yDTTB(3X2JtLDH,W3W r RO,VDPT) 
DWDD2=(HW-WGyDH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(20,DTTB( 1 )J)TTB(2)P)TTB(3),3,H,DH,W,HW,RO,VDPT) 
DWDCD=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(2045TrB(iyDTTB(2yDTTB(3X5PLDH,W3Wp.O,VDPT) 
DWDALP=(HW-WGyDH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=<BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD»DWDCD)**2 

BWALPHA=(BALPHA‘DWDALP)**2 
BWD1=(BD1*DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD 1 BD2=DWDD1 *DWDD2*BD 1 *BD2* 1 .00 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
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BTBP 1 =DWDT*DWDP1 *(BT*BP 1 )*0.00 
BDPBP1=DWDDP*DWDP1*(BDP*BP1)*0.00 
C VENTURI 20 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD1BD2+BWD1+BWD2 

&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT® 1 00. *UW/WG 
WW(1)=WG 
UNW(1)=UW 
PER( 1 )=PERCENT 
C VENTURI 20 OUTPUT 

WRITE (55,*) • PID 8801 LPFT INLET 
WRITE (55, »)' VENTURI <20>' 

WRITE (55,721) DTTB(1),DTTB(2),DTTB(3) 

WRITE (55,722) BP13DP3T 
WRITE (55,723) BCD,BALPHA,BD 1 ,BD2 
WRITE (55,724) PW 

WRITE (55,725) DWDP1 JDWDDPDWDT 
WRITE (55,726) DWDCDJDWDALPJ3WDD 1.DWDD2 
WRITE (55,727) BWP3WDP.BWT 
WRITE (55,728) BWCD3WALPHA3WD13WD2 
WRITE (55,729) WG.UW.PERCENT 

721 FORMAT(3X,Tl',4XF8.2,5X,'DP3X38.2,5X,T3X38.2) 

722 F0RMAT(3X,EP1'3X38.4,5X,EDP3X,F8.4,4X,'BT1',1X,F8.4) 

723 FORMAT(3X,’BCD' r 3X38.5,5X,EALP3X38-3 r 3X,BDl , ^X38.5, 

&6X,ED2’^X38.5) 

724 FORMAT(3X, I PW,4X,F8.4) 

725 F0RMAT(2X,T>WDP1',3X39.4,2X,'DWDDP,3X39.4,2X,'DWDT,3X39.3) 

726 FORMAT(2X,’DWDCD , r 3X39.4^X,T>WALP^XE9.3,2X,T)WDDl^X,E9.3, 
&2X,DWDD2 , ,2X39.4) 

727 FORMAT(2X,’BWP,5X39.43X'BWDP^XE9.4 r 3X,'BWDT^X39.4) 

728 FORMAT(2X,EWCD',4XE9.43X,BWALP3X39.4^X,'BWDl , ^X,E9,4, 
&4X,'BWD2'^X39.4) 

729 FORMAT(2X,’ W ^3X39.4,4X,'UW’,5X,F9.5,2X, , % UW\2X,F8.4) 

WRITE (55,*) '***»* ********* ***** •*****»• *»***• 

WRITE (55,*) 

C 

C CALCUALATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 397 
CALL FLOWOR(3973TTB(4)3TTB(5XDTTB(6)^ILDH,W,WG30,VDPT) 

C WRITE(55,730)WG 

C730 FORMAT(2X,'CCV INLET FLOWRATE - P(8818) >'39.2,' LB/S') 

BP1=0.01*DTTB(4) 

BT=0.0 1 *DTTB(6) 

BDP=O.0 1 *DTTB(5) 

BCD=0.03063 
BALPHA=5.48E-07 
BD1=0.0005 
BD2 =0.0005 
PW=0.0 

C CALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1 =DTTB(4)*H 

CALL FLOWOR(3973P13TTB(5)JDTTB(6X9ADH,W3W30,VDPT) 
DWDPI=(HW-WG)/(HP1-DTTB(4)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(5)*H 

CALL FLOWOR(397iyrrB(4)3DPJ)TTB(6)^ ADH,W 3 W 30 ,VDPT) 
DWDDP=(HW-WG)/(HDP-DTTB(5)) 

C PARITIAL DERIVATIVE FOR T1 
HT1=DTTB(6)*H 

CALL FL0W0R(3974)TTB(4)T)TTB(5XHTUTLDH,WTIW3O,VDPT) 
DWDT=(HW-WGy(HTl-DTTB(6)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(3973TTB(4)T>TTB(5XDTTB(6),l,HT)H,W,HW30,VDPT) 
DWDD1=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(3973TTB(4),DTTB(5)JDTTB(6)^4LDH,WTrW30,VDPT) 
DWDD2=(HW-WGyDH 
C PARTIAL DERIVATIVE FOR CD 

CALL FL0W0R(397T>TTB(4)JDTTB(5XDTTB(6P33H,W3W30,VDPT) 
DWDCD=(HW-WGyDH 
C PARTIAL DERIVATIVE FOR ALPHA 
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CALL FU3WOR(397,DTTB(4)J)TTB(5)J)TTB(6),4JiDH,W3W^O,VDPT) 

DWD ALP=(HW -WG)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP I *DWDP 1 )**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWD1=(BD1*DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1*DWDD2*BD1*BD2*1.0 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.00 
BDPBP1=DWDDP*DWDP1*(BDP*BP1)*0.00 
C VENTURI 397 UNCERTAINTY IS 

UW=SQRT (BWP+BWT +BWDP+BWCD+B W ALPHA+2. *BD 1 BD2+B WD 1 +B WD2 
&+2*BTBDP+2*BDPBPl+PW**2) 

PERCENT=100.*UW/WG 
WW(2)=WG 
UNW(2>=UW 
PER(2)=PERCENT 
C VENTURI 397 OUTPUT 

WRITE (55,*) ' PID 8818 CCV INLET* 

WRITE (55,*)* VENTURI <397>* 

WRITE (55,731) DTTB(4),DTTB(5),DTTB(6) 

WRITE (55,732) BP1,BDP,BT 
WRITE (55,733) BCD3ALPHAJBD13D2 
WRITE (55,734) PW 

WRITE (55,735) DWDP1JDWDDP,DWDT 
WRITE (55,736) DWDCD,DWDALP,DWDD1 JDWDD2 
WRITE (55,737) BWPJBWDP3WT 
WRITE (55,738) BWCD,BWALPHAJBWD1,BWD2 
WRITE (55,739) WG,UW,PERCENT 

731 FORMAT(3X,*Pl’,4X38.2,5X, , DF3X38.2,5X, r r3X38.2) 

732 FORMAT(3X,'BPl , ^X38.4,5X,'BDF3X38-4,4X,*BTr,lX,F8.4) 

733 FORMAT(3X,'BCD , 3XJF8.5,5X, , BALF3X38.3^X, f BDl , ^X,F8.5, 

&6X,*BD2*^X38.5) 

734 FORMAT(3X,?W*,4X38.4) 

735 FORMAT(2X,DWDPl’3X39.4^X,DWDDF r 3XJE:9.4^X,DWDr3X39.3) 

736 FORMAT(2X,DWDCD^X39.4^X,DWALP , 3X39.3^X,DWDDl^X,E9.3, 
&2X,DWDD2\2XJE9.4) 

737 FORMAT(2X,’BWP,5X39.4^X,*BWDF3X39-43X, t BWDr^X39.4) 

738 FORMAT(2X,'BWCD’,4X39.43X, , BWALF3XJE9.4,2X, r BWDl , ^X,E9.4, 
&4X3WD2*,2X39.4) 

739 FORMAT(2X,* W*^XJF9A,4X, , UW’,5X39.5^X,*%UW*^X38.4) 

WRITE (55,*) ' ***** ********* ***** ******** ****** 

WRITE (55^*) 

C 

C CALCUALATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 3541 
CALL HjOWOR(3M13TTB(7)3riB(8)3TTB(9)^33H,W,WG30,VDPT) 

C WRTTE(55,740) WG 

C740 FORMAT(2X,NOZ CLNT FLOWRATE - P(8815) >*39.2,* LB/S*) 

BP1-0.01*DTTB(7) 

BT=0.01 *DTTB(9) 

BDP=0.01 *DTTB(8) 

BCD=0.02033 
BALPHA=5.48E-07 
BD 1=0.0005 
BD2=0.0005 
PW=0.00 

C CALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(7)*H 

CALL FL0W0R(35413P13TTB(8)3TTB(9X933H,W3WJR0,VDPT) 

DWDP 1 =(HW -WG)/(HP 1 -DTTB(7)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(8)*H 
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CALL FLOWOR(3541 J>TTB(7)^P,DTTB(9X9JtLDH,W3WJRO,VDPT) 
DWDDP=(HW-WG)/(HDP-DTTB(8)) 

C PARITIAL DERIVATIVE FOR T1 
HT1 =DTTB(9)*H 

CALLFLOWOR(3541,DTTB(7)^TTB(8),HT1,9,H,DH,\V,HW^O,VDPT) 

DWDT=(HW-WG)/(HT1-DTTB(9)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(3541^rrB{7) > DTTB(8)J)TrB(9),l 3,DH,W^W^O,VDPT) 

DWDD 1 =(HW -WG)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALLFLOWOR(3541,DTTB(7),DTTB(8),DTTB(9),2,H,DH,W,HW,RO,VDPT) 
DWDD2=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(3541JDTTB(7)^TTB(8)^TTB(9)3ADH,W,HW,RO,VDPT) 
DWDCD=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(35414)TTB(7)J)TTB(8),DTTB(9),4,HJ)H,W,HW^O,VDPT) 
DWDALP=(HW-WG)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=<BALPHA*DWDALP)**2 
BWD1=<BD1*DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1*DWDD2*BD1*BD2*L00 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.00 
BDPBP1 =DWDDP*DWDP1 *(BDP*BP1)*0.00 
C VENTURI 3541 UNCERTAINTY IS 

UW=SQRT(BWP+BWT +BWDP+BWCD+B WALPHA+2. *BD 1 BD2+BWD 1 +BWD2 
&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT=100.*UW/WG 
WW(3)=WG 
UNW(3)=UW 
PER(3>=PERCENT 
C VENTURI 3541 OUTPUT 
WRITE (55,*) 

WRITE (55,*) ’PID 8815 NOZ CLNT 
WRITE (55,*)' VENTURI <3541>’ 

WRITE (55,741) DTTB(4)JDTTB(5),DTTB(6) 

WRITE (55,742) BP1,BDP,BT 
WRITE (55,743) BCD3ALPHA3DI3D2 
WRITE (55,744) PW 

WRITE (55,745) DWDP1 ,DWDDP,DWDT 

WRITE (55,746) DWDCDJ)WDALPJ3WDD1J>WDD2 

WRITE (55,747) BWP3WDP3WT 

WRITE (55,748) BWCD3WALPHAJBWD1JBWD2 

WRITE (55,749) WG,UW,PERCENT 

741 FORMAT(3X,Tr,4XJF8^,5X,DP3X38.2^X,T3X38.2) 

742 FORMAT(3X, , BPl^X38.4,5X, , BDF3X38.4,4X, , BTr,lX38.4) 

743 FORMAT(3X, , BCD , 3X38.5,5X, , BALP3X383^X, , BDr^X38.5, 

&6X/BD2’,2X38.5) 

744 FORMATtfXTW^X^M) 

745 FORMAT(2X,DWDPl , 3XJE9.4^X,DWDDP^X39.4^X, , DWDr3X393) 

746 FORMAT(2X,'DWDCD3X39.4^X,'DWALP'3X39.3^X,T)WDDl , 3X39.3, 
&2 X, t DWDD2’^X39.4) 

747 FORMAT(2X, f BWP‘,5X39.43X, , BWDF3X39.43X , BWDT^X,E9A) 

748 FORMAT(2X, T BWCD,4X,E9.43X,*BWALF3X,E9.4^X,'BWDl‘^XJE9.4, 
&4X,*BWD2’^X39.4) 

749 FORMAT(2X,' W 'JXJ9A,4X; UW,5X39.5,2X, , % UW^F8.4) 

WRITE (55,*) ' ***** ********* ***** ******** ****** 

WRITE (55**) 

C 

c 

c CALCUALATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 3541 
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CALL FLOWOR(3 542 r DTTB( 1 0),DTTB( 1 1 ),DTTB( 12),9,H,DH,W,WG,RO, VDPT) 
C WRITE(55,750)WG 

C750 FORMAT(2X,'NOZ CLNT FLOWRATE - P(8816) >’,F9.2/ LB/S') 
BP1=0.01*DTTB(10) 

BT=0.01*DTTB(12) 

BDP=0.01*DTTB(11) 

BCD-0.02033 

BALPHA=5.48E-07 

BD1-0.0005 

BD2-0.0005 

PW=0.0 

C CALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(10)*H 

CALL FLOWOR(3542 > HPI,DTTB( I l),DTTB(12),9,H,DH,W,HW,RO,VDPT) 
DWDPI=<HW.WG)/(HP1-DTTB(10)) 

C PARTIAL DERIVATIVE FOR DP 
HDP-DTTB(11)*H 

CALL FLOWOR(3542J)rrB(10),HDP > DTTB(12),93,DH,WJHW,RO,VDPT) 
DWDDP=(HW-WG)/(HDP-DTTB(1 1)) 

C PARITIAL DERIVATIVE FOR T1 
HT1=DTTB(12)*H 

CALL FLOWOR(3542J)TTB(10)J>TTB(11)3T1^3,DH,W3W,RO,VDPT) 
DWDT=(HW-WGy(HTl-DTTB(12)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(3542 > DTTB(10XDTTB(1 !)4>TTB(12),l,liDH,W,HW,RO,VDPT) 
DWDD 1 =(HW -WG)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(3542 > DTTB(l 0),DTTB(1 l)^TTB(12)^ > HJDH,W3W,RO,VDPT) 
DWDD2=(HW -WG)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(3542 > DTTB(10)^TTBCI l)J^TIB(l2)^ 3 H^H 9 W 9 HW 9 RO,VDPT) 
DWDCD=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(3542 > DTTB(10XDTTB(1 l)J)TTB(12),4ADH,W3W^O,VDPT) 
DWDALP=(HW -WG)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP-{BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWDl=(BDl*DWDDl)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1*DWDD2*BD1*BD2*1.0 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.0 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.0 
BDPBP1=DWDDP*DWDP1*(BDP*BP1)*0.0 
C VENTURI 3542 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD1BD2+BWD1+BWD2 

&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT=100.*UW/WG 
WW(4>=WG 
UNW(4>=UW 
PER(4)=PERCENT 
C VENTURI 3542 OUTPUT 
WRITE (55,*) 

WRITE (55,*)’ PID 8816 NOZ CLNT 

WRITE (55,*)’ VENTURI <3542>' 

WRITE (55,75 1) DTTB( 1 0XDTTB( 1 1XDTTB(12) 

WRITE (55,752) BP1JBDP,BT 
WRITE (55,753) BCD3ALPHA3D13D2 
WRITE (55,754) PW 

WRITE (55,755) DWDP1,DWDDPJDWDT 

WRITE (55,756) DWDCDJ)WDALPJDWDD1JDWDD2 

WRITE (55,757) BWP3WDP,BWT 

WRITE (55,758) BWCD3WALPHA3WD13WD2 
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WRITE (55,759) WG,UW PERCENT 

75 1 FORMAT(3X,Tl*,4X,F8.2,5XT)F3X,F8.2,5X,T3X,F8.2) 

752 FORMAT(3X, f BPl^X,F8.4,5X,’BDF3X,F8.4,4X, , BTlMXF8.4) 

753 FORMAT(3X, , BCD , 3X,F8.5,5X, , BALF,3XE8.3,3X, , BDr,2X,F8.5, 

&6X, r BD2',2X,F8.5) 

754 FORMATCSXTW^XJS .4) 

755 FORMAT(2X, r DWDPl’3X39.4^X,T)WDDF^X39.4^X, , DWDT,3XE9.3) 

756 FORMAT(2X, f DWDCD , 3X^9.4^X, f DWALF3X^9.3,2X,T>WDDl , ,3XE9.3, 
&2X,*DWDD2',2X£9.4) 

757 FORMAT(2X,*BWP*,5X,E9.43X, , BWDF3X,E9.43X,TBWDr^X,E9.4) 

758 FORMAT(2X, r BWCD’,4XJE9-43X,'BWALF3X39.4,2X , BWDl^X,E9.4, 
&4X, f BWD2*^CE9.4) 

759 FORMAT(2X,' W , 3XJ9.4,4X f UW’,5XJF9.5^X, , % UWV2X,F8.4) 

WRITE (55,*) ' ***** ********* ***** ******** ****** 

WRITE (55^*) 

C 

C CALCUALATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 3543 
CALL FLOWOR(3543,DTTB( 1 3),DTTB( 14),DTTB( 1 5),9,ILDH,W, WG,RO, VDPT) 

C WRITE(55,760)WG 

C760 FORMAT(2X/NOZ CLNT FLOWRATE - P(8817) > f ,F9.2/ LB/S’) 

BPl=O.01*DTTB(13) 

BT=0.O1*DTTB(15) 

BDP=0.0I*DTTB(14) 

BCD=0.02033 

BALPHA=5.48E-07 

BD1-0.0005 

BD2=0.0005 

PW=0.00 

C CALCUALTE THE PARTTIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(13)*H 

CALL FLOWOR(3543^IP14)TTB(14) J)TTB(1 5)^,H£>H, W,HW,RO,VDPT) 
DWDP1=(HW-WG)/(HP1-DTTB(13)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB( 1 4)*H 

CALL FLOWOR(35434>TTB(13XHDP y DTTB(l 5X93»DH,W JIW^.O,VDPT) 

D WDDP=(HW-WGy(HDP-DTTB( 14)) 

C PARITIAL DERIVATIVE FOR T1 
HT1=DTTB(15)*H 

CALL FLOWOR(3543^TTB(13)J)TTB(14)3Tl^ADH,W,HW^O,VDPT) 

DWDT=(HW -WG)/(HT1-DTTB(1 5)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(3543J)TTB(13)J)TTB(14y)TTB(15),lADH,W^IW^O,VDPT) 

DWDD 1 =(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(3543^TTB(13)J)TTB(14XDTTB(15)AHJ)H,W4IW,RO,VDPT) 
DWDD2=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(3543JDTTB(13)J)TTB(14y3TTB(15P3J)H,W3WJlO,VDPT) 
DWDCD=(HW-WGyDH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL ITX)WOR(3543J)TTB(13)^TTB(14),DTTB(15),4ADH,W3W^O,VDPT) 

D WD ALP=(HW - WG)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=<BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA-(BALPHA*DWDALP)**2 
BWD1 =(BD1 *DWDD1 )**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD 1 BD2=DWDD 1 *DWDD2*BD 1 *BD2* 1 .0 
BTBDP= ; DWDT*DWDDP*(BT*BDP)*0.0 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.0 
BDPBP1=DWDDP*DWDP1 *(BDP*BP1)*0.0 
C VENTURI 3543 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD1BD2+BWD1+BWD2 
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&-K2. *BTBDP+2*B DPBP1 +PW**2) 

PERCENT=100*UW/WG 
WW(5)=WG 
UNW(5)=UW 
PER(5)=PERCENT 
C VENTURI 3543 OUTPUT 
WRITE (55,*) 

WRITE (55,*)’ PED 8817 NOZ CENT 
WRITE (55,*)’ VENTURI <3543>’ 

WRITE (55,761) DTTB(13)^TTB(14),DTTB(1 5) 

WRITE (55,762) BP1,BDP,BT 
WRITE (55,763) BCD,BALPHA,BD1,BD2 
WRITE (55,764) PW 

WRITE (55,765) DWDP1,DWDDPJ)WDT 

WRITE (55,766) DWDCD,DWDALP,DWDD1,DWDD2 

WRITE (55,767) BWPJBWDP,BWT 

WRITE (55,768) BWCD,BWALPHAJBWDl,BWD2 

WRITE (55,769) WG,UW,PERCENT 

761 FORMAT(3X 9 Tl , ,4X£$.2,5X 9 f DP , £X£&.2,5X,'T£X£8.2) 

762 FORMAT(3X,BPr3XJ r 8.4,5X,'BDP , 3X,F8.4,4X, , BTlMX^8.4) 

763 FORMAT(3X,BCD , 3X^8.5,5X,’BALP3X^8.33X,*BDl , ,2X^8.5, 
&6X,BD2*,2X^8.5) 

764 FORMAT(3X,TW,4X^8.4) 

765 FORMAT(2X,’DWDPl , 3X39.4,2X, , DWDDF3X,E9.4,2X, f DWDr3X,E9.3) 

766 FORMAT(2X, f DWDCD , 3XJE9.43X,T>WALP , 3X,E9.33X,’DWDDr3X,E9.3, 
&2X,T>WDD2\2X£9.4) 

767 FORMAT(2X, f BWP , ,5XE9.43X,’BWDP3X39.43X, , BWDr3X39.4) 

768 FORMAT(2X,BWCD^4X£9.4,3X,BWALF,3X£9.4,2X, , BWDr,2X£9.4, 
&4X,'BWD2 , 3X,E9.4) 

769 FORMAT(5X , W , 3XJ^.4,4X,TJW,5XJ?9.53X, ,0 /oUW , 3XJF8.4) 

WRITE (55,*) * ***** ********* ***** ******** ****** 

WRITE (55,*) 

C 

C CALCULATE VENTURI FLOW FOR VENTURI 139 AND THE 
C ASSOCIATED VENTURI UNCERTAINTY, THE UNCERTAINTY DATA IS 
BP1=0.01*DTTB(16) 

BT=0.0 1 *DTTB(1 8) 

BDP=0.01*DTTB(17) 

BCD=0.02066 
BALPHA=6.5 IE-07 
BD1=0.0005 
BD2=0.0005 
PW=0.0 

CALL FLOWOR(1394>TTB(16)4>TTB(17)jnTB(18X93^H,W,WGmVDPT) 
C WRITE(55,770) W 

C770 FORMAT(2X,LPOT INLET FLOWRATE - P(8802) >'£9.2,' LB/S*) 

C C ALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(16)*H 

CALL FLOWOR(1393Pl^TTB(17)J>TTB(18)3ADIiHW,WG^O,VDPT) 
DWDP1=(HW-W)/(HP1-DTTB(16)) 

C WRITE(55,*) l H , JH,DH*4>H, , HW’^IW, , DWDP l'J3WDPl 

C PARTIAL DERIVATIVE FOR DP 
C H= 1.0001 

HDP=DTTB(17)*H 

CALL FLOWOR(1394)TTB(16)JHDP4>TTB(18X93U>IiHW,WG^O,VDPT) 
DWDDP=(HW-W)/(HDP-DTTB(17)) 

C WRITE(55,*) TT3,T>H , 4)H, , HW’3W,DWDDF4)WDDP 
C PARITIAL DERIVATIVE FOR T1 
C H=1.0001 

HT1~DTTB( 1 8)*H 

CALL FLOWOR(1393TTB(16)J)TIB(17)3T133T>H3W,WG,RO,VDPT) 
DWDT=(HW-Wy(HTl-DTTB(18)) 

C WRITE(55,*) TT4iT>H’JDH, , HW3W,DWDr,DWDT 
C PARTIAL DERIVATIVE FOR D1 
C H=1.0001 

CALL FLOWOR(139^TIB(16)^TTB(17y)riB(18XlADILHW,WGJRO,VDPT) 
D WDD 1 =(HW - WyDH 

C WRITE(55 ,*) TT,H,DH\DH,TIW,HW, , DWDD 1 ',DWDD 1 

C PARTIAL DERIVATIVE FOR D2 
C H=1.0001 
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CALL FLOWOR( 1 39,DTTB( 1 6),DTTB( 17),DTTB( 1 8),2,H,DH,HW,WG,RO, VDPT) 
DWDD2=(HW -W)/DH 

CALL FLOWOR( 139,DTTB( 1 6),DTTB( 1 7),DTTB( 1 8),7,H,DH,HW, WG,RO,VDPT) 
DWDCD=(HW-WyDH 
C PARTIAL DERIVATIVE FOR ALPHA 
C H=L0001 

CALI. FLOWOR(139 > DTTB(16)T)TTB(17),DTTB( 1 8),4 ; H,DH,HW, WG.RO, VDPT) 
DWDALP=(HW-W)/DH 

C WRITE(55,*) 1 H’,H,DH,DH,'HW > HW,DWDALP' ) DWDALP 

C CALCULATE THE SQUARED TERMS FOR HTE PARITIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP1»DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWD1=(BD1*DWDDI)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1 BD2=DWDD1 *DWDD2*BD1 *BD2* 1 .0 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.0 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.0 
BDPBP 1 =DWDDP*DWDP 1 *(BDP*BP 1 )*0.0 
C VENTURI 139 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+BWD1+BWD2 

&+2.*BTBDP+2*BDPBPl+2.*BDlBD2+PW*PW) 

PERCENT=100.*UW/W 
WW(6>=W 
UNW(6)=UW 
PER(6)=PERCENT 
C VENTURI 139 OUTPUT 
WRITE (55,*) 

WRITE (55, *)’PID 8802 LPOT INLET 

WRITE (55,*)' VENTURI <139>’ 

WRITE (55,771) DTTB(16),DTTB(17XDTTB(18) 

WRITE (55,772) BP13DP3T 
WRITE (55,773) BCD3ALPHA3D 1 3D2 
WRITE (55,774) PW 

WRITE (55,775) DWDP1,DWDDP,DWDT 
WRITE (55,776) DWDCD.DWDALP.DWDD13WDD2 
WRITE (55,777) BWP3WDP3WT 
WRITE (55,778) BWCD3WALPHA3WD13WD2 
WRITE (55,779) W.UW3ERCENT 

771 FORMAT(3X,Tl',4X38.2,5X,'DP'3X38.2,5X,T3X38-2) 

772 FORMAT(3X,'BPl , 3X38.4,5X,'BDP3X38.4,4X,'BTl’, 1X38.4) 

773 FORMAT(3X,’BCDtfX38.5,5X,BALP , ,3X383,3X,BDl , ,2X38.5, 

&6X,BD2’3X38.5) 

774 FORMAT(3X,TW,4X38-4) 

775 FORMAT(2X,’DWDPl’3X39-43X,'DWDDP3X39-43X,DWDT3X39.3) 

776 FORMAT(2X,’DWDCD*^X39.44X,T)WALP , 3X393^X,T>WDDl , 3X39.3, 
&2X,T>WDD2\2X39.4) 

777 FORMAT(2X,BWP,5X39.4,3X,BWDP,3X39.4,3X,BWDr,2X39.4) 

778 FORMAT(2X,BWCD,4X39-43X,BWALP3X39.4^X,BWDr^X39.4, 
&4X,BWD2'3X39.4) 

779 FOEMAT(5X,’ W ’,3X39.4, 4X,’UW’,5X39.5,2X, , % UW\2X,F8.4) 

WRITE (55,*) ' ***** *»»**«»»* »»*»* *****»«» **»*»> 

WRITE (55 ’•) 

C 

C CALCULATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 426 
CALLFLOWOR(4263TTB(19)3TTB(20)3>TTB(21)333>H,W,WG3O,VDPT) 
BP1=0.01*DTTB(19) 

BT=0.01*DTTB(20) 

BDP=0.01*DTTB(21) 

BCD=0.01964 

BALPHA=5.49E-07 

BDI=0.0005 

BD2=0.0005 

PW=0.0 

C CALCUALTE THE PARTTIAL DERIVATIVES 
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C PARTIAL DERIVATIVE FOR PI IS 
HP1=DTTB(19)*H 

CALL FLOWOR(426,HP13TTB(20),DTTB(21),9,H,DH,HW,WG,RO, VDPT) 
DWDP1=(HW-W)/(HP1-DTTB(19)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(20)»H 

CALL FLOWOR(426,DTTB( 1 9),HDP,DTTB(2 1 ),9,ILDH > HW, WG,RO,VDPT) 
DWDDP=(HW-W)/(HDP-DTTB(20)) 

C PARITIAL DERIVATIVE FOR T1 
HT1=DTTB(21)*H 

CALL FLOWOR(426JOTTB( 1 9),DTTB(20),HTl,9 t H,DH,HW,WG,RO,VDPT) 
DWDT=(HW-Wy(HTl-DTTB(2 1)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(426JDTTB( 1 9)J)TTB(20)JDTTB(2 1 ), 1 ,ILDILHW,WG,RO, VDPT) 
DWDDI=(HW-W)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(426,DTTB( 1 9)J)TTB(20),DTTB(2 1 ),2,H,DH,HW, WG,RO, VDPT) 
DWDD2=(HW-W)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(4263TTB(19)3TTB(20),DTTB(21),7,H,DH,HW,WG,RO, VDPT) 
DWDCD=(HW-WyDH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(426 J5TTB(1 9XDTTB(20)JDTTB(2 1 ),4,H,DH3W,WG,RO,VDPT) 
DWDALPHA=(HW-WyDH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 

BWP=(BP1*DWDP1)**2 

BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALJ>HA=(BALPHA*DWDALPHA)**2 
BWD1=(BD1 *DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1*DWDD2*BD1*BD2*1.0 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.0 
BTBP 1 =DWDT*DWDP 1 *(BT*BP 1 )*0.00 
BTBP1=DWDDP*DWDP1*(BDP*BP1)*0.00 
C VENTURI 426 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD1BD2+BWD1+BWD2 

&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT=100.*UW/W 
WW(7)=W 
UNW(7)=UW 
PER(7)=PERCENT 
C VENTURI 426 OUTPUT 
WRITE (55,*) 

WRITE (55,*) 'PID 8819 HPOP DISCH 1 
WRITE (55,*) 'VENTURI <426>' 

WRITE (55,781) DTTB(19)J}TrB(20)J)TTB(21) 

WRITE (55,782) BP13DP3T 
WRITE (55,783) BCD3ALPHA3D13D2 
WRITE (55,784) PW 

WRITE (55,785) DWDPIJJWDDPJJWDT 
WRITE (55,786) DWDCD3WDALP3WDD13WDD2 
WRITE (55,787) BWP3WDP3WT 
WRITE (55,788) BWCD3WALPHA3WD13WD2 
WRITE (55,789) W,UW,PERCENT 

781 FORMAT(3X,Tl’,4X38.2,5X,DF3X38.2,5X,T^X38.2) 

782 FORMAT(3X,'BPr3X38.4,5X,'BDP’^X38.4,4X,'BTl',lX,F8.4) 

783 FORMAT(3X,’BCD'3XJ8.5,5X,’BALP'^X38.33X,'BDl’^X38.5, 
&6X,*BD2'^X38.5) 

784 FORMAT(3X,’PW\4X38.4) 

785 FORMAT(2X,*DWDPl\3X39A,2X, , DWDDP , ,3X39A,2X,'DWDT’,3X39.3) 

786 FORMAT(2X,'DWDCD’3X39.4,2X,DWALP3X393^X,'DWDDl , 3X39.3, 
&2X,T>WDD2’3X39.4) 

787 FORMAT(2X,’BWF,5X39.43X,’BWDP’3X39.4,3X, , BWDT,2X,E9.4) 

788 FORMAT(2X,’BWCD’,4X39.43X,BWALP3X39.43X, , BWDl , 3X39.4, 
&4X,'BWD2’3X39.4) 
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789 FORMAT(5X,’ W , ,3X,F9.4,4X,XJW’,5X,F9.5,2X,'% UW\2X,F 8.4) 
WRITE (55,*) ' ***** ********* ***** ******** ****** 

WRITE (55!*) 


CALCULATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 268 
CALLFLOWOR(268,DTTB(22)4)TTB(23),DTTB(24)^ADH,W,WG^O,VDPT) 

BP1 =0.0 1 *DTTB(22) 

BT=0.01*DTTB(24) 

BDP=0.01*DTTB(23) 

BCD=0.02066 
BALPHA=6.7 1 E-07 
BD 1=0.0005 
BD2 =0.0005 
PW=0.0 

C CALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR Pi IS 
HP 1 =DTTB(22)*H 

CALL HX)WOR(268JIP14>T33(23>J)TTB(24X9ADH^W,WG30,VDPT) 

DWDP 1 =(HW «W)/(HP 1 -DTTB(22)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(23)*H 

CALL FLOWOR(268J>TTB(22XHDP,DTrB(24),9,H^HJIW,WG^O, VDPT) 
DWDDP=(HW-W)/(HDP“DTTB(23» 

C PARITIAL DERIVATIVE FOR T1 
HT1 =DTTB(24)*H 

CALL FLOWOR(268J3TTB(22) t DTTB(23),HT1^3^H3W,WG > RO,VDPT) 
DWDT=(HW-W)/(HT 1 -DTTB(24)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(268J)TIB(22)43TTB(23)J)TTB(24),lADILHW,WG^O,VDPT) 
DWDD 1 =(HW-W)/DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(268J3TTB(22)J)TTB(23)J)TTB(24)^ 34>ILHW,WGJtO,VDPT) 

D WDD2=(HW -W)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(268JC)TTB(22yDTTB(23)J3TTB(24),73J>H^W,WG^O,VDPT) 
DWDCEKHW.WyDH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(268»DTTB(22)Jl)TTB( 23)J)TTB(24),4,HJ)ILHW,WG,RO,VDPT) 
DWDALP=(HW-WyDH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWD1=(BD1*DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1 *DWDD2*BD1 *BD2 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
BTBP 1 =D WDT*DWDP 1 *(BT*BP 1 )*0.00 
BTBP1=DWDDP*DWDP1*(BDP*BP1)*0.00 

UW=SQRT(B WP+B WT+B WDP+BWCD+B WALPHA+2.*BD 1 BD2+B WD 1 +B WD2 
&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT=100.*UW/W 
WW(8)=W 
UNW(8)=UW 
PER(8)=PERCENT 
C VENTURI 268 OUTPUT 
WRITE (55,*) 

WRITE (55,*) ’PID 8804 OPB LOX' 

WRITE (55,*)' VENTURI <268>* 

WRITE (55,791) DTTB(22y)TTB(23XDTTB(24) 

WRITE (55,792) BP1,BL>P,BT 
WRITE (55,793) BCD JB ALPHA3D 1 ,BD2 
WRITE (55,794) PW 

WRITE (55,795) DWDP1,DWDDP,DWDT 
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WRITE (55,796) DWDCD ,DWDALP,DWDD1,DWDD2 
WRITE (55,797) BWP,BWDP,BWT 
WRITE (55,798) BWCD,BWALPHA,BWD1,BWD2 
WRITE (55,799) W,UW,PERCENT 

79 1 FORMAT(3X, r Pr,4XF8.2,5X/DP',3X,F8,2,5X,T,3X,F8.2) 

792 FORMAT(3X,'BPl , 3X,F8.4,5X,'BDF3X,F8.4,4X,*BTl’,lXJ ? 8.4) 

793 FORMAT(3X, r BCD , 3X,F8.5,5X,’BALP*,3X,E9.33X, , BDr,2X,F8.5, 

&6X,’BD2',2X,F8.5) 

794 FORMAT (3X,TW 1 ,4XJ ? 8.4) 

795 FORMAT(2X, r DWDPl t 3X,E9.43X, , DWDDF3X,E9-43X,T)WDT,3X,E9.3) 

796 FORMAT(2X,T)WDCD , ,3X,E9.43X, , DWALF,3XJE9.3,2X, , DWDDr3X,E9.3, 
&2X,T)WDD2*3X,E9.4) 

797 FORMAT(2X, , BWP,5X39.43X, , BWDF3XJE9.43X, , BWDr,2X,E9.4) 

798 FORMAT(2X, f BWCD‘,4X^:9.4 > 3X, , BWALP*3X39.43X, f BWDr3X^:9.4, 
&4X*BWD2'3X,E9.4) 

799 FORMAT(5X,' W , 3X,F9.4,4X,TJW , ,5X,F9.53X,ro UW,2X,F8.4) 

WRITE (55,*) * ***** ********* ***** ******** ****** 

WRITE (55,*) 

C 

C 

C CALCULATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 271 
CALL FLOWOR(27 1 ,DTTB(25)JDTTB(26),DTTB(27),9,H,DH,W, WG,RO, VDPT) 

C WRITE(5 5,800) WG 

C800 FORMAT(2X,'OPB FUEL FLOWRATE - W(8805) >*,F9.2, f LB/S') 

BP1 =0.0 1 *DTTB(25) 

BT=0.0 1 *DTTB(27) 

BDP=0.0 1 *DTTB(26) 

BCD=0.01992 

BALPHA=3.23E-07 

BD1=0.0005 

BD2=0.0005 

PW=0.0 

C C ALCUALTE THE PARJTIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1 =DTTB(25)*H 

CALL FLOWOR(271^1J)TTB(26y)TTB(27X9ADH,W3W^O,VDPT) 
DWDPl=(HW-WGy(HPl-DTTB(25)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(26)*H 

CALL FLOWOR(271^TTB(25)4IDPJDTTB(27)^ ADH,WJHWJlO,VDPT) 
DWDDP=(HW-WG)/(HDP-DTTB( 1 4)) 

C PARITIAL DERIVATIVE FOR T1 
HT 1 =DTTB(27)*H 

CALL FIX>WOR(271JDTTB(25y)TTB(26)3n'133J5H,W3W^O,VDPT) 

D WDT =(HW- W G)/(HT 1 *DTTB(27)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(271,DTTB(25)4)TTB(26)JDTTB(27), l,HJ>H,W,HW,RO, VDPT) 
DWDD 1 =(HW-WG>DH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(271JDTTB(25XDTTB(26XDTTB(27)33>DH,W^IW^O,VDPT) 
DWDD2=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(271JDTTB(25y)TTB(26)J)TTB(27),73J3H,W3WJlO,VDPT) 
DWDCD=(HW-WG)/DH 
C PARTIAL DERIVATIVE FOR ALPHA 

CALL FLOWOR(271JDTTB(25XDTTB(26XDTTB(27),43*DH,WJHW^O,VDPT) 
DWDALP=(HW -WG)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTIALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWD1=(BD1*DWDD1)**2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2-DWDD1 *DWDD2*BD1 *BD2*I.O 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
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BTBP1 =DWDT*DWDP1 *(BT*BP1 )*0.00 
BTBP 1 =DWDDP*DWDP 1 *(BDP*BP1)*0.00 
C VENTURI 271 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD 1 BD2+BWD 1 +BWD2 
&+2.*BTBDP+2*BDPBPl+PW**2) 

PERCENT= 1 00. *UW/WG 
WW(9)=WG 
UNW(9)=UW 
PER(9)=PERCENT 
C VENTURI 271 OUTPUT 
WRITE (55,*) 

WRITE (55,*)* PID 8805 OPB FUEL* 

WRITE (55,*)' VENTURI <271>’ 

WRITE (55,801) DTTB(25),DTTB(26),DTTB(27) 

WRITE (55,802) BP13DPJBT 
WRITE (55,803) BCD,BALPHA3D1,BD2 
WRITE (55,804) PW 

WRITE (55,805) DWDP1JDWDDP,DWDT 
WRITE (55,806) DWDCD,DWDALP,DWDD1,DWDD2 
WRITE (55,807) BWP,BWDP,BWT 
WRITE (55,808) BWCD,BWALPHA,BWD1,BWD2 
WRITE (55,809) WG,UW,PERCENT 

80 1 FORMAT(3X,Tl , ,4XJ ? 8.2,5X, , DF3X^8.2,5X,T3X^8.2) 

802 FORMAT(3X,’BPl , 3X,F8.4,5X, , BDF^X4= , 8.4,4X, , BTr,lXJ8.4) 

803 FORMAT(3X,’BCD*^X^8.5,5X, , BALF3X,E9.3^X, , BDI , ^F8.5, 

&6X/BD2’,2XJ8.5) 

804 FORMAT(3X,TW\4XJ8.4) 

805 FORMAT(2X,T3WDPl^XJE9.4^X, f DWDDF^XJE9.4^X, , DWDT3XJE9.3) 

806 FORMAT(2X,T)WDCiy3X39.4^X, , DWALP , 3XJE9.3^X, f DWDDr3X^:9.3, 
&2X/DWDD2 , ,2X39.4) 

807 FORMAT(2X, f BWF,5XJE9.4^X, , BWDF3X39.43X f BWDT^X,E9.4) 

808 FORMAT(2X, f BWCD , ,4X39.43X, r BWALF^X^9.4^X, f BWDl , ^X,E9.4, 
&4X,*BWD2’,2XJE9.4) 

809 FORMAT(5X,' W UW\2XJ8.4) 

WRITE (55,*) * ***** ********* ***** ******** ****** 

WRITE (55,*) 

C 

C CALCULATE THE FLOWRATE AND FLOWRATE UNCERTAINTY FOR VENTURI 296 
CALL FLOWOR(2964)TTB(28)J)TTB(29)JDTTB(30X9ADH,W,WG^O,VDPT) 

BP1 =0.0 1 *DTTB(28) 

BT=0.01*DTTB(30) 

BDP=0.0 1 *DTTB(29) 

BCEN0.02052 
BALPHA=6.71E-07 
BD 1=0.0005 
BD2=0.0005 
PW=0.0 

C CALCUALTE THE PARITIAL DERIVATIVES 
C PARTIAL DERIVATIVE FOR PI IS 
HP1 =DTTB(28)*H 

CALL FLOWOR(2963P1^TTB(29)4)TTB(30X9ADH,HW,WG^O,VDPT) 

D WDP 1 =(HW -Wy(HP 1 -DTTB(28)) 

C PARTIAL DERIVATIVE FOR DP 
HDP=DTTB(29)*H 

CALL FLOWOR(296,DTTB(28)jn>P4>TTB(30)^ 4U>H3W,WGJ10,VDPT) 
DWDDP=(HW-Wy(HDP-DTTB(29)) 

C PARITIAL DERIVATIVE FOR T1 
HT1 =DTTB(30)*H 

CALL FLOWOR(2964)TTB(28y3TTB(29Xirri^ADH4IW,WG^O,VDPT) 
DWDT=(HW-W)/(HT1-DTTB(30)) 

C PARTIAL DERIVATIVE FOR D1 

CALL FLOWOR(29643TTB(28)J)TTB(29)J)TTB(30X1 3,DH,HW,WGJRO,VDPT) 
DWDDl=(HW-WyDH 
C PARTIAL DERIVATIVE FOR D2 

CALL FLOWOR(296 JHTB(28)JDTIB(29y>rrB(30X23^>IUIW,WG4lO,VDPT) 
DWDD2=(HW-WyDH 
C PARTIAL DERIVATIVE FOR CD 

CALL FLOWOR(296,DTrB(28)J>TTB(29)J)TTB(30),7ADH,HW,WG,RO,VDPT) 
DWDCD=(HW-W)/DH 
P PAPTTAT r>FRTVATTVF FOP AT PFTA 

CALLFLOWOR(296J>TTB(28XDTTB{29)J>TTB(30),4,H r DtyfW,WG^O,VDPT) 
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D WD ALP=(HW - W)/DH 

C CALCULATE THE SQUARED TERMS FOR THE PARTLALS TIMES THE 
C UNCERTAINTY IN THE VARIABLES 
BWP=(BP1*DWDP1)**2 
BWT=(BT*DWDT)**2 
BWDP=(BDP*DWDDP)**2 
BWCD=(BCD*DWDCD)**2 
BWALPHA=(BALPHA*DWDALP)**2 
BWD1=(BD1*DWDD1)«2 
BWD2=(BD2*DWDD2)**2 
C CALCULATE THE FLOW UNCERTAINTY 

C CALCULATE THE CORRELATED BIAS TERMS, SWAG FOR CORRELATION 
C COEFFICIENT USED AT THIS POINT 
BD1BD2=DWDD1*DWDD2*BD1*BD2*0.0 
BTBDP=DWDT*DWDDP*(BT*BDP)*0.00 
BTBP1=DWDT*DWDP1*(BT*BP1)*0.00 
BTBP1=DWDDP*DWDP1*(BDP*BPI)*0.00 
C VENTURI 296 UNCERTAINTY IS 

UW=SQRT(BWP+BWT+BWDP+BWCD+BWALPHA+2.*BD 1 BD2+BWD 1+BWD2 
&+2.*BTBDP+2*BDPBPl +PW**2) 

PERCENT=100.*UW/W 
WW(10)=W 
UNW(10)=UW 
PER( 1 0)=PERCENT 
C VENTURI 296 OUTPUT 
WRITE (55,*) 

WRITE (55,*) ' PID 8810 FPB LOX 
WRITE (55,*) * VENTURI <296>’ 

WRITE (55,81 1) DTTB(28),DTTB(29),DTTB(30) 

WRITE (55,812) BP13DPJBT 
WRITE (55,813) BCDJBALPHAJBD 1 JBD2 
WRITE (55,814) PW 

WRITE (55,815) DWDP1,DWDDP,DWDT 
WRITE (55,816) DWDCD,DWDALP,DWDD1,DWDD2 
WRITE (55,817) BWP,BWDP,BWT 
WRITE (55,818) BWCD,BWALPHA,BWD1,BWD2 
WRITE (55,819) W,UW,PERCENT 

811 FORMAT(3X,Tr,4XJ8.2^X,DP3XJ8.2,5X,T3X38.2) 

812 FORMAT(3X,’BPl , 3X38.4,5X,'BDP3X,F8.4,4X, , BTr,lX38.4) 

813 FORMAT(3X, f BCD , 3XJ8.5,5X, , BALP , 3X38.33X, , BDl , ^X,F8.5, 
&6X3D2',2XJ8.5) 

814 FORMAT(3X,TW*,4XJS.4) 

815 FORMAT(2X,’DWDPl , 3XJE9.4^X, , DWDDP3X39-4^X, f DWDT3X39-3) 

816 FORMAT(2X,DWDCD , ^X39.4^X,’DWALP3X,E9.3^X,DWDDr3X39.3, 
&2X,T>WDD2*,2X39.4) 

817 FORMAT(2X,’BWF,5X39.4^X, , BWDF^X39.43X, f BWDT^XJE9.4) 

818 FORMAT(2X,3WCD,4X39.43X, , BWALP^X39.4^X, , BWDl , ^X39.4, 
&4X/BWD2’,2X39.4) 

819 FORMAT(5X,’ W ^XT9.4,4X, W,5XJ9.5^X, , % UW\2XJ8.4) 

WRITE (55,*) * ***** ********* ***** ******** ****** 

WRITE (55,*) 

C 

C CALL FL0W0R(2473Tre(l)3TTB(2XDTrB(3)^33H,W,WG30,VDPT) 

C WRITE(55,911)W 

C91 1 FORMAT(2X/HTX INLET FLOWRATE - P(8803) >39.2/ LB/S') 

C 

WRITE (55/) 

WRITE(55,915) NSLC3TIME3TIME 
915 FORMAT(2X,^LNUM>’45^X,*SUCE START TIME >35.1, 

$ 2X/SUCE END TIME >35.1) 

WRITE (5 5,*) 

DO 925 1=1,10 

WRITE (55,920) NPID(I),WW(I),UNW(IXPER(I) 

920 FORMAT(2X,4(I4^X37.2^X36.2^F6.2)) 

925 CONTINUE 
WRITE (55,*) 

C 

C 

GOTO 1 
C 

999 STOP 
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END 


C 

C- 


-00310099 


C AUTHOR: J. BUTAS / J. LEAHY / R. NEUMEYER - DEC 1990 00320099 

C 00330099 

C DATED REVISIONS: NONE 00340099 

C 003 50099 

C PURPOSE: TO CALCULATE FLOW VENTURI FLOWS IN LBM/SEC 00360099 

C 00370099 

C FORM OF CALL: CALL FLOWOR 00380099 

C 00390099 

C INPUT: ID FLOW ORIFICE ID NUMBER 00400099 

C ICD SWITCH TO ITERATE ON CD=F(REYNOLDS) 00401099 

C PI UPSTREAM PRESSURE PSIA 00410099 

C T1 UPSTREAM TEMPERATURE DEGR 00420099 

C DP ORIFICE DELTA PRESSURE PSIA 00430099 

C 00440099 

C OUTPUT: WL FLOWRATE, UQ CALCULATION LBM/SEC 00450099 


WG 

RHO 


FLOWRATE, GAS CALCULATION LBM/SEC 00460099 


DENSITY 


LBM/FT**3 00470099 
00480099 

00490099 

00500099 

00510099 

00520099 

00530099 

00540099 


00550099 


C 
C 
C 

C METHOD: 

C 

C RESTRICTIONS: FORT77 
C 

C NOTES: 

C 

GASEOUS FLOW EQ REQUIREMENTS : 

00560099 

DP = MEASURED DELTA PRESSURE PSID 00570099 
PI - MEASURED INLET PRESSURE PSIA 00580099 
Tl = MEASURED INLET TEMPERATURE DEGR 00590099 
D1 =* VENTURI INLET DIAMETER IN 00600099 

D2 - VENTURI THROAT DIAMETER IN 00610099 

CDFV = DISCHARGE COEFFICIENT MEASURED VIA WATER/AIR 00620099 
TESTS AT RKDN ENGINEERING DEV LAB IN CANOGA PK 00630099 
Z - COMPRESSIBILITY = F(P1,T1) 00640099 

GAM = GAMMA - F(P1,T1) 00650099 

RHO = DENSITY =F(P1,T1) LBM/FT3 00660099 
00670099 

LIQUID FLOW EQ REQUIREMENTS : 00680099 

00690099 

DP = MEASURED DELTA PRESSURE PSID 00700099 
PI = MEASURED INLET PRESSURE PSIA 00710099 
Tl = MEASURED INLET TEMPERATURE DEGR 00720099 
Dl = VENTURI INLET DIAMETER IN 00730099 
D2 = VENTURI THROAT DIAMETER IN 00740099 

CDFV = DISCHARGE COEFFICIENT MEASURED VIA WATER/AIR 00750099 
TESTS AT RKDN ENGINEERING DEV LAB IN CANOGA PK 00760099 
CDFVT= THEORETICAL DISCH. COEFFICIENT FOR ID = 247. 00770099 
GAM = GAMMA =F(P1,T1) 00780099 

RHO = DENSITY =F(P1,T1) LBM/FT3 00790099 
00800099 

00810099 
00820099 


C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

c REVISIONS: 
c NONE 


-00830099 


SUBROUTINE FLOWOR(ID4>l^)P,Tl^VAR4LI)H,WL,WG^HO,VDPT) 00840199 
IMPLICIT REAL (A-H,0-Z) 00841 199 

C0MM0N/PR0PTY/KU43W>V3UIV,S^L,SV,CV,CVL,CVV,CP,CPUCPV,GAMMA, 

X GAMMAL>GAMMAV,C,CL>CVP,MU,MUL*MUV,K4CL>KV,SIGMA^XCESJyiXCL,EXCV 
REAL BETA,C,CD,CDFVA,CDFVC,CPH2,CVH2,CPHE,CVHE,FV,Kl,K2,MACHCOR 
REAL KM,KVDI^IUH2^1UHEJ > COOUPHEJlD,TE,TEQU^C, Y COOL* YHE 
REAL WH2 203,WHEH_203>IWH2^1WHE^WMIX,CPMIX,CVMIX4CMIX,C1,C2 
REAL MU,MUL^tfUV,K^Q^V^A2^4UMIXJtfWO2>fWN2,P203RAT,RHEJlCOOL 
C 

REAL*8 PFlXJID,TFLUn},ROFLUID,VFLUID,HFLUID,SFLUID,CVFLUID, 00841399 
• CPFLUID,WFLUID,EFLUID,THFLUro,V2FLUID,ERFLUID,TSFLUID 00841499 
REALM XMW 00890099 

CHARACTER*3 FLD 00870099 

CHARACTER*2 TYPE 00880099 

C 00910099 
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poo o oo pooo po 


DATA NAMOX/2HG2/ 

DATA NAMHY/2HH2/ 

C 00920099 

C CHECK TO SEE IF ZERO. IF SO, SET WL TO VERY 00930099 

SMALL NUMBER AND EXIT. 00940099 

00950099 

WRTTE(6,*) 'VENTURI #',ID,P1,DP,T1 
WRITE (55,*) ENTERING SUBROUTINE FLOWOR’ 

WRITE (55,*) 7D , ,ID,'NVAR',NVAR 
IF (l.EQ.l) RETURN 

00950099 

IF ((PI .LE. 0) .OR. (T1 .LE. 0) .OR (DP .EQ. 0)) THEN 00960099 
WL = -10002.0 00970099 

WRJTE(86,*) INSTRUMENTATION INADEQUATE TO CALC. FLOWRATE' 

GOTO 999 00990099 

END IF 01000099 

01010099 

SET A DUMMY VALUE FOR THE RETURNED STEP SIZE 
DH-9999. 

DH IS SET UNREASONABLY HIGH TO BE VISIBLE IN OUTPUT 
TYPE = 'RO* 01020099 

0 1030099 

VENTURI CONSTANTS AS A FUNCTION OF VENTURI NUMBER (IMRL ITEM NO.) 01040099 
01050099 


MACHCOR = 1.0 

IF (ID .EQ. 20) THEN 01130099 

C *** LPFTP TURBINE DRIVE DUCT (ALL COMMENTED VARIABLES IN THIS SECTION01 140099 
C REFER TO WATER FLOW CHARACTERISTICS. ALL CURRENT VARIABLES ARE FROM 01140099 
C AIR FLOW TEST DATA (LPFT VENTURI WAS RE-CAUBRATED WITH AIR DUE TO: 01 140099 
C E0215(A1666) AND LOW REYNOLDS NUMBER H20 CAL. JC LEAHY 3-4-92 01140099 


IOP1 = 2 
FLD = *H2’ 
D1 = 1.9920 


01150099 

01160099 

01170099 


IF (NVAR.EQ.1) THEN 
D1=D1*H 


DH-D1-D1/H 

ENDIF 

D2 = 1.5225 01180099 

IF (NVAR.EQ.2) THEN 
D2=D2*H 


DH-D2-D2/H 

ENDIF 

C CDFV = 0.99610 
CDFV= 1.0125 
IF (NVAR.EQ3) THEN 
CDFV=CDFV*H 
DHCDFV-CDFV/H 
ENDIF 
XK = 0.448 
RP = 1.948 
ALPHA- 0.00000417 
IF (NVAR.EQ.5) THEN 
ALPOLD= ALPHA 
ALPHA- ALPHA* 1.01 
DH-ALPHA-ALPOLD 
ENDIF 
TYPE = ’RO’ 

C AO « 0.99522 
AO -1.0213183 
C BO - 0.000626 

BO = 0.01503 
C REDT = 2.50E+06 
REDT = 5.00E-HJ6 
C AA = 1.80 
AA = 0.20 
TREF- 480.0 
MACHCOR - 1.00246 
C 

ELSEIF(ID .EQ. 139) THEN 
C *** LPOTP TURBINE DRIVE DUCT 
IOP1 = 1 


01190099 

01190099 


01200099 

01210099 

01220099 


01230099 

01240099 

01240099 

01250099 

01250099 

01260099 

01260099 

01270099 

01270099 

01280099 

01290099 

01300099 

01310099 

01320099 
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nn ^ xxi ' 01330099 

D1 = 2.2998 01340099 

IF (NVAR.EQ.1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDEF 

D2 = 1.5743 01350099 

EF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 
ENDIF 

CDFV= 1.03332 01360099 

IF (NVAR.EQ.7) THEN 
CDOU>CDFV 
CDFV=CDFV*H 
DH=CDFV-CDOLD 

C WRITE (55,*) TD’JD^AR’^TVAR 
C WRTTE(55,*) , CDFV’,CDFV;CDOLD , ,CDOLD,T)H , ,DH 
ENDIF 


XK =0.316 
RP =2.441 
ALPHA= 0.00000651 
IF (NVAR.EQ.4) THEN 
ALPOLD= ALPHA 
ALPHA= ALPHA* 1.01 
DH=ALPHA*ALPOLD 
ENDIF 
AO =1.04232 
BO = 0.01258 
REDT = 4.00E+05 
AA= 1.00 
TREF= 190.0 
C 

ELSEIF(ID .EQ. 247) THEN 
C *** HEX OXIDIZER SUPPLY 
IOP1 = 1 
FLD = K>2* 

D1 =0.5535 
IF (NVAR.EQ.1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 
D2 =0.1745 
IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 
ENDIF 

CDFV = 0.9850 
IF (NVAR.EQ3) THEN 
CDFV=CDFV*H 
DH=CDFV-CDFV/H 
ENDEF 

ALPHA= 0.00000569 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 
DH=ALPHA*(1.-(17H)) 
ENDIF 
TREF= 185.0 
C 

ELSEIF(ID .EQ. 268) THEN 
C ***OPB OXIDIZER SUPPLY 
IOP1 = 1 
FLD =02’ 

D1 =1.0972 
IF (NVAR.EQ.1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 
D2 =0.6551 
EF (NVAR.EQ.2) THEN 
D2=D2*H 


01370099 

01380099 

01390099 


01400099 

01410099 

01420099 

01430099 

01440099 

01450099 

01460099 

« NOTE: CDFV IS AN ESTIMATE ***** 01470099 
01480099 
01490099 
01500099 


01510099 


01520099 


01530099 


01540099 

01550099 

01560099 

01570099 

01580099 

01590099 

01600099 


01610099 
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DH=D2-D2/H 


ENDIF 

CDFV =1.03717 
IF (NVAR.EQ.7) THEN 
CDOLD=CDFV 
CDFV=CDFV*H 
DHCDFV-CDOLD 
ENDIF 

ALPHA= 0.00000671 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 
DH=ALPHA*(l.-(l./H)) 
ENDIF 
AO =1.04585 
BO = 0.009103 
REDT = 6.00E+05 
AA = 1.10 
TREF = 200.0 
C 

ELSEIF(ID .EQ. 271) THEN 
C *** OPB FUEL SUPPLY 
IOP1 = 2 
FLD = H2' 

D1 =2.0060 
IF (NVAR.EQ.1) THEN 
D1=D 1*H 
DH=D1-D1/H 
ENDIF 
D2 =1.4931 
IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 


01620099 


01630099 


01640099 
01 650099 
01660099 
01670099 
01680099 
01690099 

01700099 

01710099 

01720099 

01730099 

01740099 


01750099 


ENDIF 

CDFV= 1.02162 
IF (NVAR.EQ.7) THEN 
CDOLD=CDFV 
CDFV=CDFV*H 
DH=CDFV-CDOLD 
ENDIF 
XK =0.276 
RP =1.233 
ALPHA= 0.00000323 
IF (NVAR.EQ.4) THEN 
ALPOLD= ALPHA 
ALPHA= ALPHA* 1 .0 1 
DH=ALPHA-ALPOLD 
ENDIF 
AO =1.02163 
BO = 0.003231 
REDT = 7.00E+05 
AA= 1.50 
TREF = 280.0 


01760099 


01770099 

01780099 

01790099 


01800099 

01810099 

01820099 

01830099 

01840099 

01850099 


ELSEIF(ID .EQ. 275) THEN 
C *** OPB ASI OXIDIZER SUPPLY 
IOP1 = 1 
FLD =*02' 

D1 =03030 
D2 =0.0818 
CDFV= 1.05226 
ALPHA= 0.00000549 
AO =1.05499 
BO = 0.02565 
REDT = 8.00E+05 
AA= 1.75 
TREF = 280*0 


01860099 

01870099 

01880099 

01890099 

01900099 

01910099 

01920099 

01930099 

01940099 

01950099 

01960099 

01970099 

01980099 

01990099 


ELSEIF(ID.EQ. 274) THEN 02000099 

C *** OPB ASI IGNITOR OXIDIZER SUPPLY 02010099 

IOP1 = 1 02020099 

FLD =*02' 02030099 
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D1 =0.3030 
D2 =0.0561 
CDFV = 0.99143 
ALPHA= 0.00000549 
AO = 0.99656 
BO = 0.007346 
REDT = 5.00E+05 
AA = 0.20 
TREF = 200.0 
C 

ELSEIF(ID .EQ, 277) THEN 
C ***OPBASI FUEL SUPPLY 
IOP1 = 2 
FLD =712’ 

D1 =0.4000 
D2 =0.1870 
CDFV = 0.99061 
ALPHA= 0.00000536 
AO = 0.99086 
BO = 0.002473 
REDT = 1.00E+06 
AA= 1.00 
TREF= 90.0 
C 

ELSEIF(ID .EQ. 296) THEN 
C ***FPB OXIDIZER SUPPLY 
IOP1 = 1 
FLD =’02’ 

D1 =2.0090 
IF (NVAREQ.l) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 
D2 =1.1090 
IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 
ENDIF 

CDFV= 1.02573 
IF (NVAREQ.7) THEN 
CDOLD=CDFV 
CDFV=CDFV*H 
DH=CDFV-CDOLD 
ENDIF 

ALPHA= 0.00000671 
IF (NVAREQ.4) THEN 
ALPOLD= ALPHA 
ALPHA= ALPHA* 1.01 
DH=ALPHA-ALPOLD 
ENDIF 
AO =1.02633 
BO = 0.0006843 
REDT = 1.80E+06 
AA= 1.50 
TREF = 200.0 
C 

ELSEIF(ID .EQ. 303) THEN 
C *** FPB ASI OXIDIZER SUPPLY 
IOP1 = 1 
FLD =02* 

D1 =03030 
D2 = 0.0928 
CDFV = 1.08156 
ALPHA= 0.00000549 
AO =1.08269 
BO = 0.02121 
REDT = 6.50E+05 
AA = 2.00 
TREF = 280.0 
C 

ELSEIF(ID .EQ. 302) THEN 


02040099 

02050099 

02060099 

02070099 

02080099 

02090099 

02100099 

02110099 

02120099 

02130099 

02140099 

02150099 

02160099 

02170099 

02180099 

02190099 

02200099 

02210099 

02220099 

02230099 

02240099 

02250099 

02260099 

02270099 

02280099 

02290099 

02300099 

02310099 

02320099 


02330099 


02340099 


02350099 


02360099 

02370099 

02380099 

02390099 

02400099 

02410099 

02420099 

02430099 

02440099 
02450099 
” 02460099 ~ 

02470099 

02480099 

02490099 

02500099 

02510099 

02520099 

02530099 

02540099 

02550099 

02560099 
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02570099 


* FPB ASI IGNITER OXIDIZER SUPPLY 


IOP1 = 1 02580099 

FLD =*02* 02590099 

D1 - 0.3030 02600099 

D2 =0.0502 02610099 

CDFV = . 98288 02620099 

ALPHA= 0.00000549 02630099 

AO = 0.98296 02640099 

BO = 0.00161 02650099 

REDT = 4.00E+05 02660099 

AA= 1.50 02670099 

TREF = 200.0 02680099 

02690099 

ELSEIF(ID .EQ. 305) THEN 02700099 

*** FPB ASI FUEL SUPPLY 02710099 

IOP1 = 2 02720099 

FLD =112' 02730099 

D1 =0.4000 02740099 

D2 =0.1811 02750099 

CDFV = 0.96581 02760099 

ALPHA= 0.00000536 02770099 

AO = 0.96835 02780099 

BO = 0.008258 02790099 

REDT= 1.00E+06 02800099 

AA = 0.50 02810099 

TREF = 90.0 02820099 

02830099 

ELSEIF(ID .EQ. 3541) THEN 02840099 

*** FLIGHT NOZZLE COOLANT SUPPLY NO. 1 


02890099 


02900099 


02910099 

02920099 

02930099 


IOP1 = 2 02860099 

FLD = *H2* 02870099 

D1 = 1.5241 02880099 

IF (NVAR.EQ.1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 

D2 = 0.9726 02890099 

IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 
ENDIF 

CDFV = 1.01644 02900099 

IF (NVAR.EQ.3) THEN 
CDFV=CDFV*H 
DH=CDFV-CDFV/H 
ENDIF 

XK =1.123 02910099 

RP = 1.656 02920099 

ALPHA= 0.00000548 02930099 

IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 
DH= ALPHA*( 1 .*( 1 7H)) 

ENDIF 

AO =1.01645 02940099 

BO = 0.002959 02950099 

REDT = 5.00E+05 02960099 

AA= 1.50 02970099 

TREF = 90.0 02980099 

C 02990099 

ELSEIF(ID .EQ. 3542) THEN 03000099 

C *** FLIGHT NOZZLE COOLANT SUPPLY NO. 2 
IOP1 = 2 03020099 

FLD=H2' 03030099 

D1 = 1.5250 03040099 

IF (NVAR.EQ.1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 

D2 = 0.9760 03050099 

IF (NVAR.EQ.2) THEN 
D2=D2*H 


02850099 


03010099 
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DH=D2-D2/H 


ENDIF 

CDFV = 0.99590 
IF (NVAR.EQ.3) THEN 
CDFV=CDFV*H 
DH=CDFV-CDFV/H 
ENDIF 
XK =1.141 
RP =1.702 
ALPHA= 0.00000548 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 
DH= ALPHA*( 1 .-( 17H)) 
ENDIF 
AO = 0.9959 
BO = 0.002620 
REDT = 5.00E+05 
AA= 1.50 
TREF= 90.0 


03060099 


03070099 

03080099 

03090099 


03100099 

03110099 

03120099 

03130099 

03140099 

03150099 


C 


C 


ELSEIF(ID .EQ. 3543) THEN 
♦** FUGHT NOZZLE COOLANT SUPPLY NO. 3 


03160099 


IOPl = 2 
FLD =W2' 

D1 =1.5205 
IF (NVAR-EQ.l)THEN 
D1=D1*H 


03180099 

03190099 

03200099 


DH=D1-D1/H 

ENDIF 

D2 = 0.9730 03210099 

IF (NVAR.EQ.2) THEN 
D2=D2*H 


DH=D2-D2/H 


ENDIF 

CDFV = 0.99918 
IF (NVAR.EQ.3) THEN 
CDFV=CDFV*H 
DH=CDFV-CDFV/H 
ENDIF 
XK =1.217 
RP =1.795 
ALPHA= 0.00000548 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 


03220099 


03230099 

03240099 

03250099 


DH=ALPHA*(1.-(17H)) 
ENDIF 
AO = 0.99919 
60 = 0.002990 
REDT = 5.00E+05 
AA= 1.50 
TREF= 90.0 


03260099 

03270099 

03280099 

03290099 

03300099 

03310099 


ELSEIF(ID .EQ. 397) THEN 
C •** CCV FUGHT NOZZLE COOLANT SUPPLY 


03320099 


IOPl =2 
FLD = H2' 

D1 =2.4382 
IF (NVAR.EQ. 1) THEN 
DI=DI*H 
DH=D1-D1/H 
ENDIF 
D2 =1.6939 
IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 


03340099 

03350099 

03360099 


03370099 


ENDIF 

CDFV = 1.02091 03380099 

IF (NVAREQ.3) THEN 
CDFV=CDFV*H 
DH=CDFV-CDFV/H 
ENDIF 


03170099 


03330099 
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XK =0.196 
RP =1.913 
ALPHA= 0.00000548 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 
DH=ALPHA*( 1 .-(1 ./H)) 
ENDIF 
AO - 1.02106 
BO = 0.01 152 
REDT = 5.00E+05 
AA = 1.00 
TREF = 90.0 
C 

ELSEIF(ID .EQ. 426) THEN 
C *** HPOTP PUMP DISCHARGE 
IOP1 = 1 
FLD =02’ 

D1 =3.9910 
IF (NVAR.EQ. 1) THEN 
D1=D1*H 
DH=D1-D1/H 
ENDIF 
D2 =3.2090 
IF (NVAR.EQ.2) THEN 
D2=D2*H 
DH=D2-D2/H 
ENDIF 

CDFV = 0.98195 
IF (NVAR.EQ.7) THEN 
CDOLD=CDFV 
CDFV=CDFV*H 
DH=CDFV-CDOLD 
ENDIF 

ALPHA= 0.00000549 
IF (NVAR.EQ.4) THEN 
ALPHA=ALPHA*H 


03390099 

03400099 

03410099 


03420099 

03430099 

03440099 

03450099 

03460099 

03470099 

03480099 

03490099 

03500099 

03510099 

03520099 


03530099 


03540099 


03550099 


DH=ALPHA*(1.-(17H)) 
ENDIF 
AO = 0.98091 
BO = 7340.0 
REDT = 1.00E+05 
AA= 5.00 
TREF = 200.0 


03560099 

03570099 

03580099 

03590099 

03600099 

03610099 

03620099 

03630099 

03640099 

03761199 


ELSE 
GOTO 999 
ENDIF 

IF (IOP1 .EQ. 1)THEN 
WRITE (6,*) *FLOW ORIFICE OUTPUT-UQUIiy 03761299 

WRITE (6,*) TO*, ID,T>l , JPl,Tr,Tl,T>PJDP 03761399 

WRITE (6,*) TLD^FLD 03761499 

WRITE (6 *) T) 1*,D 1 ,*D2’,D2,'CDFV\CDFV,'XIv ,XK 03763099 

WRITE (6,*) f RF4H > , , ALPHA , ALPHA,’AO^O, , BO , 30 03763199 

WRITE (6,*) REDTVREDT/AA',AA,TREP,TREF 03763199 

WRITE (6,*) NVAR , ^TVAR, , H , 3, , DH , t DH 
WRITE (6,*) 

WRITE (55,*) TDJD,NVAR , ^VAR, , H3, , DH’J)H 
ELSE 03764099 

WRITE (6 *) FLOW ORIFICE OUTPUT-GAS* 03770099 

WRITE (6 *) ID’, ID,*Pr,P 1 ,T 1 7n,*DP,DP 03770199 

WRITE (6,*) TLiyjLD 03761499 


WRITE (6,*) D1*,D1,*D2*,D2, , CDFV*,CDFV, , XK , ,XK 
WRITE (6,*) , RF,RP,'ALPHA , r ALPHA, , AO' r AO, , BO , ,BO 
WRITE (6,*) *REDT^EDT, , AA'AA,TREF,TREF 
WRITE (6,*) , NVAR’^rVAR, f H*A'DH*4)H 
WRITE (55,*) TD^D,'NVAR*^VAR, , H , 3, , DH , J)H 
WRITE (6,*) 

ENDIF 03781099 

03790099 

WRTTE(6,*) *BEGIN PROPERTY LOOKUP 
02770099 


03763099 

03763199 

03763199 
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C PROPERTY LOOK-UPS FOR GAMMA AND DENSITY 02780099 

C 02790099 

IF (FLD .EQ. ’H2 f ) THEN 02800099 

XMW = 2.016 02810099 

WRITE(6,*) *CALL - SETUP(NAMHY)’ 

CALL SETUP(NAMHY) 

WRITE(6,*) 'RETURN - SETUP(NAMHY)’ 

KU = 3 

C *** KS = 1 RETURNS DENSITY, KP = 1+2+4+8=15 RETURNS CP,CV AND MU... 
WRITE(6,*) 'CALL - GASP2 H2' 

CALL GASP2(1,15,T1,P1,DH2,HH2,3) 

WRITE(6,*) 'RETURN - GASP2 H2' 

RHO = DH2 
CPI - CP 
CV1 = CV 
VU = MU 

GAM = CP1/CV1 02920099 

ELSEIF(FLD .EQ. 02’) THEN 02940099 

XMW = 32.00 02950099 

WRITE(6,*) ’CALL - SETUP(NAMOX)' 

CALL SETUP(NAMOX) 

WRITE(6,*) 'RETURN - SETUP(NAMOX)' 

KU = 3 

C *** KS = 1 RETURNS DENSITY, KP = 1+2+4+8=15 RETURNS CP,CV AND MU- 
WRITE^,*) 'CALL -GASP2 02’ 

CALL GASP2(1, 1 5,T1,P1,D02302,3) 

WRITER,*) 'RETURN - GASP2 02' 

RHO = D02 
CPI = CP 
CV1=CV 
VU = MU 

GAM = CP1/CV1 02920099 

ENDIF 03070099 

WRITE(6,*)'END PROPERTY LOOKUP 
C 03080099 

WRITE(6*) 'BEGIN FLOW CALCULATIONS’ 

C 03090099 

C FLOW CALCULATIONS 03100099 

C 03110099 

C •** B IS F/M BETA RATIO (DIAMETER RATIO) WHERE D1 IS INLET DIA. & *** 03120099 
C *** D2 IS THE THROAT DIAMETER *** 03130099 


C *•* TE IS THE THERMAL EXPANSION FACTOR *** 
C •** CW IS THE CALIBRATION COEFFICIENT *** 


03140099 

03150099 


= D2/D1 031i 

= D2**2 *3.14159/4. 

= B**2 031' 

= B«4 0311 

= TREF - 528. 0: 

= 1.0 + (2 • ALPHA * (T1 - 528.)) 

= 1.0 + (2 * (ALPHA / T4) * (T1 - 528.)) 
= D2**2.*CDFV / (SQRT(1. - B4)) 


03160099 

03161099 

03170099 

03180099 

03180099 

03190099 

IS.)) 03190099 

03200099 


XP = 1 - DP/Pl 
P2 = PI - DP 
P20P1 =P2/P1 
G2 =27GAM 
GAM1 - GAM/(GAM-1) 
GAM2 = (GAM-1 yGAM 
GAM3 = (GAM+1)/GAM 


03201099 

03201199 

03202099 

03203099 

03204099 

03205099 

03205199 


YA =(XP**G2*GAMl)*((l-XP**GAM2y(l-XP))*((l-B4y(l-B4*XP**G2)) 03206099 

YA = SQRT(YA) 03207099 

CONST » 2. *32. 174/144. 03208099 

AREAC - 3.14159/4.00 03210099 

CDSAVE = CDFV 0321 1099 

WRTTE(6,*) END FLOW CALCULATIONS' 

C 03220099 

C LIQUID EQUATION OPTION 03230099 

C 03240099 

C *** THE CONSTANT 0.525019 IS A UNIT CONVERSION FACTOR BASED ON: 03250099 

C •** WL IN LB/S, DELTA P IN PSID, RHO IS LB/FT3, AND D2 IN INCHES *** 03260099 
C **• K = 3.14159/4*SQRT(2*GC/144) = 0.5250204 «• 03270099 

IF (IOP1 .EQ. 1) THEN 03280099 
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WRITE(6 ,*) 'BEGIN UQ EQUATION OPTION 1 
ENERGY = CONST*RHO*DP 03290099 

WL = MACHCOR * TE * CW * AREAC * SQRT(ENERGY) 03300099 

C ******** CURVE FIT FOR ID = 247 •***♦**♦*♦** 03320099 


C IB32 = 1 
C GO TO 301 
C300 CONTINUE 

CW = D2**2.*CDFV / (SQRT(I. - B4» 
ENERGY = CONST*RHO*DP 


03340099 

03341099 

03350099 

03350199 

03350299 


03350399 


WL « MACHCOR *TE*CW* AREAC * SQRT(ENERGY) 

C301 CONTINUE 03381099 

C RE = 48.0 * WL/ (3.14159 * D2 * VU) 03390099 

IF (ID .NE. 247) THEN 03392099 

CDFVT = AO - ABS(BO*(REDT/RE)**AA) 03391099 

ELSE 

CDFVT = .97253 + 3.3476E-08 * RE - 2.862E-14 * RE**2 + 03400099 

X 8.6492E-21 * RE**3 - 7.048E-28 * RE**4 03410099 

ENDIF 03411099 

CALL PSLP5 (CDFV,CDFV,CDFVT331,B32,0.0001,IB314B32^0300) 03420099 
IF (IB3 1 .EQ. 1) GO TO 300 03430099 

WLG = WL 03441099 

WRITE(6,*) ’END UQ EQUATION OPTION’ 

C 03450099 

C 03480099 

C 03490099 

C GAS EQUATION OPTION - IOP1 = 2 


03500099 


ELSE 


03510099 

03520099 


03740099 


03744399 


WRJTE(6,*) *BEGIN GAS EQUATION OPTIOFT 
ENERGY = CONST*RHO*DP 03521099 

WG = MACHCOR * YA*TE*CW* AREAC * SQRT(ENERGY) 

C IB42 = 1 03742099 

C GO TO 601 03743099 

C600 CONTINUE 03744099 

CW = D2**2.*CDFV / (SQRT(1. - B4)> 03744199 

ENERGY = CONST*RHO*DP 03744299 

WG - MACHCOR * YA *TE * CW * AREAC * SQRT(ENERGY) 

C60I CONTINUE 03748099 

C RE = 48.0 *WG/ (3.14159 *D2 * VU) 037 49099 

C CDFVT = AO -ABS(BO*(REDT/RE)**AA) 03749199 

C CALL PSLP5 (CDFV,CDFV,CDFVT341342,0.0001JB4UB42^0,600) 03749399 
C IF (IB41 .EQ. 1) GO TO 600 03749499 

WLG = WG 03749699 

WRITE(6,*) ’END GAS EQUATION OPTION’ 

ENDIF 03750099 

C *** VENTURI RESISTANCE CALCULATION *** 03751099 

IF (ID .EQ. 296) THEN 03751299 

VDPT = DP * (0.436 - 0.86 * B + 0.59 * B**2.) 03751399 

ELSE 03752199 

VDPT = DP * (0.218 - 0.42 * B + 0.38 * B**2.) 03752299 

ENDIF 03754099 

C 03760099 


IF (ID .EQ. 275) ID = 2741 
IF (ED .EQ. 274) ID = 2742 
IF (ID .EQ. 303) ID = 3021 
IF(ID.EQ. 302) ID = 3022 

WRITE (86,*) WENTURI ED NO:’JD 03761099 

IF (IOP1 JEQ. 1) THEN 03761199 

WRITE (86,*) TLOW ORIFICE OUTPUT-UQUID 03761299 

WRITE (86,*) TO 1 , ID,T1’J > 1,T1' > T1,’DP , 43P 03761399 

WRITE (86,*) f WL',WL,'RHO’,RHO,'GAMMA’,GAM,’VISCO',VU 03761499 

WRITE (86,*) *RE NO. , ^E,'CDFV’,CDFV,’CDFV THEO*,CDFVT 03763099 
WRITE (86,*) •GAMMA’,GAM,’VDPT,VDPT,’MACH COR’,MACHCOR 03763 199 
WRITE (86,*) , NVAR’JTVAR,’H’3, , DH’ J>H 
C WRITE (55,*) f NVAR , ,NVAR,’H’,H,’DH,DH 
ELSE 03764099 

WRITE (86,*) TLOW ORIFICE OUTPUT-GAS’ 03770099 

WRITE (86,*) TUy ID^l’JPLTlVTLDP^P 03770199 

WRITE (86,*) ’W G,WG,’RHO’^HO,’GAMMA , ,GAM,’VISCO’, VU 03771099 

WRITE (86,*) XF^P,’YA , ,YA,’XK’^K,'CW’,CW 03780099 

WRITE (86,*) ’RE NO.’J1E,'CDFV’,CDFV,ODFV THEO',CDFVT 03780199 
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03780299 


WRITE (86,*) 'VDPT.VDPT.'MACH COR'MACHCOR 
WRITE (86,*) NVAR\NVAR,W,H,DH\DH 
C WRITE (55,*) 'NVAR’.NVAR.'H'ADH’JJH 

ENDIF 03781099 

03790099 

WRITE (55,*) "LEAVING SUBROUTINE FLOWOR' 03802099 

999 RETURN 03810099 

END 03820099 

{THE THERMOPHYSICAL PROPERTY SUBROUTINE FOLLOWS 
AT THIS POINT IN THE PROGRAM, HOWEVER 
IT IS NOT INCLUDED IN THIS PRINTOUT} 


99 


